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Introduction: 


Background:  Estrogen  is  essential  for  normal  growth  and  differentiation  of  the  mammary  gland.  It 
also  supports  growth  of  70%  of  primary  breast  cancers.  The  biological  effects  of  estrogen  are 
mediated  by  two  estrogen  receptors,  ERa  and  ER(3.  Both  are  transcription  factors  regulating  target 
gene  expression  by  formation  of  dimers  at  the  genomic  sites.  While  ERa  is  the  prominent  target  in 
breast  cancer  therapy,  and  its  expression  strongly  predicts  response  to  endocrine  therapy,  ER|3  is 
neither  a  prognostic  marker  nor  a  therapeutic  target.  Expression  of  ERp  partially  overlaps  with  that  of 
ERa  in  breast  cancer  cells.  Among  ERa-positive  breast  cancer,  60%  of  cells  co-express  ERp.  ERp 
also  uniquely  expresses  in  50-80%  of  triple-negative  breast  cancer  (TNBC),  which  do  not  express 
ERa,  progesterone  receptor  and  Her2.  TNBCs  have  a  worse  prognosis,  and  targeted  therapies  are 
desperately  needed  for  this  group  of  women.  Our  idea  is  that  novel  treatment  strategies  should  be 
developed  based  on  the  function  of  different  ER  dimers  in  context  of  breast  cancer  subtypes. 

Although  dimerization  of  ERs  is  a  prerequisite  to  transcriptional  activation,  the  biology  of  ER  homo- 
and  hetero-dimers  (i.e.  ERa/a,  ERp/p,  ERa/p)  in  physiological  and  pathological  processes  are  poorly 
understood.  In  order  to  probe  for  ER  dimer  function,  we  developed  Bioluminescence  Resonance 
Energy  Transfer  (BRET)  assays  to  monitor  ER  dimerization  in  living  cells  and  identified  compounds 
selectively  induce  ERa/a,  ERp/p  homodimers  and  ERa/p-heterodimers.  Dimer-selective  estrogen 
receptor  modulators  will  be  explored  for  breast  cancer  prevention  and  treatment. 

Objective/Hypothesis:  We  hypothesize  that  ER  homo-  and  hetero-dimers  are  targets  for  different 
subtypes  of  breast  cancer.  Three  hypotheses  were  proposed:  1 )  For  ERa-positive  breast  cancer 
which  failed  conventional  treatment,  we  propose  to  target  CARM1,  a  histone  arginine 
methyltransferase  which  is  an  emerging  epigenetic  target  for  differentiation  therapy;  2)  We  have 
found  that  activating  ERp  activity  (agonism)  provides  anti-cancer  effects.  Therefore,  we  hypothesize 
that  drugs  that  activate  ERp  will  be  effective  for  treating  TNBC;  3)  Using  ERa/p-heterodimer  selective, 
natural  estrogenic  compounds  as  tools,  we  reveal  that  ERa/p  is  anti-proliferative  in  cell-based 
models.  We  hypothesize  that  ERa/p  heterodimers  may  be  targeted  in  60%  of  ERa-positive  tumors 
which  co-express  ERp,  as  well  as  served  as  cancer  preventive  targets  for  women  who  take  natural 
products  enriched  in  ERa/p  heterodimer  inducing  activities.  Our  overall  objectives  are  to  target  ERa+ 
tumors  for  differentiation  in  patients  who  failed  conventional  treatment,  to  target  ERp  for  treatment  of 
TNBC,  and  to  explore  the  biology  of  ERa/p,  the  most  mysterious  dimers  in  cancer  prevention  and 
treatment. 


PROGRESS  SINCE  LAST  REPORT: 

We  have  made  considerable  progress  towards  each  project  proposed  in  the  original  grant. 

AIM1:  To  explore  CARM1  as  epigenetic  target  for  differentiation  therapy  in  ERa-positive  breast 
cancer.  We  propose  to  develop  assays  for  identification  of  CARM1  methyltransferase  activators, 
which  may  be  developed  as  tumor  differentiation  drugs. 

We  have  demonstrated  that  overexpressing  CARM1  or  activating  CARM1  could  lead  to  differentiation  of 
ERa-positive  breast  cancer  cells.  In  collaboration  with  scientists  at  Life  Technology,  Madison,  Wl,  we 
have  developed  time-resolved  fluorescence  resonance  energy  transfer  (TR-FRET)  assays  that  can  be 
adapted  for  high  throughput  screening  to  identify  chemical  activator  of  CARM1  for  breast  cancer 
differentiation  therapy.  This  assay  is  based  on  methylation  of  poly(A)  binding  protein  1  (PABP1),  a 
known  CARM1  substrate,  by  cellular  CARM1 .  We  have  shown  that  PABP1  methylation  level  is 
proportional  to  CARM1  level  and  activity.  We  are  collaborating  with  Dr.  Hoffmann,  director  of  UW  Small 
Molecule  Screen  and  Synthesis  Facility,  to  optimize  the  assay  for  high  throughput  screening.  A 
manuscript  describing  the  assay  characterization  is  in  preparation  for  publication. 


a)  We  have  demonstrated  that  CARM1  is  a  unique  coactivator  of  ERa  and  over-expressing 
CARM1  induces  cell  differentiation  in  ERa-positive  breast  cancer  cells  (Al-Dhaheri,  M.,  et 
al.  (201 1 )  CARM1  Is  an  Important  Determinant  of  ERa-Dependent  Breast  Cancer  Cell 
Differentiation  and  Proliferation  in  Breast  Cancer  Cells.  Cancer  Res.,  71:  2118-2128).  At  the 
mechanism  level,  we  have  identified  PAF1  complex  as  “effector”  protein  to  be  recruited  to 
H3R17me2  mark  modified  by  CARM1  (Wu,  J.,  (2012),  Histone  H3R17me2a  mark  recruits 
human  PAF  complex  to  activate  transcription.  Proc.  Natl.  Acad.  Sci.  USA,  109:  5675-80). 

b)  TR-FRET  assay  have  been  established  for  identification  of  allosteric  activator  of  CARM1  for 
breast  cancer  differentiation  therapy.  A  manuscript  describing  the  assay  is  under  preparation 
(Zeng  H.,  Bi,  K.,  Hoffmann  M.,  Xu,  W.  Establishment  of  assays  to  monitor  cellular  activity  of  a 
histone  arginine  methyltransferase  CARM1). 

AIM2:  To  target  ER(3  for  treatment  of  triple-negative  breast  cancer.  We  will  evaluate  ER|3-selective 
agonists  as  imaging  tools  and  for  TNBC  treatment  and  identify  ERp  target  genes  as  ancillary 
biomarker  for  ongoing  TNBC  trial.  We  will  use  xenograft  mouse  model  to  develop  18Fluoro-ERp 
agonists  as  imaging  agents  for  PET  imaging  of  ERp+  human  tumors.  Genomic  approaches  will  be 
taken  to  identify  direct  ERp  target  genes.  These  effectors  will  be  tested  to  determine  if  they  are 
responsible  for  the  anti-cancer  effect  of  activating  ERp.  Also,  these  target  genes  will  be  tested  in 
tumor  samples  from  clinical  trial  to  determine  whether  estrogen  treatment  activated  ERp  in  patients, 
and  whether  TNBC  response  to  estrogen  correlate  with  ERp  expression  in  the  phase  II  clinical  trial. 

We  have  developed  isogenic  MDA-mb-468  cell  line  inducibly  express  either  ERa  or  ERp.  These  cell 
lines  have  been  validated  to  form  xenograft  tumors  in  nude  mice.  We  have  begun  to  assess  whether 
18Fluoro-17p-estradiol  (FES)  could  detect  ERp-positive  breast  tumors  in  vivo  and  thus  suitable  for  in  vivo 
PET  imaging  to  identify  ERp-expressing  breast  tumors.  We  have  optimized  conditions  to  use  a 
commercial  ERp  antibody  for  immunohistochemistry  (IHC).  Using  Tissue  Microarray  (TMA),  we 
demonstrated  that  -25%  of  triple  negative  breast  tumors  express  significant  amount  of  ERp.  Using 
MDA-mb-468-ERp  cell  line,  we  performed  RNA-sequencing  and  identified  novel  ERp  target  genes. 
Expression  of  these  target  genes  will  be  tested  in  tumor  samples  from  clinical  trial  to  determine  whether 
estrogen  treatment  activated  ERp  in  patients.  The  phase  II  clinical  trial  is  ongoing.  Ten  patients  have 
enrolled  in  the  clinical  trial.  One  patient  showed  regression  of  liver  metastasis  after  6-week  treatment  of 
high  dose  of  E2.  This  result  proves  the  principle  that  triple-negative  breast  cancer  patient  could  have 
response  to  endocrine  therapy. 

We  have  established  two  isogenic  cell  lines  for  screening  ERp-selective  ligands  (Shanle  E.,  Hawse  J., 
Xu,  W.  (2012),  Biochemical  Pharmacology,  82:1940-1949). 

AIM3:  To  delineate  the  biological  function  of  ERa/p  heterodimer  in  cultured  cells  and  in  vivo.  We  plan 
to  examine  whether  two  phytoestrogens  identified  to  be  ERa/p  heterodimer  selective  could  prevent 
estrogen-induced  breast  cancer  in  an  ACI  rat  model  through  collaboration  with  Jim  Shull’s  laboratory. 

We  will  also  identify  genes  regulated  by  ERa/p  heterodimers  in  normal  mammary  epithelial  cell  lines 
and  epithelial  cells  isolated  from  rat  mammary  gland. 

The  previously  established  BRET  assays  identified  two  naturally  occurring  phytoestrogens  that  induce 
ERa/p  heterodimer  formation  in  cell-based  assay.  We  hypothesize  that  compounds  inducing  ERa/p 
heterodimer  formation  in  vivo  could  have  preventive  and  therapeutic  effects  for  breast  cancer.  In 
collaboration  with  Dr.  Jim  Shull  in  my  department,  we  have  learned  to  implant  one  of  the  ERa/p 
heterodimer-selective  phytoestrogens,  cosmosiin,  as  pellet  in  the  back  of  ACI  female  rat.  We  have 
determined  pharmacokinetics  of  cosmosiin  in  ACI  rat  model  and  showed  that  serum  concentration  of 
cosmosiin  could  reach  pM  that  is  sufficient  to  induce  ERa/p  heterodimer  in  vivo  after  2-wks  cosmosiin 
implant.  We  have  begun  to  do  a  longer  time  pharmacokinetics  experiment  to  optimize  conditions  for 
breast  cancer  prevention  experiment.  Meanwhile,  we  are  modifying  BRET  assays  to  BRET3  for  in  vivo 


application  as  compounds  exhibiting  heterodimer  specificity  in  cell-based  assays  may  not  maintain 
heterodimer-specificity  in  vivo. 


Using  BRET  assays,  we  identified  two  naturally  occurring  phytoestrogens,  cosmosiin  and  angolensin, 
which  induce  ERa/p  heterodimer  formation  in  cells  (Powell,  E.  et.  Al.  (2012)  Plos  One,  7(2):e30993). 


Key  Research  Accomplishments: 

1 )  We  have  demonstrated  that  CARM1  is  a  unique  coactivator  of  ERa  and  over-expressing  CARM1 
induces  cell  differentiation  in  ERa-positive  breast  cancer  cells  (Al-Dhaheri,  M.,  et  al.  (201 1 ) 

CARM1  Is  an  Important  Determinant  of  ERa-Dependent  Breast  Cancer  Cell  Differentiation  and 
Proliferation  in  Breast  Cancer  Cells.  Cancer  Res.,  71 :  2118-2128).  At  the  mechanism  level,  we 
have  identified  PAF1  complex  as  “effector”  protein  to  be  recruited  to  H3R17me2  mark  modified  by 
CARM1  (Wu,  J.,  (2012),  Histone  H3R17me2a  mark  recruits  human  PAF  complex  to  activate 
transcription.  Proc.  Natl.  Acad.  Sci.  USA,  109:  5675-80). 

2)  TR-FRET  assay  have  been  established  for  identification  of  allosteric  activator  of  CARM1  for  breast 
cancer  differentiation  therapy.  A  manuscript  describing  the  assay  is  under  preparation  (Zeng  H.,  Bi, 
K.,  Hoffmann  M.,  Xu,  W.  Establishment  of  assays  to  monitor  cellular  activity  of  a  histone  arginine 
methyltransferase  CARM1). 

3)  We  have  established  two  isogenic  cell  lines  for  screening  ERp-selective  ligands  (Shanle  E.,  Hawse 
J.,  Xu,  W.  (2012),  Biochemical  Pharmacology,  82:1940-1949). 

4)  Using  BRET  assays,  we  identified  two  naturally  occurring  phytoestrogens,  cosmosiin  and 
angolensin,  which  induce  ERa/p  heterodimer  formation  in  cells  (Powell,  E.  et.  Al.  (2012)  Plos  One, 
7(2):e30993). 


Reportable  outcomes: 

See  Appendix  Materials  for  published  research.  We  have  developed  isogenic  ERa  or  ERp  expressing 
Hs578T  and  MDA-mb-468  cell  lines.  We  have  developed  cell-based  BRET  assay  for  high  throughput 
screening  of  naturally  occurring  estrogenic  compounds. 


Presentations: 

Department  of  Human  Oncology,  University  of  Wisconsin-Madison,  March  15,  2012.  Title:  Epigenetic 
control  of  breast  cancer. 

Oregon  Health  &  Science  University,  March  20,  2012.  Title:  Regulation  of  estrogen  signaling  in  breast 
cancer. 

World  Cancer  Congress,  Beijing,  May  18-21 , 2012.  Title:  Targeting  estrogen  receptor  coactivator  for 
breast  cancer  treatment. 

National  Cheng  Kung  University,  Taiwan,  June  7-9,  2012.  Title:  Transcriptional  control  of  estrogen 
signaling  in  breast  cancer. 


Recent  Funding  Applications: 


NIH,  submitted  6/2012:  Please  add  “Cell-based  HTS  for  Allosteric  Activators  of  the  CARM1  Arginine 
Methyltransferase”. 

We  propose  to  develop  and  optimize  a  cell-based  high  throughput  screening  assay  to  identify  small 
molecule  activators  of  CARM-1 ,  using  a  time-resolved  fluorescence  method  to  quantify  the  CARM1 
methylation  of  specific  arginines  on  poly  A  binding  protein  1  (PABP1 ).  We  hope  these  screens  will 
lead  to  more  facile  and  rigorous  approaches  to  understanding  the  roles  of  CARM1  activity  in  cancer 
cells  in  vitro  and  in  vivo,  as  well  as  provide  leads  for  anticancer  drugs  with  a  novel  mechanism  of 
action.  We  will  use  a  library  of  -50,000  compounds  available  at  the  UW-Madison  Small  Molecule 
Screening  &  Synthesis  Facility.  The  best  compound(s)  identified  from  these  assays  will  be  used  to 
interrogate  the  global  effects  of  CARM1  activation  on  re-programming  the  epigenome  through  RNA- 
seq  and  chromatin  immunoprecipitation  approaches.  Successful  completion  of  the  this  proposal  will 
provide  1 )  a  new  HTS  assay  for  an  important  molecular  target  in  epigenetic  regulation  that  can  be 
implemented  at  many  other  academic  or  NIH  screening  centers,  2)  an  assay  strategy  that  could  be 
adopted  for  HTS  with  other  protein  methylases,  and  3)  molecular  probes  for  testing  hypotheses  about 
CARM1  in  human  disease,  including  the  hypothesis  that  small  molecules  enhancing  CARM1  activity 
may  be  useful  as  “epigenetic”  drugs  for  ERa+  breast  cancers  that  fail  conventional  endocrine  therapy. 


Conclusions: 

We  have  demonstrated  that  CARM1  is  a  valid  epigenetic  target  for  breast  cancer  differentiation  therapy 
TR-FRET  assay  has  been  established  to  screen  for  small  molecule  activators  of  CARM1  to  be 
developed  as  epi-drugs.  To  our  knowledge,  assays  for  screening  “activators”  of  epigenetic  enzymes  do 
not  exist  and  epigenetic  therapy  for  breast  cancer  is  not  yet  widely  applicable.  We  have  provided  solid 
evidence  suggesting  that  ER(3  should  be  activated  in  triple-negative  breast  cancer,  both  in  vitro  and  in 
the  ongoing  clinical  trial.  Using  our  highly  optimized  isogenic  ERa  or  ERp  expressing  cell  lines,  animal 
models  have  been  developed  to  validate  FES  for  in  vivo  PET  imaging  to  screen  patients  whose  tumors 
express  ERp.  Finally,  two  ERa/p  heterodimer  inducing  phytoestrogens  have  been  identified.  One 
compound  cosmosiin  has  been  shown  to  reach  pM  serum  concentration  in  the  ACI  rat  model  in  a 
pharmacokinetics  experiment.  This  is  proof  of  principle  that  cosmosiin  might  have  chemo-preventative 
effects  in  vivo. 
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Abstract 

Breast  cancers  with  estrogen  receptor  a  (ERa)  expression  are  often  more  differentiated  histologically  than 
ERa-negative  tumors,  but  the  reasons  for  this  difference  are  poorly  understood.  One  possible  explanation  is  that 
transcriptional  cofactors  associated  with  ERa  determine  the  expression  of  genes  which  promote  a  more 
differentiated  phenotype.  In  this  study,  we  identify  one  such  cofactor  as  coactivator-associated  arginine 
methyltransferase  1  (CARM1),  a  unique  coactivator  of  ERa  that  can  simultaneously  block  cell  proliferation 
and  induce  differentiation  through  global  regulation  of  ERa-regulated  genes.  CARM1  was  evidenced  as  an  ERa 
coactivator  in  cell-based  assays,  gene  expression  microarrays,  and  mouse  xenograft  models.  In  human  breast 
tumors,  CARM1  expression  positively  correlated  with  ERa  levels  in  ER-positive  tumors  but  was  inversely 
correlated  with  tumor  grade.  Our  findings  suggest  that  coexpression  of  CARM1  and  ERa  may  provide  a  better 
biomarker  of  well-differentiated  breast  cancer.  Furthermore,  our  findings  define  an  important  functional  role  of 
this  histone  arginine  methyltransferase  in  reprogramming  ERa-regulated  cellular  processes,  implicating  CARM1 
as  a  putative  epigenetic  target  in  ER-positive  breast  cancers.  Cancer  Res;  71(6);  2118-28.  ©2011  AACR. 


Introduction 

In  normal  breast  tissue,  estrogen  receptor  a  (ERa)  regulates 
growth  and  development  of  the  mammary  gland  by  regulating 
the  balance  between  cell  proliferation  and  differentiation.  This 
balance  is  deregulated  in  cancer.  Enhanced  ERa  proliferative 
action  contributes  to  the  initiation  and  progression  of  breast 
cancer  (1)  by  promoting  cell-cycle  progression,  in  particular  S- 
phase  entry  (2, 3).  Microarray  analyses  using  breast  cancer  cell 
lines  have  revealed  that  a  majority  of  ERa  target  genes  are 
involved  in  metabolism  and  cell-cycle  regulation  (4-6).  ERa  is 
expressed  in  nearly  70%  of  breast  cancers.  Interestingly,  ERa- 
positive  tumors  are  more  histologically  well-differentiated  (7). 
ERa  decreases  in  high-grade  tumors  (8),  and  its  presence 
serves  as  a  hallmark  of  differentiation  and  predictor  of  low 
aggressiveness  and  favorable  disease-free  survival  (9,  10).  The 
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protective  effect  of  ERa  raises  the  possibility  that  ERa  func¬ 
tions  to  regulate  both  proliferation  and  differentiation  in 
breast  cancer  cells,  albeit  with  the  balance  tilted  toward 
proliferation.  Cell  proliferation  and  differentiation  are  2 
mutually  exclusive  processes.  Forced  differentiation  of  pri¬ 
mary  tumors  with  therapeutic  compounds  can  inhibit  pro¬ 
liferation  (11).  Differentiation  therapy  such  as  all -trans 
retinoic  acid  was  successfully  used  in  treating  acute  promye- 
locytic  leukemia.  However,  this  strategy  is  not  widely  applied 
to  breast  carcinomas  because  breast  tumors  are  more  hetero¬ 
geneous.  Moreover,  how  ERa  regulates  the  balance  of  pro¬ 
liferation  and  differentiation  is  not  well  understood. 

ERa  regulates  transcription  through  recruitment  of  multiple 
cofactors  (12).  Although  ERa  coactivators  share  the  common 
feature  of  activating  ERa  in  transcriptional  assays,  to  date,  no 
ERa  coactivator  has  been  reported  to  promote  differentiation 
in  breast  cancer  cells.  Coactivator-associated  arginine  methyl¬ 
transferase  1  (CARM1)  was  originally  identified  as  a  steroid 
receptor  coactivator  which  activates  transcription  of  ERa 
target  genes  (13, 14).  Furthermore,  loss  of  CARM1  in  the  mouse 
embryo  leads  to  abrogation  of  the  estrogen  response  and 
reduced  expression  of  some  ER  target  genes  (15),  highlighting 
the  significance  of  CARM1  in  ER-regulated  processes.  CARM1 
is  a  multifunctional  protein  engaged  in  a  variety  of  cellular 
processes  including  gene  expression  (16),  coupling  of  tran¬ 
scription  and  mRNA  processing  (17),  regulating  protein  stabi- 
lity  (18),  and  tissue  development  (15).  However,  the  function  of 
CARM1  in  regulating  cell  differentiation  or  proliferation  is 
contradictory  and  seems  to  be  context  dependent.  CARM1 
is  required  for  differentiation  of  adipocytes  (19),  myocytes  (20), 
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pulmonary  epithelial  cells  (21),  and  early  thymocyte  progenitor 
cells  (22).  In  contrast,  CARM1  was  implicated  in  cancer  cell 
proliferation  and  was  shown  to  regulate  the  expression  of  E2F 1 
and  cyclin  El,  factors  promoting  cell-cycle  progression  (16, 23). 
Thus,  functions  of  CARM1  in  ERa-dependent  breast  cancer 
require  further  elucidation. 

Here  we  report  an  extensive  study  of  the  biological  function 
of  CARM1  in  ERa-regulated  processes  in  breast  cancer  cells, 
using  both  gain-of-function  and  loss-of-function  approaches. 

Materials  and  Methods 

Cell  culture  maintenance  and  construction 

All  cell  lines  were  purchased  from  American  Type  Culture 
Collection  (ATCC)  and  used  within  6  months.  MCF7-Tet-on- 
shCARMl  was  generated  by  2  steps:  First,  we  synthesized 
shRNA-encoding  oligo  DNA  "CGGCGAGATCCAGCGGCAC" 
targeting  human  CARM1  and  cloned  the  sequence  into 
pLVTHM  plasmid  (24).  MCF7  cells  (ATCC  HTB-22)  were 
sequentially  infected  with  the  lenti-KRAB  and  pLVTHM- 
CARM1  shRNA  vectors,  followed  by  selection  of  single  clones 
by  Western  blotting.  MCF7-Tet-on-CARMl  was  generated  by 
cotransfecting  pTRE-tight  CARM1  plasmid  with  pBabe-pur- 
omycin  vectors  at  a  ratio  of  10:1  into  MCF7-Tet-on  cells 
(Clontech  cat  #630918,  generous  gift  of  Elaine  Alarid),  followed 
by  selection  with  puromycin.  MCF7-CARM1  has  been 
described  previously  (25).  MDA-MB-231  (ATCC  HTB-26)  or 
ZR-75  (ATCC  CRL-1500)  cells  were  infected  with  retrovirus 
derived  from  pLNCX  (Clontech),  pLNCX-CARMl,  and  pSIREN- 
Q-shCARMl  vectors  to  obtain  pooled  clones.  The  shRNA- 
encoding  oligo  DNA  "CAGCGTCCTCATCCAGTTC"  targeting 
human  CARM1  was  cloned  into  pSIREN-Q  (Clontech)  vector. 
All  cells  were  maintained  in  DMEM  containing  10%  FBS  and 
seeded  in  phenol  red-free  media  with  5%  of  stripped  FBS  for 
experiments. 

Immunohistochemistry  and  tissue  microarrays  of 
human  breast  tumors 

All  primary  invasive  breast  tumors  used  in  this  study  were 
obtained  from  the  Manitoba  Breast  Tumor  Bank  [MBTB; 
Cancer  Care  Manitoba  and  University  of  Manitoba  (http:// 
www.umanitoba.ca/institutes/manitoba_institute_cell_biol- 
ogy/MBTB/Index4.htm)].  The  MBTB  operates  with  approval 
from  the  Research  Ethics  Board  of  the  Faculty  of  Medicine, 
University  of  Manitoba.  Tissue  collection  and  selection  of 
samples  for  constructing  tissue  microarrays  (TMA)  have  been 
reported  before  (26,  27).  Although  450  cases  were  represented 
on  the  original  ER-positive  TMAs,  due  to  exhaustion  of  tumor 
cores  from  previous  use  of  the  TMAs,  or  incomplete  data  for 
some  cases,  the  number  (n)  of  tumors  analyzed  for  some  of  the 
markers  was  less  than  450.  Immunohistochemistry  scores 
(IHC  scores)  derive  from  assessment  of  both  staining  intensity 
(scale:  0,  1,  2,  3)  and  percentage  of  positive  cancer  cells  (0%- 
100%).  These  2  scores  were  multiplied  to  generate  an  IHC  or 
H-score  with  a  range  of  0  to  300  as  previously  described  (27). 
Only  nuclear  staining  was  scored.  Briefly,  serial  sections 
(5  pm)  of  the  TMAs  were  stained  with  anti-CARMf  antibodies, 
using  an  automated  tissue  immunostainer  (Discovery  Staining 


Module;  Ventana  Medical  Systems).  The  dilution  of  the  pri¬ 
mary  CARM1  antibody  (Ab)  applied  initially  to  the  slides  was 
1:150,  incubated  for  1  hour  at  42°C  on  the  Ventana  Discovery 
Staining  module,  using  CC1  buffer  and  the  Mild  and  Standard 
protocols  for  Antigen  Retrieval  (AR).  The  VIEW  DAB  kit  and 
reagents  from  Ventana  Medical  Systems  were  used.  Slides 
were  viewed  and  scored  using  standard  light  microscopy. 

Additional  assays  and  methods 

Cell  growth  assays,  quantitative  real-time  PCR  (qRT-PCR), 
cell-cycle  profiling,  CARM1  Ab  characterization,  mouse  xeno¬ 
graft  experiments,  and  microarray  gene  expression  analyses 
are  described  in  Supplementary  Methods  and  Materials.  The 
microarray  GEO  numbers  are  GSE26454  and  GSE26259. 

Statistical  analysis 

One-way  ANOVA  and  ANOVA  single-factor  analysis  was 
applied  and  a  value  of  P  <  0.05  was  regarded  as  statistically 
significant. 

Results 

Overexpression  of  t  ARM  I  in  MCF7  breast  cancer  cell 
line  inhibits  estrogen-dependent  cell  proliferation  and 
anchorage-independent  growth  and  stimulates  the 
expression  of  CDK  inhibitors  p21c,pl  and  p27ki|>l 

MCF7  breast  cancer  cells  stably  overexpressing  CARM1, 
MCF7-CARM1,  were  generated  (25).  MCF7-CARM1  grows  at  a 
slower  rate  than  parental  cells  carrying  empty  vector  (MCF7- 
vector)  as  measured  by  MTT  assays  (Fig.  1A).  The  altered 
growth  rate  was  not  due  to  the  superphysiologic  amounts  of 
CARM1  in  MCF7-CARM1,  because  compared  with  MCF7- 
vector  controls,  MCF7-CARM1  cells  have  only  a  2-fold  increase 
in  CARM1  expression  (Fig.  IB).  Consistent  with  the  growth 
phenotype,  the  expression  of  the  CDK  inhibitors  p21clpl  and 
p27kipl  was  elevated  in  MCF7-CARM1  treated  with  ^beta- 
estradiol  (E2)  (Fig.  IB).  Also  the  expression  of  p21clpl  and 
P27kipl  was  stimulated  by  E2  in  a  time-dependent  manner  in 
MCF7-CARM1  cells  (Supplementary  Fig.  1)  but  not  in  parental 
MCF7  cells,  suggesting  that  ERa  and  CARM1  are  involved  in 
regulating  their  expression.  The  effect  of  CARMf  on  ancho¬ 
rage-independent  cell  growth  was  determined  using  soft-agar 
assays.  E2  stimulates  colony  formation  of  MCF7-vector  cells; 
in  contrast,  no  colonies  were  formed  in  soft  agar  with  MCF7- 
CARM1  (Fig.  1C).  This  result  suggests  that  overexpressing 
CARM1  in  MCF7  may  inhibit  anchorage-independent  growth. 
In  contrast  to  MCF7,  no  growth  effects  were  detected  by 
overexpressing  or  knocking  down  CARM1  in  MDA-MB-231, 
an  ERa-negative  breast  cancer  cell  line  (Fig.  ID  and  E). 
Consistent  with  it  being  ERa  negative,  the  growth  rate  of 
MDA-MB-231  was  E2  independent  (Fig.  ID).  Similarly,  over¬ 
expressing  CARM1  exhibits  no  growth  effect  on  MDA-MB-468, 
another  ERa-negative  breast  cancer  cell  line  (Supplementary 
Fig.  2),  supporting  the  notion  that  the  growth  inhibitory  effect 
of  CARM1  in  MCF7  is  ERa  dependent.  The  growth  inhibitory 
effect  of  CARM1  was  further  validated  in  another  ERa-positive 
breast  cancer  cell  line  ZR-75  (Supplementary  Fig.  3).  p21clpl 
has  been  reported  to  induce  both  cell  cycle  arrest  and  cell 
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Figure  1.  Overexpression  of 
CARM1  inhibits  growth  and 
colony  formation  of  MCF7  while 
exhibiting  no  effect  on  MDA-MB- 
231.  A,  the  CARM1- 
overexpressing  MCF7  cell  line 
MCF7-CARM1  grew  at  a  slower 
rate  than  the  MCF7-vector  control 
as  measured  by  MTT.  Error  bars, 
SD  (n  =  3).  B,  p21cip1  and  p27kip1 
were  expressed  at  higher  levels  in 
MCF7-CARM1  than  in  MCF7- 
vector  cells  in  the  presence  of  1 0 
nmol/L  E2.  C,  overexpression  of 
CARM1  inhibited  colony  formation 
of  MCF7-CARM1  in  soft  agar  in 
the  presence  of  1 0  nmol/L  E2.  The 
inset  shows  colonies  under  higher 
amplification  in  MCF7-vector 
cells.  D,  neither  overexpressing 
nor  knocking  down  CARM1  has 
growth  effects  on  MDA-MB-231. 
Each  assay  was  conducted  in 
triplicate  and  the  error  bars 
represent  SD.  E,  Western  blotting 
showed  that  CARM1  was 
overexpressed  in  MDA-MB-231  - 
CARM1  and  knocked  down  in 
MDA-MB-231 -shCARMI  cells. 


differentiation  in  various  carcinomas  (28,  29).  The  findings 
that  p21clpl  expression  is  increased  by  E2  in  the  presence  of 
exogenous  CARM1  (Supplementary  Fig.  1)  raises  the  possibi¬ 
lity  that  CARM1  may  inhibit  breast  cancer  growth  by  mod¬ 
ulating  key  ERa  target  genes  involved  in  cell-cycle  control  and 
differentiation. 

CIA K\1 1  decreases  estrogen-dependent  breast  cancer 
cell  growth  and  S-phase  entry 

To  eliminate  the  possibility  that  the  growth  effects  of 
CARM1  in  MCF7-CARM1  cells  could  be  attributed  to  addi¬ 
tional  changes  during  retroviral  integration  events,  we  gen¬ 
erated  2  inducible  MCF7  stable  cell  lines:  one  overexpresses 
CARM1  (MCF7-Tet-on-CARMl)  and  the  other  expresses 
CARM1  shRNA  (MCF7-Tet-on-shCARMl)  under  the  control 
of  a  tetracycline-inducible  promoter.  These  stable  cell  lines 
serve  as  gain-of-function  and  loss-of-function  cell  culture 
models  for  studying  the  effects  of  CARM1  in  estrogen-depen- 
dent  breast  cancer  growth.  Cells  were  preincubated  with 
doxycycline  (Dox)  for  4  days  to  induce  or  knockdown  CARM1 
expression,  followed  by  E2  treatment  for  24  hours.  With  either 
cell  line,  E2  alone  has  no  significant  effect  on  CARM1  expres¬ 
sion  at  both  mRNA  and  protein  levels  (Fig.  2A  and  B).  Dox 
could  increase  CARM1  expression  in  MCF7-Tet-on-CARMl 
cells  by  2-fold  (Fig.  2A)  and  reduce  CARM1  to  greater  than  90% 
in  MCF7-Tet-on-shCARMl  cells  (Fig.  2B).  E2  has  no  additional 
effect  on  CARM1  expression  compared  with  Dox  alone  when 


both  are  present.  The  2  cell  lines  were  employed  to  measure 
cell  growth  by  MTT  assays  under  4  treatment  conditions: 
vehicle,  E2,  Dox,  or  combination  of  Dox  and  E2  for  4  time 
points  (24,  48,  72,  and  96  hours).  As  expected,  E2  treatment 
significantly  increases  MCF7  cell  growth  starting  from  day  2  ( P 
value  <0.001;  Fig.  2C).  Overexpression  of  CARM1  by  Dox 
treatment  alone  decreased  MCF7  cell  growth  (Fig.  2C).  Sta¬ 
tistical  analysis  of  3  independent  experiments  suggested  that 
overexpression  of  CARM1  by  Dox  treatment  significantly 
repressed  E2-induced  cell  growth  in  2  individual  clones,  clone 
7  (Fig.  2C)  and  clone  13  (Supplementary  Fig.  4).  This  is  in 
contrast  to  the  nonstatistically  significant  effect  of  Dox  upon 
E2-induced  cell  growth  in  MCF7-Tet-on-shCARMl  cells  (P  > 
0.05;  Fig.  2D)  and  a  CARM1  stable  knockdown  MCF7  (MCF7- 
shCARMl)  cell  line  expressing  shRNA  targeting  a  different 
sequence  of  human  CARM1  (Supplementary  Fig.  5,  P  =  0.04). 

The  main  proliferative  action  of  E2  in  breast  cancer  is  to 
promote  cell-cycle  progression  during  G ,  / S  transition  (3). 
Since  CARM1  can  induce  expression  of  p21clpl  and  p27Kipl, 
which  are  negative  regulators  of  the  cell  cycle,  and  inhibit  E2- 
dependent  growth,  we  determined  whether  CARM1  would 
interfere  with  E2-induced  cell  cycle  progression.  MCF7-Tet- 
on-CARMl  cells  were  preincubated  with  Dox  for  4  days, 
followed  by  E2  treatment  for  24  hours.  Fluorescence-activated 
cell-sorting  analysis  of  MCF7-tet-on-CARMl,  using  propidium 
iodide  labeling,  showed  that  E2-induced  S-phase  entry  was 
inhibited  by  overexpressing  CARM1  (Fig.  2E).  This  result  was 
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Figure  2.  CARM1  inhibits  E2-dependent  MCF7  cell  growth  and  S-phase  entry.  MCF7-Tet-on-CARM1  (A)  and  MCF7-Tet-on-shCARM1  (B)  tetracycline- 
inducible  cell  lines  were  constructed  as  gain-of-function  and  loss-of-function  cell-based  models  to  modulate  the  endogenous  level  of  CARM1 .  CARM1  was 
detected  by  Western  blotting  and  mRNA  was  determined  by  qRT-PCR.  *,  statistically  significant  P  value.  C,  inhibition  of  E2-dependent  cell  growth  by 
overexpressing  CARM1  (+Dox)  in  MCF7-Tet-on-CARM1 .  Bars,  SD  (n  =  9).  D,  knocking  down  CARM1  in  MCF7-Tet-on-shCARM1  (+Dox)  did  not  affect  the 
basal  and  E2-dependent  growth  of  MCF7  cells  in  MTT  assays.  Bars,  SD  (n  =  9).  E,  E2  increased  the  S-phase  entry,  an  event  inhibited  by  CARM1 
overexpression  in  MCF7-Tet-on-CARM1  cells.  F,  the  effects  of  CARM1  overexpression  or  knockdown  on  E2-dependent  S-phase  entry  (P>  0.05)  measured  by 
BrdU  labeling.  Error  bars,  SD  (n  =  3). 


validated  by  bromodeoxyuridine  (BrdU)  labeling  (Fig.  2F). 
Although  E2  and  E2  +  Dox  both  increased  S-phase  entry  as 
compared  with  that  of  the  vehicle  (P  value  <0.001  and  0.0015, 
respectively),  results  from  3  independent  experiments  showed 
that  the  percentage  of  S-phase  entry  induced  by  Dox  +  E2  was 
significantly  decreased  compared  with  E2  treatment  alone 
(P  =  0.0013),  indicating  that  overexpression  of  CARM1 
decreased  E2  induction  of  S-phase  entry.  In  contrast,  in 
MCF7-Tet-on-shCARMl,  Dox  +  E2  treatment  displayed  no 
difference  in  S-phase  entry  compared  with  E2  alone  and  both 
treatment  groups  induced  S-phase  entry  compared  with  the 
vehicle  treatment  (P  =  0.0014).  In  either  MCF7-Tet-on-CARMl 
or  MCF7-Tet-on-shCARMl  cells,  Dox  alone  had  no  significant 
effect  on  S-phase  entry  (Fig.  2F).  These  data  suggest  that 
overexpression  of  CABM1  can  inhibit  E2-stimulated  cell 
growth  through  modulating  cell  cycle,  while  loss  of  CARM1 
could  not  further  accelerate  E2-stimulated  growth  within 
4  days  of  treatment. 


Changes  of  cell  morphology  and  differentiation  marker 
expression  by  increasing  CAKM1  level  in  MCF7  cells 

In  addition  to  the  growth  inhibitory  effects  of  CARM1,  we 
noticed  that  MCF7  cells  stably  overexpressing  CARM1  dis¬ 
played  a  distinct  cell  morphology  from  that  of  MCF7-vector 
cells  (Fig.  3A)  and  exhibited  increased  cell  adhesion 
(requires  longer  trypsin  treatment  time).  Next  we  investi¬ 
gated  desmoplakin  1  (DSP1)  expression,  a  known  differen¬ 
tiation  marker  of  epithelial  cells  that  plays  an  essential  role 
in  maintaining  cell  adhesion  and  differentiation  (30,  31),  and 
a  CARM1  target  gene  identified  in  this  study.  Three  inde¬ 
pendent  experiments  showed  that  E2  significantly 
decreased  DSP1  mRNA,  which  was  reversed  by  overexpres¬ 
sing  CARM1  in  MCF7-Tet-on-CARMl  cells  (Fig.  3B).  In 
addition,  induction  of  2  additional  differentiation  markers, 
GATA-3  and  E-cadherin,  by  overexpressing  CARM1  was 
observed  in  MCF7-Tet-on-CARMl  (Fig.  3C)  by  Western 
blotting.  These  data  suggested  that  growth  inhibitory 
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Figure  3.  Overexpression  of 
CARM1  induces  morphologic 
changes  in  MCF7  characteristic  of 
a  differentiated  phenotype.  A, 
MCF7-CARM1  and  MCF7-vector 
cells  exhibit  different  morphology. 
B,  DSP1  mRNA  level  was 
significantly  repressed  by  E2 
treatment  and  overexpression  of 
CARM1  could  reverse  E2- 
mediated  DSP1  repression.  Error 
bars,  SD  from  triplicate 
experiments.  C,  Western  blotting 
showed  that  overexpression  of 
CARM1  restored  protein  levels  of 
E2-repressed  E-cadherin  and 
GATA-3  at  12  hours  after  E2 
treatment. 


function  of  CARM1  may  be  accompanied  by  the  induction  of 
cell  differentiation. 

CARM1  levels  in  MCF7  cells  modulate  the  ERa  gene 
signature 

Since  CARM1  inhibits  E2-dependent  growth  of  MCF7  cells 
and  induces  a  morphology  change,  we  determined  the  global 
effect  of  CARM1  on  E2-dependent  ERa  gene  signature  by 
microarray  analyses  of  CARM1  gain-of-function  and  loss-of- 
function  cell  lines  treated  with  vehicle  or  E2.  MCF7-inducible 
cells  were  treated  under  4  conditions:  DMSO,  Dox,  E2,  and  E2 
plus  Dox.  The  gene  signature  as  calculated  by  fold  change  was 
normalized  to  vehicle  control  (Fig.  4A).  Microarray  analysis  of 
MCF7-Tet-on-CARMl  reveals  that  E2  upregulated  expression 
of  313  genes  and  downregulated  157  genes  ( P  <  0.05,  fold 
change  >1.6).  Overexpression  of  CARM1  drastically  altered 
E2-regulated  gene  signatures.  Approximately  16%  of  E2- 
induced  genes  including  cell-cycle  regulators  (e.g.,  c-Myc; 
Fig.  4B)  were  inhibited.  The  most  profound  effect  of  CARM1 
overexpression  on  E2-dependent  signature  was  to  relieve  the 
repression  of  approximately  56%  of  E2-repressed  genes 
(Fig.  4C;  P  <  0.05,  fold  change  <0.6  compared  with  vehicle). 
To  our  knowledge,  CARM1  is  the  only  coactivator  which 
affects  global  expression  of  E2-repressed  genes.  Interestingly, 
gene  ontology  (GO)  of  the  affected  genes  suggested  that  most 
E2-repressed,  CARMl-activated  genes  are  involved  in  cell 
differentiation  and  development  (Fig.  4C).  The  ability  of 
CARM1  both  to  inhibit  E2-dependent  growth  and  S-phase 
entry  and  to  modulate  E2-dependent  genes  involved  in  cell- 


cycle  progression,  cell  differentiation,  and  development  sup¬ 
ports  a  role  of  CARM1  in  modulating  the  programming  of  E2- 
dependent  cellular  processes  (i.e„  regulating  the  balance 
between  cell  differentiation  and  proliferation). 

Since  CARM1  has  putative  effects  on  ERa-dependent  pro¬ 
liferation  and  differentiation,  we  applied  qRT-PCR  to  validate 
the  effect  of  CARM1  overexpression  on  6  differentially 
expressed  genes  identified  by  microarray.  p21clpl  and  p27ldpl 
are  known  to  inhibit  breast  cancer  growth  (32).  Cyclin  G2  is  an 
ERa  target  gene  and  a  negative  regulator  of  cell  cycle  (33). 
Among  genes  involved  in  cell  differentiation,  GATA-3  is  an  ERa 
target  gene  and  prodifferentiation  marker  of  breast  cancer  (8, 
34).  MAZ  is  a  transcriptional  factor  (35),  and  KRTAP10.12  is  a 
potential  prodifferentiation  marker.  As  shown  in  Figure  4D,  E2 
alone  significantly  decreased  cyclin  G2  and  KRTAP10.12  mRNA 
but  not  p21cip,  p27kipl,  MAZ,  and  GATA-3  mRNA  after  4-hour 
treatment.  However,  overexpression  of  CARM1  relieved  E2 
repression  of  cyclin  G2  and  KRTAP10.12  and  significantly 
induced  p21clpl,  p27kipl,  MAZ,  and  GATA-3  regardless  of  E2 
(Fig.  4D)  at  mRNA  level.  Consistently,  the  protein  levels  of 
GATA-3,  E-cadherin  (Fig.  3C),  and  P21cipl  and  p27kipl  (Fig.  IB) 
were  also  increased  by  CARM1  overexpression  and  E2  treat¬ 
ment.  These  results  validate  our  microarray  data  and  reinforce 
the  hypothesis  that  CARM1  may  antagonize  the  proliferative 
action  of  estrogen  in  breast  cancer  cells  by  activating  multiple 
cell-cycle-negative  regulators  and  prodifferentiation  genes.  It 
is  worth  noting  that  p21cipl  induction  requires  both  CARM1 
overexpression  and  E2  treatment.  In  contrast,  overexpressing 
CARM1  alone  is  sufficient  to  induce  genes  such  as  p27kipl. 
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Figure  4.  Overexpression  of  CARM1  modulates  E2-dependent  gene  signature.  A,  the  ratios  of  normalized  intensities  for  Dox,  E2,  or  Dox  +  E2-treated  samples 
versus  that  of  samples  treated  with  control  vehicle  (Dox  vs.  DMSO,  E2  vs.  DMSO,  and  DOX  +  E2  vs.  DMSO)  were  used  to  show  the  activation  or  repression. 
Heat  map  of  gene  expression  calculated  as  fold  changes  compared  with  vehicle  indicated  that  CARM1  -induced  genes  (+Dox)  are  largely  nonoverlapping  with 
E2-activated  genes.  Among  CARM1  -repressed  genes,  many  are  activated  by  E2  (see  blowup  of  the  heat  map),  indicating  that  overexpressing  CARM1  can 
inhibit  some  E2-activated  genes.  B,  pie  graph  shows  that  among  all  E2-upregulated  genes,  16%  of  them  are  downregulated  by  CARM1  overexpression. 
The  bottom  of  chart  shows,  among  all  CARM1  -downregulated,  E2-upregulated  genes,  the  percentage  of  genes  in  each  molecular  function  category.  GO  of 
the  affected  genes  was  determined  by  GOStat  tool  online  (http://gostat.wehi. edu.au/cgi-bin/goStat. pi), "goa_human"  database,  where  %  represent  the 
percentage  of  the  affected  genes  that  belong  to  each  represented  category.  C,  pie  graph  shows  that  among  all  E2-downregulated  genes,  56%  of  them  are 
upregulated  by  CARM1  overexpression.  The  bottom  chart  shows,  among  all  CARMI-upregulated,  E2-downregulated  genes,  the  percentage  of  genes  in 
each  molecular  function  category.  D,  Q-RT  PCR  analyses  of  p21clp1,  p27klp1,  cyclin  G2,  MAZ,  KRTAP10.12,  and  GATA-3  expression  in  MCF7-Tet-on-CARM1 . 
Error  bars,  SD  from  3  independent  experiments. 


suggesting  that  CARM1  may  regulate  some  genes  in  hormone- 
deprived  conditions. 

The  global  effects  of  CARM1  on  ER  target  genes  were  next 
examined  in  the  loss-of-function  model,  MCF7-Tet-on- 
shCARMl  under  aforementioned  conditions.  The  heat 
map  of  the  fold-change  gene  signature  relative  to  vehicle 
control  is  shown  in  Figure  5A  and  additional  description  can 
be  found  in  Supplementary  Methods.  Using  Agilent  array 
platform  herein,  CARM1  shRNA  expression  upregulated  62 
genes  and  downregulated  2,122  genes  ( P  <  0.05,  fold  change 
>1.6  and  <0.6  compared  with  vehicle,  respectively).  E2 
treatment  upregulated  780  genes  and  downregulated  5,099 
genes  (P  <  0.05,  fold  change  >1.6).  Interestingly,  the  genes 
affected  by  loss-of-CARMl  largely  overlapped  with  those 
affected  by  E2  in  wild-type  cells.  Further  microarray  analysis 


showed  that  65%  of  genes  activated  by  knocking  down 
CARM1  are  also  activated  by  E2  (Fig.  5B)  and  75%  of  genes 
repressed  by  CARM1  knockdown  are  also  repressed  by  E2 
(Fig.  5C).  GO  of  genes  affected  by  Dox  and  E2  treatment  also 
overlap  (Fig.  5B  and  C).  Among  these  genes,  a  majority  are 
involved  in  metabolism,  development,  protein  binding,  and 
gene  expression  (Fig.  5B  and  C).  These  data  further  support 
the  notions  that  CARM1  is  a  global  regulator  of  E2-respon- 
sive  genes  in  breast  cancer  cells  and  profoundly  affects 
estrogen-mediated  processes.  We  also  validated  the  effect 
of  loss  of  CARM1  on  p21cipl,  p27kipl,  cyclin  G2,  MAZ,  GATA- 
3,  and  KRTAP10.12  mRNA  expression.  Loss  of  CARM1  sig¬ 
nificantly  repressed  p21clpl,  p27kipl,  cyclin  G2,  MAZ,  GATA-3, 
KRTAP10.12,  and  DSP1  at  mRNA  levels,  similar  to  the  effect 
of  E2  in  MCF7-Tet-on-shCARMl  (Fig.  5D  and  Supplementary 
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Figure  5.  Knocking  down  CARM1  shared  common  gene  signatures  with  that  of  E2  treatment  in  MCF7-Tet-on-shCARM1 .  A,  heat  map  showed  that  CARM1 
shRNA  (+Dox)-activated  and  -repressed  genes  were  largely  overlapping  with  those  regulated  by  E2.  A  blowup  of  the  heat  map  illustrated  that  many  Dox- 
activated  genes  (knocking  down  CARM1)  were  also  activated  by  E2.  B,  pie  graph  shows  that  approximately  65%  of  CARM1  shRNA-activated  genes  were 
also  activated  by  E2;  approximately  35%  of  knocking  down  CARM1 -activated  genes  was  E2-nonresponsive  (labeled  as  others).  Gene  ontology  illustrated 
the  top  affected  gene  categories  by  expressing  CARM1  shRNA  or  E2  treatment.  C,  pie  graph  shows  that  among  all  downregulated  genes  by  CARM1 
knockdown,  75%  of  them  are  E2-repressed  genes  and  25%  of  them  are  E2-nonresponsive  (labeled  as  others).  The  bottom  chart  shows,  among  all 
downregulated  genes  by  E2  treatment  or  CARM1  knockdown,  the  percentage  of  genes  in  each  molecular  function  category.  D,  qRT-PCR  analyses  of  ER 
target  genes  in  MCF7-Tet-on-shCARM1 .  Error  bars,  SD. 


Fig.  6A).  In  agreement  with  the  mRNA  results,  cyclin  G2, 
GATA-3,  and  E-cadherin  (Supplementary  Fig.  6B)  were 
decreased  at  protein  levels  with  the  loss  of  CARM1.  Since 
both  cyclin  G2  (33)  and  GATA-3  (36)  are  ERa  target  genes, 
CARM1  may  antagonize  E2  action  via  ERa  during  repro¬ 
gramming  of  ERa-dependent  differentiation  and  prolifera¬ 
tion  processes.  Fold  changes  of  key  cell-cycle  regulators  and 
genes  involved  in  cell  differentiation  in  MCF7-Tet-on- 
shCARMl  are  listed  in  Supplementary  Table  SI.  Overall, 
our  data  suggest  that  loss  of  CARM1  induces  gene  signatures 
resembling  those  affected  by  E2  and  CARM1  is  a  regulator  of 
E2-dependent,  key  cell  cycle  progression,  and  differentiation 
genes.  Collectively,  the  microarray  analyses  using  CARM1 
gain-of-function  and  loss-of-function  cell  models  reveal  that 
CARM1  is  a  unique  ER  coactivator  that  profoundly  affects 
the  balance  of  genes  involved  in  cellular  differentiation  and 
proliferation  (i.e.,  inhibit  cell  growth  and  promote  cell 
differentiation). 


Knocking  down  of  CARM1  increased  E2-dependent 
tumor  growth  in  an  MCF7  xenograft  mouse  model 

To  examine  the  effects  of  CARM1  in  vivo,  we  transplanted 
MCF7-Tet-on-CARMl  shRNA  cells  in  nude  mice.  The  design  of 
the  xenograft  experiment  is  shown  in  Figure  6A,  representing 
one  of  triplicate  experiments.  We  first  validated  that  the 
growth  of  xenografted  tumors  was  E2-dependent  because 
no  growth  or  only  tiny  tumors  developed  in  the  negative 
control  group  not  receiving  estrogen.  Tumors  collected  from 
mice  engrafted  with  MCF7-Tet-on-shCARMl  cells  and  receiv¬ 
ing  Dox  showed  a  reduction  of  CARM1  expression  at  the 
mRNA  and  protein  levels  (Fig.  6B).  Knocking  down  CARM1 
increased  the  size  of  E2-induced  tumors  (Fig.  6C)  and  was 
associated  with  a  modest  increase  in  BrdU  labeling.  The 
differential  rate  of  BrdU  labeling  for  xenografted  tumors 
was  further  increased  in  mice  receiving  a  higher  dose  of  E2 
and  that  was  associated  with  higher  mitotic  index  (Fig.  6D  and 
E).  All  the  data  suggest  that  knocking  down  CARM1  enhances 
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Figure  6.  Knocking  down  CARM1  increased  E2-dependent  tumor  growth  in  xenografted  mice.  A,  schematic  design  of  the  xenograft  study  using  MCF7-Tet-on- 
shCARMI  cell  line.  B,  CARM1  mRNA  and  protein  levels  were  decreased  in  tumors  from  Dox  recipient  mice.  CARM1  mRNA  in  tumors  was  quantified 
using  qRT-PCR  (vehicle:  n  =  7;  Dox:  n  —  10);  protein  was  visualized  by  IHC.  C,  increased  tumor  volume  in  CARM1  knocked  down  mice.  A  representative 
experiment  showed  a  higher  tumor  volume  with  knocking  down  CARM1  (n  =  10)  than  E2-alone-induced  tumors  (n  =  8).  The  red  arrows  point  to  tumors 
in  2  representative  mice.  D,  BrdU  and  hematoxylin  and  eosin  staining  of  representative  tumor  samples  from  vehicle  and  Dox-treated  mice.  E,  mitotic  index 
of  representative  tumor  samples  from  vehicle  and  Dox-treated  mice  implanted  with  a  high-dose  E2  pellets.  F,  correlation  of  p21olp1  and  E-cadherin  protein 
level  by  IHC. 


E2-dependent  proliferation  of  breast  cancer  cells  in  vivo.  Since 
CARM1  inhibits  E2-dependent  growth  by  modulating  negative 
cell-cycle  regulators  p21clpl,  p27kipl,  and  cyclin  G2  and  pro¬ 
differentiation  genes,  we  examined  the  relationship  between 
p21clpl  and  E-cadherin,  a  differentiation  marker,  in  E2- 
induced  xenografted  tumors.  A  direct  correlation  was 
observed  between  p21clp  and  E-cadherin  expression  in 
tumors  derived  from  xenografts  (Fig.  6F),  suggesting  inhibi¬ 
tion  of  cell  growth  and  induction  of  differentiation  are  coher¬ 
ent  processes  in  ERa-positive  tumors. 

(  ARM  I  expression  in  human  breast  tumor  biopsy 
samples 

Our  rabbit  polyclonal  CARM1  Ab  was  determined  to  be 
specific  because  it  detects  both  nuclear  and  cytoplasmic 


CARM1  in  normal  breast  tissues  and  breast  tumors  while 
exhibiting  no  activity  toward  mouse  embryonic  fibroblasts 
derived  from  CARM1  knockout  mice  (MF.F  ;  Supplemen¬ 
tary  Fig.  7).  CARM1  expression  was  determined  by  IHC  in  ER+ 
breast  tumor  TMAs  available  in  the  MBTB  (26, 37).  Statistically 
significant  correlations  between  ERa  expression  as  deter¬ 
mined  by  IHC  ( n  =  310,  Spearman  r  =  0.324,  P  <  0.0001) 
and  tumor  grade  ( n  =  328,  Spearman  r  =  —0.159,  P  =  0.004) 
were  found.  Significantly  higher  CARM1  expression  as  deter¬ 
mined  by  IHC  score  was  found  in  tumors  with  higher  ERa 
expression  than  in  those  with  lower  ERa  expression  (Fig.  7A). 
Significantly  higher  CARM1  expression  was  found  in  lower- 
grade  (3, 4)  tumors  as  well  (Fig.  7B  and  C).  In  addition,  CARM1 
expression  was  positively  correlated  with  ERa  levels  in  ER- 
positive,  node-negative  human  breast  tumors  (P  <  0.0001; 
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Figure  7.  The  expression  level  of 
CARM1  positively  correlates  with 
ERoc  level  and  inversely  correlates 
with  tumor  grade  in  ER-positive 
breast  tumors.  A,  CARM1 
expression  level  is  directly 
correlated  with  ERoc  expression 
level  in  more  than  300  human 
breast  tumor  samples,  using 
TMAs.  B  and  C,  CARM1 
expression  level  is  higher  in  ER- 
positive,  low-grade  tumors  than  in 
high-grade  tumors,  supporting  a 
potential  link  between  CARM1 
and  the  differentiation  status  of 
ERoc-dependent  breast  tumors. 
ERa+ve,  grade  3,  4  (low),  IHC 
score  270,  magnification  x  500; 
ERa+ve,  grade  8,  9  (HIGH),  IHC 
score  0,  magnification  x  500. 


Supplementary  Fig.  8A).  We  also  found  an  inverse  correlation 
between  CARM1  expression  and  tumor  grade  in  ER-positive, 
node-negative  human  breast  tumors  ( P  <  0.0398;  Supplemen¬ 
tary  Fig.  8B).  Collectively,  the  findings  from  clinical  samples 
support  a  role  of  CARM1  in  regulating  ERoc-dependent  differ¬ 
entiation  in  ERa-positive  tumors. 

Discussion 

In  most  cases,  proliferation  and  differentiation  are  inversely 
coupled:  repression  of  proliferation  is  a  prerequisite  for  initia¬ 
tion  of  differentiation  (11).  In  many  cell  types,  however,  cell- 
cycle  arrest  is  necessary  but  not  sufficient  for  differentiation. 
CARM1  seems  to  be  a  unique  ER  coactivator  regulating  both 
processes.  Overexpression  of  CARM1  in  MCF7  cells  results  in 
inhibition  of  E2-dependent  growth  through  inhibition  of  the 
G0/Gx  transition  to  S  phase.  This  is  in  part  due  to  upregulation 
of  key  negative  cell-cycle  regulators  such  as  p21clpl,  p27kipl, 
and  cyclin  G2.  Inhibition  of  E2-dependent  cell  growth  by 
CARM1  is  accompanied  by  morphologic  changes  character¬ 
istic  of  a  more  differentiated  phenotype  and  induction  of 
multiple  differentiation  markers  such  as  GATA-3  and  MAZ. 
This  finding  is  supported  by  previous  reports  that  CARM1  can 
promote  cell  differentiation  in  other  systems  (19-21).  None¬ 
theless,  regulation  of  cell  differentiation  by  CARM1  seems  to 
be  cell-type  and  context  dependent.  In  mouse  embryo  and 
embryonic  stem  cells,  CARM1  was  shown  to  elevate  expres¬ 
sion  of  key  pluripotency  genes  and  delay  their  response  to 
differentiation  signals  (38). 

In  contrast  to  growth  inhibition  by  CARM1  overexpression, 
knocking  down  CARM1  in  MCF7  did  not  alter  E2-dependent 
cell  growth  in  cell  culture,  nor  did  it  affect  E2-induced  S-phase 
entry.  This  observation  contradicts  the  conclusion  by  Frietze 
and  colleagues  that  CARM1  increases  growth  of  MCF7  cells. 
The  discrepancies  may  be  due  to  the  transient  transfection  of 
CARM1  siRNA  throughout  the  cell-cycle  study  by  Frietze  and 
colleagues  (16).  Moreover,  the  authors  measured  the  percen¬ 


tage  of  cells  in  S  4-  G2  +  M  phase  without  distinguishing  the 
percentage  of  cells  in  S  phase.  Also,  in  consistent  with  the 
observation  of  O’Brien  and  colleagues  (21),  we  did  not  observe 
change  of  E2F1  with  CARM1  knockdown  in  contrast  to  Frietze 
and  colleagues  (16).  Interestingly,  and  in  contrast  to  cells 
grown  in  culture,  knocking  down  CARM1  enhanced  E2- 
induced  xenograft  tumors.  This  may  be  due  to  increased 
breast  cancer  cell  interaction  with  the  microenvironment 
which  plays  essential  roles  in  promoting  tumor  growth  in 
animals. 

The  growth  inhibitory  effect  of  CARM1  is  unique  from  that 
of  SRCs.  Knocking  down  SRC2  and  SRC3  but  not  SRC1  inhibits 
growth  of  MCF7  cells  and  decreases  cyclin  D1  expression  (39). 
Overexpression  of  SRC3  also  increases  breast  cancer  cell 
proliferation  and  invasiveness.  Similarly,  SRC-1  promotes 
breast  tumor  metastasis  and  inhibits  tumor  cell  differentia¬ 
tion  (40).  Thus,  the  ERa-dependent,  growth  inhibitory  effect  of 
CARM1  is  unlikely  to  be  mediated  through  SRC1,  2,  and  3. 

Cell-cycle  genes  that  are  regulated  by  E2  or  loss  of  CARM1 
include  cyclin  Dl,  c-Myc,  cyclin  G2,  cyclin  LI,  cyclin  T2, 
p21clpl,  p27ldpl,  pl30,  and  Rb  (Supplementary  Table  SI).  E2 
treatment  alone  significantly  represses  cyclin  G2  (33),  which  is 
reversed  by  overexpressing  CARM1.  Cyclin  Dl  is  a  well-known, 
E2-induced  ERoc  target  gene;  however,  its  expression  is  not 
affected  by  overexpression  of  CARM1  in  the  presence  of  E2,  yet 
knocking  down  CARM1  upregulates  cyclin  Dl  in  MCF7  cells 
(Supplementary  Table  SI).  c-Myc  is  upregulated  by  E2  alone  or 
loss  of  CARM1  (Supplementary  Table  SI)  but  is  not  affected  by 
depletion  of  any  of  the  pl60  coactivators  in  MCF7  cells  (39). 
Thus,  the  mechanism  of  CARM1  regulation  of  cell-cycle 
regulators  is  complex  and  only  partially  depends  on  the 
pl60  coactivators. 

Microarray  gene  expression  analyses  reveal  that  approxi¬ 
mately  16%  of  E2-activated  genes  were  repressed  by  CARM1, 
consistent  with  the  repressive  effects  of  CARM1  on  some  ER 
target  genes  (41).  The  mechanism  of  CARM1 -mediated  repres¬ 
sion  is  unclear.  The  major  effect  of  CARM1  overexpression  was 
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to  relieve  E2-repressed  genes.  CARM1  methyltransferase  activ¬ 
ity  may  be  responsible  for  the  activation  because  we  observed 
increased  H3R17Me2  mark  on  p21clpl  promoter  upon  CARM1 
induction  in  MCF7-Tet-on-CARMl  (data  not  shown),  consis¬ 
tent  with  a  recent  publication  that  CARM1  is  recruited  to 
p21clpl  promoter  (42).  Whether  CARM1  regulates  ERa  target 
genes  via  an  epigenetic  mechanism  remains  to  be  determined. 
Nonetheless,  global  ERa  transcriptional  regulation  by  CARM1 
leads  to  the  induction  of  many  E2-repressed  genes  associated 
with  differentiation. 

Consistent  with  this  finding,  knocking  down  CARM1  shares 
more  than  65%  of  the  E2  gene  signature.  Most  CARM1-  and  E2- 
regulated  genes  are  involved  in  gene  expression,  metabolism, 
cell  cycle,  and  differentiation.  Knocking  down  CARM1  leads  to 
upregulation  of  positive  cell-cycle  regulators  (e.g.,  c-Myc)  and 
downregulation  of  negative  cell-cycle  regulators  (e.g.,  cyclin 
G2).  This  result  suggests  that  loss  of  CARM1  function  may  lead 
to  the  acquisition  of  a  proliferative  phenotype  resembling 
estrogen  stimulation  of  breast  cancer.  Furthermore,  knocking 
down  CARM 1  also  modulates  genes  involved  in  cell  differ¬ 
entiation.  For  example,  combination  of  CARM1  shRNA  and  E2 
treatment  significantly  reduced  the  level  of  PPARy,  which 
induces  terminal  differentiation  of  breast  cancer  (43).  Loss 
of  CARM1  also  significantly  decreases  KRTAP10.12,  an  E2- 
repressed  gene  involved  in  keratin  filament  formation  and 
potentially  in  cell  differentiation  processes  (44,  45).  Collec¬ 
tively,  either  loss  of  CARM1  or  E2  treatment  significantly 
inhibits  expression  of  various  differentiation  markers  (Sup¬ 
plementary  Table  SI).  Overall,  the  gene  expression  data  from 
CARM1  gain-of-function  and  loss-of-function  models  suggest 
that  CARM1  plays  an  important  role  in  regulating  ERa  target 
genes  in  differentiation  and  proliferation. 

Evidence  for  a  functional  interplay  of  ERa  and  CARM  I  was 
explored  in  human  breast  cancer  specimens.  A  direct  correla¬ 
tion  was  observed  between  CARM1  and  ERa  in  ER-positive 
tumors.  Higher  ERa  expression  is  associated  with  less  aggres¬ 
sive  and  more  differentiated  tumors,  and  ER  status  is  known 
to  inversely  correlate  with  histologic  grade  (46).  Our  observa¬ 
tion  contradicts  an  earlier  report  that  CARM1  is  overex¬ 
pressed  in  grade  III  breast  tumors  (23).  The  difference 
could  result  from  analysis  of  RNA  versus  protein  and  the 
sample  size.  In  the  study  by  El  Messaoudi  and  colleagues, 
CARM1  was  analyzed  only  at  the  RNA  level  in  81  human  breast 
tumors  whereas  we  analyzed  CARM1  protein  level  in  more 
than  300  human  breast  tumors. 
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Estrogen  signaling  is  mediated  by  two  estrogen  receptors  (ERs),  ERa  and  ERP,  which  have  unique  roles  in 
the  regulation  of  breast  cancer  cell  proliferation.  ERa  induces  proliferation  in  response  to  estrogen  and 
ER3  inhibits  proliferation  in  breast  cancer  cells,  suggesting  that  ERp  selective  ligands  may  be  beneficial 
for  promoting  the  anti-proliferative  action  of  ERp.  Subtype  selective  ligands  can  be  identified  using 
transcriptional  assays,  but  cell  lines  in  which  ERa  or  ERP  are  independently  expressed  are  required.  Of 
the  available  reporter  cell  lines,  none  have  been  generated  in  breast  cancer  cells  to  identify  subtype 
selective  ligands.  Here  we  describe  the  generation  of  two  isogenic  breast  cancer  cell  lines,  Hs578T- 
ERaLuc  and  Hs578T-ER3Luc,  with  stable  integration  of  an  estrogen  responsive  luciferase  reporter  gene. 
Hs578T-ERaLuc  and  Hs578T-ER3Luc  cell  lines  are  highly  sensitive  to  estrogenic  chemicals  and  ER 
subtype  selective  ligands,  providing  a  tool  to  characterize  the  transcriptional  potency  and  subtype 
selectivity  of  estrogenic  ligands  in  the  context  of  breast  cancer  cells.  In  addition  to  measuring  reporter 
activity,  ERP  target  gene  expression  and  growth  inhibitory  effects  of  ERp  selective  ligands  can  be 
determined  as  biological  endpoints.  The  finding  that  activation  of  ERP  by  estrogen  or  ERP  selective 
natural  phytoestrogens  inhibits  the  growth  of  Hs578T-ERP  cells  implies  therapeutic  potential  for  ERP 
selective  ligands  in  breast  cancer  cells  that  express  ERp. 

©  201 1  Elsevier  Inc.  All  rights  reserved. 


1.  Introduction 

Estrogens  regulate  mammary  gland  growth  and  differentiation, 
ovary  and  uterus  maturation,  and  bone  homeostasis  [1].  The 
physiological  effects  of  estrogens  are  primarily  mediated  by  two 
estrogen  receptors  (ERs),  ERa  and  ERp.  Because  of  the  broad  range 
of  ER  target  tissues  and  the  ligand  dependent  activity  of  the 
receptors,  synthetic  and  natural  estrogens  hold  therapeutic 
promise  in  selectively  targeting  ERs.  Therapies  aimed  at  preventing 
ERa  transcriptional  activation  are  currently  used  for  breast  cancer 
treatment  and  osteoporosis  prevention  [2].  Though  ER(3  is  not 
currently  a  therapeutic  target,  accumulating  evidence  suggests  an 
anti-proliferative  role  for  ERP  in  breast  cancer  [3],  In  the  mammary 
gland,  ERa  and  ERp  play  opposing  roles  in  regulating  growth  and 


Abbreviations:  BRET,  bioluminescence  resonance  energy  transfer;  Cos,  cosmosiin; 
Dox,  doxycycline;  DPN,  diarylpropionitrile;  E2,  1 7fl-esLiadiol;  ER,  estrogen 
receptor;  PPT,  propyl  pyrazole  triol;  ERE,  estrogen  response  element;  ICI,  IC1 
182,780;  Liq,  liquiritigenin. 

*  Corresponding  author  at:  McArdle  Laboratory  for  Cancer  Research,  1400 
University  Ave,  Madison,  Wl  53706,  USA.  Tel.:  +1  608  265  5540; 
fax:  +1  608  262  2824. 

E-mail  address:  wxu@oncology.wisc.edu  (W.  Xu). 

0006-2952/$  -  see  front  matter  ©  201 1  Elsevier  Inc.  All  rights  reserved, 
doi:  1 0. 1 01 6/j.bcp.201 1 .08.026 


differentiation  in  response  to  estrogens;  ERa  promotes  prolifera¬ 
tion  while  ERP  inhibits  ERa-mediated  proliferation  [4-6].  Because 
the  anti-proliferative  action  of  ERp  may  be  enhanced  by  ligand- 
dependent  activation,  the  paradigm  of  ER  targeted  therapies  is 
expanding  towards  the  development  of  ER  subtype  selective 
ligands  [7], 

Though  ERa  and  ERp  share  many  structural  and  transcriptional 
features,  ligands  can  display  subtype  selectivity.  In  classical  ligand 
dependent  transcriptional  activation,  the  receptors  dimerize  upon 
ligand  binding  and  undergo  conformational  changes  to  allow 
cofactor  recruitment.  The  receptors  directly  bind  DNA  most  often 
at  estrogen  response  elements  (EREs),  consisting  of  a  consensus 
GGTCAnnnTGACC  sequence.  ERa  and  ERp  have  97%  identity 
within  the  DNA  binding  domains,  and  the  receptors  bind  similar 
DNA  sequences  with  high  affinity.  Genome  wide  binding  studies  in 
MCF7  breast  cancer  cells  expressing  ERa  or  ERp  independently 
have  shown  that  ERa  and  ERP  bind  similar  sites  in  response  to 
17P-estradiol  (E2);  ~60%  of  ER  binding  sites  contain  full  EREs  and 
~25%  contain  half  EREs  [8], 

The  ligand  binding  pockets  of  ERa  and  ERP  are  relatively 
large,  and  the  receptors  bind  a  wide  array  of  chemicals.  The 
ligand  binding  domains  of  ERa  and  ERP  have  59%  identity,  and 
the  receptors  bind  E2  with  similar  affinities.  Despite  similarities 
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in  their  ligand  binding  domains,  several  ligands  have  modest 
selectivity  for  ERot  or  ER(3  [9],  and  some  synthetic  ligands 
maintain  high  selectivity.  For  example,  propyl  pyrazole  triol 
(PPT)  is  an  ERa  selective  agonist  that  displays  a  400-fold  higher 
binding  affinity  for  ERa  compared  to  ER[3  [10].  Estrogenic 
chemicals  produced  in  plants,  known  as  phytoestrogens,  often 
display  subtype  selectivity  for  ERp.  For  example,  liquiritigenin  is 
a  flavanone  derived  from  Clycyrrhizae  uralensis  that  has  been 
shown  to  have  20-fold  higher  binding  affinity  for  ERP  and  even 
greater  selectivity  in  transcriptional  assays  [11].  Compounds 
such  as  liquiritigenin  often  show  low  binding  affinities  relative  to 
E2,  and  ERP  selective  ligands  with  higher  affinity  and  greater 
selectivity  are  needed  to  fully  elucidate  the  anti-proliferative  role 
of  ERP  in  breast  cancer. 

Mammalian  cell  lines  have  been  developed  to  enable  screening 
for  subtype  selective  ligands.  HeLa  cervical  carcinoma  cells  have 
been  used  to  create  HELN-ERa  and  HELN-ERp,  two  cell  lines  in 
which  ERa  or  ERp,  respectively,  are  constitutively  expressed  with 
stable  integration  of  a  luciferase  reporter  downstream  of  an  ERE 
[12 j.  Human  embryonic  kidney  cells,  HEK293,  have  also  been 
created  using  a  similar  strategy  in  which  ERa  or  ERP  are 
constitutively  expressed  and  human  placental  alkaline  phospha¬ 
tase  downstream  of  the  vitellogenin  ERE  is  stably  integrated  [13]. 
The  only  available  breast  cancer  reporter  cell  line  is  T47D-KBLuc  in 
which  three  tandem  EREs  upstream  of  a  luciferase  reporter  have 
been  stably  integrated  [14].  However,  identification  of  subtype 
selective  ligands  is  prohibited  because  T47D  cells  express  both  ERa 
and  ERp. 

Here,  we  describe  the  generation  of  two  isogenic  reporter 
cell  lines,  Hs578T-ERaLuc  and  Hs578T-ERPLuc,  that  provide  a 
tool  to  characterize  the  transcriptional  potencies  and  subtype 
selectivity  of  estrogenic  compounds  in  the  context  of  breast 
cancer  cells.  These  cell  lines  are  highly  sensitive  to  estrogenic 
ligands  and  can  be  used  to  validate  ER  transcriptional  activation 
by  analysis  of  endpoints  such  as  endogenous  target  gene 
regulation.  Further,  ERp  selective  ligands  induce  ERP-mediated 
reporter  gene  expression,  endogenous  gene  regulation,  and 
growth  inhibition,  suggesting  that  Hs578T-ERPLuc  cells  may  be 
used  to  isolate  ERP  selective  ligands  with  desired  biological 
effects. 

2.  Materials  and  methods 

2.1.  Cell  lines  and  reagents 

Cosmosiin  (apigenin  7-glucoside),  dimethyl  sulfoxide 
(DMSO),  E2,  and  diethylstilbestrol  (DES)  were  obtained  from 
Sigma  (St.  Louis,  MO);  DPN,  PPT,  and  IC1  182,780  were  obtained 
from  Tocris  (Ellinsville,  MO);  liquiritigenin  was  obtained  from 
Chromadex  (Irvine,  CA).  Doxycycline  (Dox)  was  obtained 
from  Clontech.  Hygromycin  B,  blasticidin  S,  zeocin,  NaCl,  sodium 
dodecyl  sulfate  (SDS),  and  dithiothreitol  (DTT)  were  obtained 
from  Research  Products  International  (Mount  Prospect,  IL). 
Triton  X-100  was  obtained  from  Fisher  (Fair  Lawn,  NJ);  protease 
inhibitors  were  obtained  from  Roche  Scientific  (Basel, 
Switzerland);  benzonase  was  obtained  from  Novagen  (San 
Diego,  CA).  All  other  chemicals  were  obtained  from  Sigma  (St. 
Louis,  MO). 

Cell  culture  media  were  obtained  from  Invitrogen  (Carlsbad, 
CA).  MCF7  and  HEK293  cells  were  cultured  in  DMEM  +  10%  fetal 
bovine  serum  (FBS;  Gemini  Bio  Products,  West  Sacramento,  CA) 
at  37  °C  and  5%  C02.  Hs578T-ERa  and  Hs578T-ERP  were 
previously  created  by  Secreto  and  coworkers  [15].  These  cells 
were  cultured  at  37  °C  and  5%  C02  in  DMEM/F12  supplemented 
with  L-glutamine,  10%  Tet-system  approved  FBS  (Clontech 
Mountain  View,  CA),  500  mg/L  zeocin  and  5  mg/L  blasticidin  S. 


2.2.  Generation  of  Hs578T-ERaLuc  and  Hs578T-ERfiLuc  reporter  cell 
lines 

Stable  reporter  cell  lines  were  created  using  a  modified  pGL4.32 
reporter  (Promega,  Madison,  WI)  which  contains  the  luc2P  reporter 
and  hygromycin  resistance.  The  pGL4.32  vector  was  digested  with 
Nhe  I  and  Hind  III  (New  England  Biolabs,  Ipswich,  MA)  and  three 
consensus  EREs  spaced  by  three  nucleotides  were  cloned  upstream 
of  luc2P  using  the  following  oligonucleotides:  5'-CTA  GCG  GTC  ACA 
GTG  ACC  TGC  GAG  GTC  ACA  GTG  ACC  TGC  GAG  GTC  ACA  GTG  ACC 
TGC  GA-3'  and  5'-AGC  TTC  GCA  GGT  CAC  TGT  GAC  CTC  GCA  GGT 
CAC  TGT  GAC  CTC  GCA  GGT  CAC  TGT  GAC  CG-3'.  Successful  cloning 
was  verified  by  complete  sequencing  and  the  vector  was 
designated  pGL4.3xERE.  Estrogen  responsiveness  was  validated 
by  batch  transfecting  HEK293  cells  with  2  ng  of  CMX-ERa  or  CMX- 
ERp,  45  ng  pGL4.3xERE  vector,  and  40  ng  CMX-P-galactosidase  per 
well  of  a  48  well  plate.  Cells  were  incubated  24  hr  to  allow  protein 
expression  before  the  addition  of  the  indicated  ligands.  After  24  hr 
of  ligand  treatment,  cells  were  lysed,  firefly  luciferase  substrate 
(Promega)  was  added,  and  luminescence  was  measured  on  a  Victor 
X5  microplate  reader  (Perkin  Elmer,  Waltham,  MA)  using 
luminescence  detection  and  a  700  nm  filter.  To  normalize  data 
for  transfection  efficiency,  (3-galactosidase  expression  was  ana¬ 
lyzed  using  the  Tropix  [3-galactosidase  detection  kit  (Applied 
Biosystems,  Foster  City,  CA).  Luciferase  counts  were  normalized  to 
(3-gal  counts  in  each  well. 

After  characterizing  the  pGL4.3xERE  stable  reporter  vector, 
Hs578T-ERa  and  Hs578T-ER(3  cells  were  transfected  with  10  pig  of 
the  vector  and  selected  in  125  pig/mL  hygromycin  B  for  4  weeks. 
Individual  colonies  were  selected  using  3  mm  cloning  discs, 
expanded,  and  screened  for  estrogen  induced  luciferase  expres¬ 
sion.  One  clone  from  each  cell  line  was  selected  for  further 
characterization,  referred  to  here  as  Hs578T-ERaLuc  and  Hs578T- 
ERPLuc. 

2.3.  Quantitative  western  blots  and  ligand  binding  assays 

For  quantitative  western  blots,  cells  were  split  in  phenol  red  free 
DMEM/F12  +  5%  SFS  and  treated  with  50  ng/mL  Dox  or  vehicle 
(water)  24  hr  later.  After  48  hr  treatment,  cells  were  collected  by 
trypsinization,  washed  with  Dulbecco’s  phosphate  buffer  saline 
(Invitrogen),  and  lysed  by  suspension  in  lysis  buffer  (50  mM  Tris  pH 
8.0,  400  mM  NaCl,  10%  glycerol,  0.5%  triton  X-100,  protease 
inhibitors,  and  benzonase).  After  centrifugation,  total  protein  was 
quantified  using  BioRad  Protein  Assay  (BioRad),  and  40  pug  of  protein 
was  resolved  using  SDS-PAGE  and  8%  polyacrylamide  gels.  Proteins 
were  transferred  to  a  nitrocellulose  membrane  for  1.5  hr  at  0.35  A. 
Membranes  were  blocked  with  5%  nonfat  milk  and  incubated 
overnight  with  1:1000  anti-FLAG-M2  antibody  (Sigma)  or  1:5000 
anti-(3-Actin  (Sigma)  at  4  °C.  Membranes  were  then  incubated  with 
IRDye  800CW  goat-anti-mouse  IgG  secondary  antibody  (Licor 
Biosciences,  Lincoln,  NE)  for  1  hr  at  room  temperature  and  visualized 
on  a  Licor  Odyssey  near-infrared  gel  reader  (Licor  Biosciences). 

For  ligand  binding  assays,  Hs578T-ERaLuc  and  Hs578T-ER[3Luc 
cells  were  cultured  in  phenol  red  free  DMEM/F12  +  10%  6x 
charcoal  stripped  FBS  (SFS)  for  3  days  prior  to  the  assay  to  remove 
residual  estrogens  from  the  cells.  At  90%  confluence,  cells  were 
collected,  resuspended  in  phenol  red  free  DMEM/F12  +  5%  SFS,  and 
plated  at  a  density  of  105  cells/well  on  a  24  well  plate  in  the 
presence  or  absence  of  50  ng/mL  Dox.  After  24  hr,  cells  were 
labeled  in  triplicate  with  20  nM  [3H]-E2  (89.2  Ci/mmol  specific 
activity,  Perkin  Elmer)  in  the  presence  or  absence  of  450  p,M  DES 
cold  competitor  for  2  hr  at  37  °C  and  5%  C02.  Labeled  cells  were 
washed  3  times  with  cold  PBS  +  0.1%  BSA  and  lysed  with  500  piL 
SDS  lysis  buffer  (0.5%  SDS,  0.05  M  Tris-HCl  pH  8.0,  and  1  mM  DTT). 
Total  cell  lysate  (400  p,L)  was  mixed  with  5  mL  liquid  scintillation 
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cocktail  and  [3H]  bound  radioactivity  was  liquid  scintillation 
counted  for  5  min.  Two  additional  wells  of  each  condition  were 
used  to  count  the  cell  number  and  determine  the  total  protein 
using  RC  DC  protein  assay  (BioRad,  Hercules,  CA). 

2.4.  Luciferase  assays 

Hs578T-ERaLuc  and  Hs578T-ER(3Luc  cells  were  cultured  in 
phenol  red  free  DMEM/F12  +  10%  SFS  for  3  days  prior  to  the  assay 
to  remove  residual  estrogens  from  the  cells.  Cells  were  seeded  in 
triplicate  at  a  density  of  104  cells/well  on  white  96  well  tissue 
culture  plates  (Fisher)  in  phenol  red  free  DMEM/F12  +  5%  SFS 
treated  with  50  ng/mL  Dox.  After  24  hr  of  Dox  treatment,  media 
were  replaced  with  treated  media  containing  vehicle  (0.15% 
DMSO)  or  a  range  of  serially  diluted  ligands.  All  treatments  were 
conducted  in  the  presence  and  absence  of  100  nM  IC1  182,780. 
After  treatment  for  24  hr,  cells  were  washed  with  PBS  and  lysed 
with  35  p,L  lysis  buffer  (100  mM  K2HP04,  0.2%  Triton  X-100,  pH 
7.8).  Lysate  (30  p,L)  was  mixed  1:1  with  luciferase  substrate 
(Promega)  and  luminescence  was  measured  on  a  Victor  X5 
microplate  reader  (Perkin  Elmer,  Waltham,  MA)  using  lumines¬ 
cence  detection  and  a  700  nm  filter.  Total  protein  (5  p,L)  was 
quantified  using  BioRad  Protein  Assay  (BioRad).  EC50  values  were 
calculated  using  GraphPad  Prism  Software  (Version  5.04,  Graph- 
Pad  Software  Inc.,  San  Diego,  CA)  and  a  three  parameter  log  versus 
response  nonlinear  regression.  Two  tailed  t-tests  performed  with 
GraphPad  Prism  Software  were  used  to  determine  statistically 
significant  differences  from  control  treatments. 

2.5.  Gene  expression  analysis 

For  analysis  of  reporter  induction  by  cosmosiin,  Hs578T- 
ERaLuc  and  Hs578T-ER(3Luc  cells  were  split  in  phenol  red  free 
DMEM/F12  +  5%  SFS  and  treated  with  50  ng/mL  Dox  for  48  hr 
followed  by  treatment  with  DMSO  (0.1%),  1  nM  E2,  or  1  p,M 
cosmosiin  for  4  or  24  hr.  Total  RNA  was  extracted  using  RNEasy 
Plus  Kit  according  to  manufacturer  protocol  (Qiagen,  Valencia,  CA). 
RNA  (2  p,g)  was  reverse  transcribed  using  Superscript  II  RT 
according  to  manufacturer  protocol  (Invitrogen),  and  firefly 
luciferase  (FLuc)  expression  was  determined  by  reverse-transcrip¬ 
tion  polymerase  chain  reaction  using  primers  shown  in  Table  1. 

For  quantitative  real-time  PCR  analysis  of  endogenous  target 
gene  expression,  Hs578T-ERa  and  Hs578T-ER(3  cells  were  cultured 
in  phenol  red  free  DMEM/F12  +  10%  SFS  for  3  days  prior  to  the 
assay  to  remove  residual  estrogens  from  the  cells.  Cells  were  split 
in  phenol  red  free  DMEM/F12  +  5%  SFS  and  treated  with  50  ng/mL 
Dox  for  48  hr  prior  to  ligand  treatment.  Cells  were  treated  with  Dox 
and  ligands  or  vehicle  (0.1%  DMSO)  for  24  hr,  and  total  RNA  was 


Table  1 

Primer  and  probe  sequences. 


RPL13A 

Primer  1 
Primer  2 
Probe 

5'-TGT  TTG  ACG  GCA  TCC  CAC-3' 

5'-CTG  TCA  CTG  CCT  GGT  ACT  TC-3' 

5'-CTT  CAG  ACG  CAC  GAC  CTT  GAG  GG-3' 

C3 

Primer  1 
Primer  2 
Probe 

5'-AAC  TAC  ATC  ACA  GAG  CTG  CG-3' 

5'-AAG  TCC  TCA  ACG  TTC  CAC  AG-3' 

5'-CGT  TTC  CCG  AAG  TGA  GTT  CCC  AGA-3' 

JAG1 

Primer  1 
Primer  2 
Probe 

5'-GGA  CTA  TGA  GGG  CAA  GAA  CTG-3' 
5'-AAA  TAT  ACC  GCA  CCC  CTT  CAG-3' 

5'-TCA  CAC  CTG  AAA  GAC  CAC  TGC  CG-3' 

ITGA6 

Primer  1 
Primer  2 
Probe 

5'-ACC  CGA  GAA  GGA  AAT  CAA  GAC-3' 
5'-CGC  CAT  CTT  TTG  TGG  GAT  TC-3' 

5'-TGG  GTT  GGA  AGG  GCT  GTT  TGT  CA-3' 

FLuc 

Primer  1 
Primer  2 

5'-GGC  TGA  ATA  CAA  ACC  ATC  GG-3' 

5'  Ctt  TCT  TGC  TCA  CGA  ATA  CGA-3' 

extracted  using  RNEasy  Plus  Kit  according  to  manufacturer 
protocol  (Qiagen).  RNA  (2  pig)  was  reverse  transcribed  as  above, 
and  quantitative  PCR  was  performed  using  TaqMan  Prime  Time 
custom  designed  assays  (IDT,  Coralville,  IA),  FastStart  Universal 
Probe  Master  Mix  (Roche  Scientific),  and  a  CFX96  instrument 
(BioRad).  Primer  and  probe  sequences  are  shown  in  Table  1.  Data 
were  analyzed  using  the  AACq  method  calculated  by  the  CFX 
Manager  Software  (BioRad).  Two  tailed  t-tests  performed  with 
GraphPad  Prism  Software  were  used  to  determine  statistically 
significant  differences  from  control  treatments  using  data  from 
three  biological  replicates. 

2.6.  Cell  counting  assays 

Hs578T-ERa  and  Hs578T-ER(3  cells  were  cultured  in  phenol  red 
free  DMEM/F12  +  10%  SFS  for  3  days  prior  to  the  assay  to  remove 
residual  estrogens  from  the  cells.  Cells  were  seeded  at  a  density  of 
15,000  cells/well  in  phenol  red  free  DMEM/F12  +  5%  SFS  in 
triplicate  in  6  well  tissue  culture  dishes  in  the  presence  or  absence 
of  50  ng/mL  Dox.  After  24  hr,  the  cells  were  treated  with  DMSO 
(0.1%)  or  compound  in  the  presence  or  absence  of  50  ng/mL  Dox. 
Media  were  refreshed  every  48  hr,  and  cells  were  counted  after 
trypan  blue  exclusion  using  an  automated  cell  counter  (BioRad) 
according  to  manufacturer  protocol. 

3.  Results 

3.1.  Generation  of  Hs578T-ERaLuc  and  Hs578T-ERfiLuc  reporter  cell 
lines 

In  order  to  generate  stable  reporter  breast  cancer  cell  lines,  we 
first  cloned  a  construct  encoding  a  selection  marker  and  a  luciferase 
reporter  linked  to  EREs.  The  pGL4.32  vector  (Promega)  contains 
the  luc2P  gene  and  was  modified  to  contain  3  tandem  consensus 
EREs  upstream  of  the  minimal  promoter  (pGL4.3xERE,  Fig.  1  A).  Upon 
complete  sequencing,  the  estrogen  responsiveness  of  the  vector  was 
validated  in  ER-negative  HEK293  cells  transfected  with  full  length 
ERa  (Fig.  IB)  or  ER(3  (Fig.  1C).  The  pGL4.3xERE  reporter  showed 
extremely  low  background  with  a  65-fold  induction  in  cells 
transfected  with  ERct.  The  ER  antagonist  ICI  182,780  abolished 
estrogen  induced  expression,  reducing  the  luciferase  signal  to  that  of 
vehicle  treated  cells.  Cells  transfected  with  ER(3  showed  a  1 5-fold 
induction  of  luciferase  upon  E2  treatment;  ICI  182,780  inhibited 
luciferase  expression  in  both  vehicle  and  estrogen  treated  cells.  The 
minimal  background  luciferase  expression  and  the  selection  marker 
conferred  by  the  pGL4.3xERE  vector  made  the  vector  suitable  for 
creating  stable  reporter  cells  lines  for  the  identification  and 
characterization  of  ER  selective  agonists. 

In  order  to  create  stable  ER  reporter  breast  cancer  cell  lines,  an 
ER  negative  breast  cancer  cell  line  engineered  to  express  either 
ERa  or  ER(3  was  necessary.  Previously,  Secreto  and  coworkers 
created  such  lines  using  Hs578T  cells  [15],  a  triple  negative  breast 
cancer  cell  line  with  a  basal-like  gene  expression  profile  [16]. 
Hs578T  cells  lack  expression  of  ERa  and  ER(3  providing  a  clean 
background  in  which  to  express  ERa  or  ER(3.  Using  the  tetracycline 
inducible  system,  two  cell  lines  were  created  in  which  ERa  or 
ER[3  are  inducibly  expressed  (Hs578T-ERa  and  Hs578T-ER[3 
cells,  respectively)  [15].  Hs578T-ERa  and  Hs578T-ER[3  cells 
were  transfected  with  the  pGL4.3xERE  vector,  and  individual 
clones  were  isolated  after  hygromycin  selection.  Over  20  clones 
were  screened  for  estrogen  induced  luciferase  expression  (data  not 
shown).  One  clone  from  each  cell  line  was  selected  for  further 
characterization,  referred  to  here  as  Hs578T-ERaLuc  and  Hs578T- 
ERfBLuc.  Additional  ERa  and  ER(3  reporter  clones  were  used  to 
verify  reporter  data  obtained  from  Hs578T-ERaLuc  and  Hs578T- 
ER[3Luc  cells. 
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Fig.  1.  The  pGL4.3xERE  reporter  construct  is  estrogen  responsive.  (A)  Three  tandem 
EREs  were  inserted  upstream  of  the  luc2P  gene  in  the  pGL4.32  luciferase  reporter 
construct.  HEK293  cells  were  batch  transfected  with  the  pGL4.3xERE  reporter 
construct,  a  (3-galactosidase  construct,  and  full  length  ERa  (B)  or  ER(3  (C).  After 
allowing  24  hr  for  protein  expression,  cells  were  treated  in  triplicate  with  vehicle 
(DMSO)  or  1  nM  E2  and  vehicle  or  100  nM  ICI  182,780  (0.15%  final  DMSO 
concentration)  for  an  additional  24  hr.  Raw  luciferase  units  (RLUs)  were  normalized 
to  (3-galactosidase  to  normalize  for  transfection  efficiency.  Error  bars  represent 
standard  deviations. 


Hs578T-ERaLuc  and  Hs578T-ERpLuc  cells  were  first  character¬ 
ized  by  assessing  luciferase  induction  by  ER  ligands  in  the  presence 
or  absence  of  the  full  antagonist  ICI  182,780  (Fig.  2).  Cells  were 
treated  with  vehicle,  1  nM  E2, 1 0  nM  DPN  (a  reported  ER(3  selective 
agonist),  or  10  nM  PPT  (a  reported  ERa  selective  agonist).  PPT 
selectively  activated  luciferase  expression  in  Hs578T-ERaLuc, 
but  DPN  activated  the  reporter  in  both  Hs578T-ERaLuc  and 
Hs578T-ERflLuc  cells,  though  to  a  lesser  extent  in  Hs578T-ERaLuc 


Hs578T-ERal_uc 

A  ■  Veh  I  100  nM  ICI  182,780 


Hs578T-ERpLuc 


B  iVeh  i  100  nM  ICI  182,780 

•tetnnn,  * 


Fig.  2.  ER  subtype  selective  ligands  selectively  induce  luciferase  in  Hs578T-ERctLuc 
and  Hs578T-ER(3Luc  cells.  Hs578T-ERctLuc  (A)  and  Hs578T-ER|3Luc  (B)  cells  were 
seeded  in  triplicate  on  96  well  plates  in  the  presence  of  50  ng/mL  Dox  to  induce  ER 
expression.  After  24  hr,  cells  were  treated  with  vehicle  (DMSO),  1  nM  E2,  10  nM 
DPN,  or  10  nM  PPT  in  the  presence  or  absence  of  100  nM  ICI  182,780  (0.15%  final 
DMSO  concentration).  Cells  were  lysed  24  hr  after  ligand  treatment  and  raw 
luciferase  units  were  counted.  Error  bars  represent  standard  deviations.  *p 
values  <  0.05. 


cells.  Co-treatment  with  ICI  182,780  blocked  luciferase  induction 
in  both  cell  lines  (Fig.  2),  and  luciferase  was  not  induced  in  the 
absence  of  Dox  treatment  (data  not  shown). 

Basal  and  E2-induced  luciferase  signals  were  much  higher  in 
Hs578T-ERaLuc  cells  when  compared  to  Hs578T-ERPLuc  cells,  a 
trend  observed  in  all  luciferase  assays.  On  average,  Hs578T- 
ERPLuc  cells  expressed  630  luciferase  units  per  mg  protein  and 
Hs578T-ERaLuc  expressed  2900  luciferase  units  per  mg  protein  at 
saturating  E2  concentrations  (0.1  nM  or  greater).  A  range  of 
luciferase  signals  was  observed  among  the  clones  screened  (data 
not  shown),  suggesting  the  accessibility  of  the  reporter  in  the 
chromatin  may  be  responsible  for  differences  in  luciferase 
expression.  In  order  to  verify  Hs578T-ERaLuc  and  Hs578T-ER[3Luc 
cells  had  similar  ER  expression  levels  at  the  Dox  concentration 
used  throughout  the  study  (50  ng/mL),  quantitative  western  blots 
were  used  to  compare  ER  expression  in  the  parent  cell  lines  and 
reporter  cell  lines  (Fig.  3A).  Western  blots  with  FLAG  antibody 
demonstrated  similar  ER  expression  in  Hs578T-ERaLuc  and 
Hs578T-ERpLuc  cells  and  also  confirmed  expression  levels  similar 
to  the  parent  cell  lines.  In  addition,  whole  cell  ligand  binding  assays 
were  used  to  quantify  the  active  receptor  in  each  cell  line  (Fig.  3B). 
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Fig.  3.  Hs578T-ERaLuc  and  Hs578T-ER(3Luc  cells  express  similar  levels  of  ER.  (A) 
Quantitative  western  blot  with  Hs578T-ERa  (ERa),  Hs578T-ERaLuc  (ERaLuc), 
Hs578T-ERp  (ERpi),  and  Hs578T-ER(3Luc  (ERpLuc)  treated  with  vehicle  (-Dox)  or 
50  ng/mL  Dox  (+Dox).  ER  expression  was  detected  using  FLAG  antibody  and 
quantified  by  normalizing  to  beta-actin  using  the  Licor  Odyssey  near-infrared  gel 
reader.  The  normalized  integrated  intensity  for  the  FLAG  signal  is  shown  below  the 
images.  (B)  Ligand  binding  assays  confirmed  the  quantitative  western  blots. 
Hs578T-ERaLuc  and  Hs578T-ERpLuc  cells  were  seeded  in  triplicate  and  treated 
with  vehicle  or  50  ng/mL  Dox  for  24  hr.  Cells  were  labeled  with  20  nM  [3H]-E2  in  the 
presence  or  absence  of  cold  competitor  for  2  hr,  washed,  and  total  cell  lysate  was 
assessed  for  bound  radioactivity  as  described  in  Section  2.  MCF7  cells  were  included 
for  comparison.  Two  additional  wells  of  each  cell  line  and  condition  were  used  to 
determine  the  cell  number  and  the  numbers  of  receptors  per  cell  were  calculated 
based  on  a  1 : 1  molar  ratio  of  ligand  to  receptor.  The  average  and  standard  deviation 
of  three  independent  experiments  are  shown.  (C)  The  morphology  of  Hs578T- 
ERaLuc  and  Hs578T-ERpLuc  was  similar  to  that  of  the  parent  Hs578T-ERa  and 
Hs578T-ERp  cell  lines.  Representative  phase-contrast  microscopy  images  of  each 
cell  line  (lOOx  magnification). 

ERa  positive  MCF7  breast  cancer  cells  expressed  ~1 50,000 
receptors/cell  which  was  very  similar  to  reported  values  [17], 
Both  Hs578T-ERaLuc  and  Hs578T-ERPLuc  cells  expressed 
~120,000  receptors/cell  after  50  ng/mL  Dox  treatment.  The 


comparable  number  of  ERs  per  cell  suggests  that  differences  in 
ER  expression  do  not  account  for  the  higher  luciferase  signal 
observed  Hs578T-ERaLuc  cells.  Higher  luciferase  expression  in 
Hs578T-ERaLuc  cells  may  be  due  to  the  accessibility  of  the 
reporter  in  the  chromatin  or  the  enhanced  transcriptional  activity 
of  ERa,  in  agreement  with  previous  findings  that  the  transcrip¬ 
tional  activity  of  ERa  is  greater  than  that  of  ERP  on  ERE-containing 
reporters  [18].  Finally,  the  reporter  cell  lines  did  not  have  an 
altered  morphological  phenotype  compared  to  the  parent  cell  lines 
(Fig.  3C),  and  no  other  phenotypic  changes  due  to  the  integration 
of  the  luciferase  reporter  were  observed  in  Hs578T-ERaLuc  and 
Hs578T-ERfJLuc  cells. 

3.2.  Ligand  selectivity  of  Hs578T-ERaLuc  and  Hs578T-ER$Luc 
reporter  cell  lines 

We  next  assessed  ligand  subtype  selectivity  using  these 
isogenic  reporter  cell  lines.  All  luciferase  data  were  normalized 
to  the  luciferase  signal  induced  by  a  saturating  concentration  of  E2 
(0.1  nM)  and  expressed  as  the  percent  transactivation  relative  to 
0.1  nM  E2.  Dose-response  curves  were  obtained  for  E2,  DPN,  and 
PPT  to  characterize  the  sensitivity  of  the  reporter  cells  to  ER 
ligands  (Fig.  4).  Cells  were  treated  with  1 0-fold  dilutions  of  ligands 
and  approximate  EC50  concentrations  for  each  ligand  were 
calculated  from  3  independent  experiments  (Table  2).  The  ratios 
of  EC50  values  obtained  from  Hs578T-ERaLuc  cells  and  Hs578T- 
ER[3Luc  cells  are  also  presented  in  Table  2  and  provide  a  measure 
of  the  selectivity  of  the  ligands.  Higher  a/P  ratios  indicate 
selectivity  for  ERp. 

Both  cell  lines  were  highly  sensitive  to  estrogen  (Fig.  4A). 
Hs578T-ERaLuc  cells  showed  EC50  values  near  1  pM;  four 
additional  Hs578T-ERaLuc  clones  showed  similar  sensitivities 
(data  not  shown).  Hs578T-ER(3Luc  cells  also  showed  EC50  values 
for  estrogen  in  the  pM  range,  though  the  average  EC50  was  6.5- 
fold  higher  than  that  of  Hs578T-ERaLuc  cells.  Similar  differences 
in  estrogen  sensitivities  have  been  observed  in  other  ERE- 
luciferase  reporter  cell  lines  expressing  ERa  or  ERp  [12-14], 
suggesting  the  difference  in  E2  sensitivity  between  Hs578T- 
ERaLuc  and  Hs578T-ER(3Luc  cells  is  due  to  differences  in  the 
transactivation  of  ERa  and  ER(3. 

Next,  dose  responses  to  two  highly  selective  ERa  and  ERP 
agonists,  PPT  and  DPN  respectively,  were  analyzed  using  Hs578T- 
ERaLuc  and  Hs578T-ER[3Luc  cells.  PPT  showed  nearly  1000-fold 
selectivity  for  ERa  (Fig.  4B).  Surprisingly,  PPT  could  activate 
reporter  expression  in  Hs578T-ER[3Luc  cells  at  concentrations 
greater  than  100  nM,  although  it  could  not  induce  luciferase 
expression  to  the  same  extent  as  E2.  It  has  been  reported  that  PPT 
was  unable  to  induce  an  estrogen  responsive  reporter  in  HEC-1 
cells  transfected  with  ER(3  [10]  or  in  HELN-ERa  cells  [12].  DPN  was 
not  as  selective  as  PPT  and  could  maximally  activate  luciferase 
expression  Hs578T-ERaLuc  cells  at  100  nM  (Fig.  4C).  DPN  fully 
activated  ER(B  at  10  nM.  Though  DPN  has  been  shown  to  have  a  50 
to  70-fold  higher  binding  affinity  for  ERP  [12,19],  comparison  of 
EC50  values  showed  approximately  30-fold  selectivity  for  ERP  in 
these  reporter  assays. 

Next,  the  subtype  selectivity  of  two  natural  phytoestrogens, 
liquiritigenin  and  cosmosiin,  were  analyzed  using  Hs578T-ERaLuc 
and  Hs578T-ERPLuc  cells  (Fig.  5).  Liquiritigenin  is  a  phytoestrogen 
derived  from  Glycyrrhizae  uralensis  and  is  the  most  active 
estrogenic  component  of  MF101,  an  herbal  supplement  with 
therapeutic  potential  [11],  In  the  initial  characterization  of 
liquiritigenin,  Mersereau  and  coworkers  found  liquiritigenin 
showed  minimal  activation  of  ERa  at  concentrations  up  to 
2.5  p,M  in  transcriptional  assays  in  U20S,  HeLa,  or  WAR5  prostate 
cancer  cells  transfected  with  ERa  [11],  Binding  assays  demon¬ 
strated  that  liquiritigenin  had  a  20-fold  higher  affinity  for  ERp  and 
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Fig.  4.  Hs578T-ERaLuc  and  Hs578T-ER(3Luc  show  subtype  selective  activation. 
Dose  response  curves  of  E2  (A),  PPT  (B),  and  DPN  (C).  Hs578T-ERaLuc  and  Hs578T- 
ER(3Luc  were  seeded  in  triplicate  and  treated  with  50  ng/mL  Dox  for  24  hr.  Cells 
were  then  treated  with  a  range  of  ligand  concentrations  (0.15%  final  DMSO 
concentration)  for  24  hr.  Each  plate  contained  DMSO,  0.1  nM  E2,  and  100  nM  ICI 
182,780  for  controls.  Luciferase  signal  was  normalized  to  total  protein  in  each  well 
and  expressed  as  a  percent  transactivation  relative  to  signal  obtained  from 
saturating  E2  treatment  (0.1  nM).  Each  dose  response  experiment  was  conducted  at 
least  3  times;  data  shown  are  from  one  representative  experiment.  EC50  values  are 
shown  in  Table  2. 


Table  2 

Average  EC50  values  for  ER  ligands  (M  x  10_ 
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0.001 

Liquiritigenin 

100  (40) 

28  (2) 

3.6 

A 


Hs578T-ERaLuc 

Hs578T-ERpLuc 


B 


[Liq]  M 


[Cos]  M 


Time  (hr) 


Hs578TERaLuc 

VEH  .  E2  | COS 
4  24  4  24  4  24 


Hs578TER(?Luc 

VEH  |  E2  |  COS 
4  24  4  24  4  24 


FLuc  — - 

RPL13A  —  —  —  —  ■ 


Fig.  5.  Liquiritigenin  (Liq)  and  cosmosiin  (Cos)  induce  reporter  expression  in 
Hs578T-ERaLuc  and  Hs578T-ER(3Luc.  Dose  response  curves  of  liquiritigenin  (A)  and 
cosmosiin  (B).  Hs578T-ERaLuc  and  Hs578T-ER(3Luc  were  seeded  in  triplicate  and 
treated  with  50  ng/mL  Dox  for  24  hr.  Cells  were  then  treated  with  a  range  of  ligand 
concentrations  (0.15%  final  DMSO  concentration)  for  24  hr.  Each  plate  contained 
DMSO,  0.1  nM  E2,  and  100  nM  ICI  182,780  for  controls.  Luciferase  signal  was 
normalized  to  total  protein  in  each  well  and  expressed  as  a  percent  transactivation 
relative  to  signal  obtained  from  saturating  E2  treatment  (0.1  nM).  Each  dose 
response  experiment  was  conducted  at  least  3  times;  data  shown  are  from  one 
representative  experiment.  EC50  values  are  shown  in  Table  2.  EC50  values  for 
cosmosiin  could  not  be  determined  because  of  supramaximal  reporter  induction. 
The  supramaximal  induction  by  cosmosiin  was  not  due  to  supramaximal 
transcription  of  the  luciferase  reporter  (C).  Hs578T-ERaLuc  and  Hs578T-ER(3Luc 
cells  were  treated  with  50  ng/mL  Dox  for  48  hr  followed  by  treatment  with  DMSO 
(0.1%),  1  nM  E2,  or  1  pM  cosmosiin  for  4  or  24  hr.  Firefly  luciferase  (FLuc) 
expression  was  determined  by  RT-PCR.  RPL13A  expression  was  used  to  ensure 
equal  loading. 


selectivity  was  proposed  to  be  due  to  selective  recruitment  of 
co-activators  to  ER(3,  namely  SRC-2  [11].  Comparison  of  EC50 
values  showed  liquiritigenin  had  a  3.6-fold  selectivity  for  ER[3,  and 
maximal  reporter  induction  was  obtained  by  100  nM  liquiritigenin 
in  Hs578T-ER(3Luc  cells  and  1  p,M  in  Hs578T-ERaLuc  (Fig.  5A  and 
Table  2). 

Cosmosiin,  or  apigenin  7-glucoside,  is  a  flavone  found  in 
chamomile  [20]  that  was  identified  as  an  ER  agonist  that  selectively 
induces  ERct/[3  and  ER(3/(3  dimers  as  measured  by  bioluminescence 
resonance  energy  transfer  (BRET)  assays  (unpublished  data).  It  has  a 
3-fold  higher  binding  affinity  for  ER(3  as  measured  by  competitive 
ligand  binding  assays  (IC50  ERct  15.9  p,M,  IC50  ER(3  3.3  p,M, 
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unpublished  data).  Interestingly,  cosmosiin  induced  luciferase 
expression  to  a  much  greater  extent  than  E2,  an  effect  described 
as  supramaximal  induction  [21],  Even  at  concentrations  upto 
10  p,M,  cosmosiin  did  not  saturate  the  luciferase  output,  and  EC50 
values  could  not  be  reasonably  calculated  (Fig.  5B).  Another  Hs578T- 
ER(3Luc  clone  treated  with  cosmosiin  also  showed  supramaximal 
induction  (data  not  shown).  Cosmosiin  did  not  induce  luciferase 
expression  in  Dox-treated  cells  co-treated  with  1C1 182,780  or  cells 
not  treated  with  Dox  (data  not  shown),  suggesting  the  supramax¬ 
imal  induction  was  due  to  ERp  activation.  To  determine  if  the 
supramaximal  induction  truly  represented  enhanced  transcription¬ 
al  activation,  the  transcript  levels  of  luciferase  were  assessed  after  4 
and  24  hr  treatments  of  E2  and  cosmosiin  (Fig.  5C).  Cosmosiin  did 
not  induce  luciferase  expression  to  a  greater  extent  than  E2  in  either 
Hs578T-ERaLuc  or  Hs578T-EKpLuc  cells,  indicating  alternative 
mechanisms  are  responsible  for  the  supramaximal  effect. 

3.3.  Selective  regulation  of  ERa  and  ERP  target  genes  by  ERp  selective 
ligands 

We  next  sought  to  validate  the  subtype  selectivity  of  DPN,  PPT, 
liquiritigenin,  and  cosmosiin  by  assessing  regulation  of  endoge¬ 
nous  ER  target  genes.  Estrogen  responsive  target  genes  of  ERa  and 
ERP  were  previously  identified  in  Hs578T-ERa  and  Hs578T-ERp 
cells  [15],  and  two  ERP  target  genes  and  one  ERa  target  gene  were 
selected  for  analysis.  Cells  were  treated  with  50  ng/mL  Dox  for 
48  hr  to  induce  expression  of  the  receptors  and  further  treated 
with  the  corresponding  ligands  for  24  hr.  Complement  component 
3  (C3,  NM_000064)  was  up-regulated  in  Hs578T-ERP  cells  upon 
E2  treatment  (Fig.  6A).  DPN  and  liquiritigenin  were  capable  of 
inducing  C3  expression  to  a  comparable  level  as  E2  at  concentra¬ 
tions  that  fully  activate  ERP  with  minimal  ERa  activation,  as 
measured  by  reporter  assays  (Fig.  6A).  Cosmosiin  induced  C3 
expression  at  1  p,M,  but  not  to  the  same  extent  as  E2, 
demonstrating  cosmosiin  does  not  fully  activate  the  receptor  at 
this  concentration.  PPT  slightly  induced  C3  expression  compared  to 
DMSO  in  Hs578T-ERP  cells,  although  PPT  induced  expression  of  C3 
to  a  much  lesser  degree  compared  to  E2.  Repression  of  the  ERP 
target  gene  Jagged  1  (JAG1,  NM_000214)  occurred  to  a  similar 
degree  by  E2,  DPN,  liquiritigenin,  and  cosmosiin,  although  100  nM 
liquiritigenin  and  1  p.M  cosmosiin  do  not  fully  repress  JAGt 
expression  compared  to  E2,  DPN,  or  1  p,M  liquiritigenin  (Fig.  6B). 
Although  the  ERa  selective  agonist  PPT  slightly  induced  C3 
expression  in  Hs578T-ERp  cells,  it  had  no  effect  on  JAG 1  repression, 
demonstrating  incomplete  ERp  activation  by  PPT.  To  further 
validate  the  subtype  selectivity  observed  in  reporter  assays, 
expression  of  the  ERa  target  gene  alpha-6  integrin  ( ITGA6 , 
N1VL000210)  was  determined  after  treatment  of  Hs578T-ERa 
cells  with  E2,  DPN,  PPT,  liquiritigenin,  and  cosmosiin.  As  shown  in 
Fig.  6C,  ITGA6  was  up-regulated  by  E2  and  PPT  treatment,  but  DPN 
and  liquiritigenin  did  not  fully  activate  its  expression  at 
concentrations  that  showed  selectivity  in  reporter  assays  (10  nM 
and  100  nM,  respectively).  At  1  p,M,  liquiritigenin  and  cosmosiin 
were  capable  of  activating  ERa,  and  ITGA6  expression  was  induced 
in  Hs578T-ERa  cells.  Therefore,  the  subtype  selectivity  of  DPN  and 
liquiritigenin  observed  in  reporter  cell  lines  was  validated  by 
subtype  selective  regulation  of  endogenous  target  genes.  Cosmo¬ 
siin,  however,  activated  expression  of  an  Hs578T-ERa  endogenous 
gene  target  at  concentrations  that  only  slightly  activated  luciferase 
reporter  expression  in  Hs578T-ERaLuc  cells. 

3.4.  Growth  inhibition  of  Hs578T-ERfi  cells  by  liquiritigenin  and 
cosmosiin 

We  next  characterized  the  growth  effects  of  liquiritigenin  and 
cosmosiin  in  Hs578T-ERa  and  Hs578T-ERP  cells.  It  was  previously 
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Fig.  6.  ERP  selective  ligands  selectively  regulate  ER  target  genes.  Hs578T-ERa  and 
Hs578T-ERp  cells  were  treated  with  50  ng/mL  Dox  for  48  hr  to  induce  ER 
expression  followed  by  treatment  with  the  corresponding  ligands  for  24  hr.  Total 
RNA  was  assayed  for  expression  of  the  ERp  target  genes  C3  and  JAG  1  in  Hs578T-ERp 
cells  (A  and  B,  respectively)  and  the  ERa  target  gene  ITGA6  in  Hs578T-ERa  (C)  cells 
by  quantitative  reverse-transcription  polymerase  chain  reaction.  Target  gene 
expression  was  calculated  using  the  AACq  method  by  normalizing  to  the  ribosomal 
protein  RPU3A.  Data  represent  the  average  and  standard  deviation  of  three 
biological  replicates.  *p  values  <  0.05  compared  to  DMSO  control;  #p  values  <  0.05 
compared  to  E2  treatment. 


shown  that  E2  inhibits  the  growth  of  Hs578T-ER[3  cells  [15], 
supporting  the  notion  that  the  anti-proliferative  action  of  ER(3  may 
be  activated  by  estrogenic  ligands.  We  tested  whether  100  nM 
liquiritigenin,  a  concentration  at  which  ERP  was  selectively 
activated,  and  1  p,M  cosmosiin  could  also  inhibit  the  growth  of 
Hs578T-ERP  cells.  Hs578T-ERa  and  Hs578T-ER(3  cells  were 
treated  with  vehicle  (DMSO),  1  nM  E2,  100  nM  liquiritigenin  or 
1  p,M  cosmosiin  in  the  presence  or  absence  of  50  ng/mL  Dox  (with 
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Fig.  7.  Cosmosiin  (Cos)  and  liquiritigenin  (Liq)  inhibit  the  growth  of  Hs578T-ER(3  cells.  Hs578T-ERa  (A,  C,  and  E)  and  Hs578T-ER(3  cells  (B,  D,  F)  were  seeded  in  6-well  plates 
and  treated  with  vehicle  (A  and  B)  or  50  ng/mL  Dox  (C  and  D).  After  24  hr,  the  cells  were  treated  with  vehicle  (0.1%  DMSO)  or  the  indicated  ligands,  and  treatments  were 
refreshed  every  48  hr.  Cells  were  counted  at  the  times  indicated  using  trypan  blue  exclusion.  Comparisons  of  the  cell  number  on  day  5  are  represented  in  panels  E  (Hs578T- 
ERa)  and  F  (Hs578T-ER(3).  Data  represent  two  independent  experiments.  *p  values  <  0.05. 


or  without  ER,  respectively)  for  a  total  of  5  days.  When  ERa  and 
ERP  were  not  expressed  (-Dox),  the  compounds  had  no  effect  on 
the  growth  of  the  cells  (Fig.  7A  and  B).  In  contrast,  E2,  liquiritigenin, 
and  cosmosiin  inhibited  the  growth  of  Hs578T-ERf5  cells  when  ERp 
was  expressed  (+Dox,  Fig.  7D),  and  there  was  an  approximately 
50%  reduction  in  the  number  of  cells  after  5  days  of  treatment  with 
all  three  compounds  (Fig.  7F).  Hs578T-ERa  cells  showed  slight 
inhibition  with  E2  and  liquiritigenin  treatment  when  ERa  was 
expressed  (Fig.  7C),  but  there  was  not  a  statistically  significant 
effect  after  5  days  of  treatment  as  measured  by  2  independent 
experiments  (Fig.  7E).  However,  ERa  expression  in  ER  negative 
cells  often  leads  to  growth  inhibition  [22,23],  and  it  is  likely  that 
activation  of  ERa  inhibits  the  growth  of  Hs578T-ERa  cells.  This 
suggests  that  100  nM  liquiritigenin  partially  activates  ERa  despite 
minimal  regulation  of  ITGA6  at  this  concentration. 

4.  Discussion 

ERa  is  an  established  therapeutic  target  for  breast  cancer 
treatment,  but  the  development  of  subtype  selective  estrogenic 
ligands  has  gained  interest  with  the  identification  of  ERP  [1  ].  ERp 
opposes  the  actions  of  ERa  suggesting  that  it  may  be  a  potential 
therapeutic  target.  Exogenous  ERp  expression  in  ERa  positive 
breast  cancer  cells  impaired  E2  stimulated  proliferation  [24]  and 
tumor  growth  in  xenografts  [25].  In  support  of  the  anti¬ 
proliferative  role  of  ERp,  MCF7  cells  were  more  proliferative 
when  ERP  was  knocked  down  [6[.  Activation  of  ERP  by  subtype 
selective  ligands  may  enhance  ERP  growth  repression  without 
stimulating  proliferation  through  ERa;  indeed  ERP  selective 
ligands  inhibited  growth  of  HC11  mouse  mammary  cells  [5[.  Here, 
we  have  also  shown  that  ERp  ligands  can  inhibit  the  growth  of 
breast  cancer  cells  when  ERP  is  expressed.  In  breast  cancer, 
however,  ERP  expression  is  thought  to  decline  during  progression 
[26-28]  so  ligands  aimed  at  targeting  ERp  must  be  highly  selective 
and  used  only  in  patients  that  lack  ERa  or  those  with  low  ERa:ERP 
ratios  of  expression.  The  rate  of  ERp  positivity  in  breast  cancer  has 
been  reported  to  range  from  13%  to  83%  [29-32].  In  order  to 


effectively  target  ERP  for  cancer  treatment,  there  is  an  imminent 
need  to:  (a)  identify  ERP  selective  ligands  with  minimal  side 
effects  and  better  in  vivo  efficacy  and  selectivity,  and  (b)  design 
clinical  trials  to  recruit  patients  with  low  ERa: ERP  ratios  in  earlier 
stages  of  disease  progression. 

Although  ERP  selective  ligands  have  not  yet  been  used  for 
cancer  treatment,  the  therapeutic  value  of  ERp  has  been  assessed 
in  other  diseases.  Two  of  the  most  promising  ERp  selective 
therapies  are  the  ERP  selective  ligand  ERB-041  and  the  herbal 
extract  MF-101  [33].  Clinical  trials  have  been  completed  to 
determine  the  efficacy  of  ERB-041  for  treatment  of  Crohn’s 
disease,  endometriosis,  interstitial  cystitis,  and  rheumatoid 
arthritis.  Although  results  have  not  been  published  for  most  of 
the  clinical  trials,  results  of  the  rheumatoid  arthritis  trial  showed 
ERB-041  was  well  tolerated  but  did  not  improve  arthritis 
symptoms  [34].  MF-101  also  showed  a  relatively  safe  profile 
and  reduced  the  frequency  of  hot  flashes  in  a  phase  II  clinical  trial 
for  treatment  of  post-menopausal  symptoms  [35].  Liquiritigenin 
is  the  most  active  estrogenic  component  of  MF-101  [11], 
suggesting  ERP  selective  ligands  may  prove  useful  for  treating 
post-menopausal  symptoms. 

Strategies  to  identify  ER  subtype  selective  ligands  include 
competitive  ligand  binding,  dimerization,  transcriptional  reporter, 
and  proliferation  assays  [21,36].  Competitive  ligand  binding  assays 
provide  insight  into  binding  affinities  and  are  useful  for  high 
throughput  small  molecule  screening  [37],  but  they  are  limited 
because  ligands  can  act  as  agonists  or  antagonists  and  binding 
affinity  does  not  often  reflect  transcriptional  potency.  Subtype 
selective  ligands  have  been  identified  using  BRET  assays  that 
measure  receptor  dimerization  [38],  but  dimerization  assays 
cannot  differentiate  between  agonists  or  antagonists  [39].  Agonists 
can  be  characterized  using  MCF7  cell  proliferation  assays,  which 
are  highly  sensitive  and  provide  a  biologically  relevant  endpoint  in 
the  context  of  estrogen-sensitive  cells  [40].  However,  these  assays 
are  limited  by  a  lack  of  specificity,  as  non-estrogenic  mitogens  can 
stimulate  proliferation,  and  cannot  be  used  to  detect  subtype 
selective  agonists. 
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Transcriptional  assays  can  differentiate  between  agonists  and 
antagonists,  overcoming  limitations  of  binding  and  dimerization 
assays.  Mammalian  reporter  cell  lines  useful  for  identifying 
subtype  selective  ligands  have  been  created  from  HeLa  cervical 
carcinoma  cells  [12]  and  HEK293  kidney  cells  [13].  HELN-ERa  and 
HELN-ERJ3  were  generated  from  HeLa  cells  in  two  steps:  (1 )  stable 
integration  of  ERE-luciferase  to  generate  HELN  cells,  (2)  stable 
expression  of  ERa  or  ER(3  to  generate  HELN-ERa  and  HELN-ER(3 

[12] .  293/hERa  and  293/hER[3  cells  were  generated  by  a  similar 
strategy.  Only  one  breast  cancer  reporter  cell  line,  T47D-KBLuc,  is 
available  to  characterize  agonists  in  the  context  of  breast  cancer 
cells  [14],  but  both  ERa  and  ER(3  are  expressed,  preventing 
identification  of  subtype  selective  ligands. 

In  this  report,  we  described  the  development  of  a  new  set  of 
breast  cancer  reporter  cell  lines  to  characterize  subtype  selective 
estrogenic  ligands.  Hs578T-ERaLuc  and  Hs578T-ER(3Luc  cells  were 
highly  sensitive  to  E2  with  EC50  values  of  1  pM  and  6.5  pM, 
respectively  (Fig.  4A).  Similar  E2  sensitivity  was  observed  in  T47D- 
KBLuc  cells,  which  showed  an  approximate  EC50  of  3  pM  [14]. 
Hs578T-ERaLuc  and  Hs578T-ERpLuc  cells  were  more  sensitive  to  E2 
than  HELN-ER  and  293/ER  reporter  cells,  but  all  reporter  cell  lines 
showed  greater  E2  sensitivity  in  ERa  expressing  cells.  HELN-ERa 
cells  were  approximately  3  times  more  sensitive  to  E2  than  HELN- 
ER(3  cells  (EC50  of  0.017  nM  and  0.068  nM,  respectively)  [12]  and 
293/hERa  cells  were  approximately  4  times  more  sensitive  to  E2 
than  cells  expressing  ER(3  (EC50  of  50  pM  and  200  pM,  respectively) 

[13] .  Although  Hs578T-ERaLuc  and  Hs578T-ER[3Luc  cells  were  not 
created  using  the  same  strategy  as  HELN-ER  or  293/hER  reporter 
cells  and  likely  have  unique  genomic  integration  of  the  reporter, 
similar  sensitivities  observed  in  all  reporter  cell  lines  suggest  that 
this  does  not  inhibit  comparison  of  subtype  selectivity. 

Reporter  assays  with  two  ER  subtype  selective  ligands  confirmed 
that  Hs578T-ERaLuc  and  Hs578T-ER(3Luc  cells  could  be  used  to 
differentiate  between  ERa  and  ER(3  selective  ligands.  The  ER(3 
selective  agonist  DPN  maintained  33-fold  selectivity  in  Hs578T- 
ERLuc  cells  (EC50  of  0.26  nM  for  ER(3  and  8.5  nM  for  ERa,  Table  2). 
Dose  response  assays  with  the  ERa  selective  agonist  PPT  revealed 
the  sensitivity  of  Hs578T-ER[3Luc  cells  (Fig.  4B).  Although  PPT  was 
unable  to  activate  reporter  expression  in  HEC-1  cells  transfected 
with  ER(3  [12],  PPT  did  activate  reporter  expression  in  Hs578T- 
ER(3Luc  cells  at  high  concentrations,  although  not  to  the  full  extent 
induced  by  E2.  PPT  reporter  activation  was  blocked  by  1CI  182,780 
co-treatment  (Fig.  2A)  and  did  not  occur  in  the  absence  of  Dox 
treatment  (data  not  shown),  verifying  reporter  activation  was 
mediated  by  ER(3.  Despite  activation  of  ER(3  at  high  concentrations, 
PPT  could  not  fully  activate  reporter  expression  in  Hs578T-ER(3Luc 
cells  and  maintained  1000-fold  selectivity  for  ERa. 

Subtype  selectivity  of  two  natural  phytoestrogens,  cosmosiin 
and  liquiritigenin,  was  also  assessed  in  Hs578T-ERaLuc  and 
Hs578T-ER(3Luc  cells.  Liquiritigenin  maintained  selectivity  for 
ER[3  but  to  a  lesser  extent  than  expected,  as  it  has  been  shown  to 
minimally  activate  ERa  in  other  cell  lines  [11].  The  discrepancy  in 
the  selectivity  of  liquiritigenin  may  be  due  to  the  enhanced 
sensitivity  of  Hs578T-ERaLuc  cells,  differences  in  cofactor 
expression  in  Hs578T  cells,  or  purity  of  the  compound  (our  studies 
utilized  commercially  available  liquiritigenin  and  Mersereau  and 
coworkers  [11]  used  extract  from  G.  uralensis).  The  selectivity  of 
cosmosiin  could  not  be  assessed  using  luciferase  assays  due  to 
supramaximal  induction  (Fig.  5B).  Supramaximal  activation  of 
estrogen  responsive  reporters  have  been  described  in  many 
systems  [21].  Here,  we  showed  that  supramaximal  induction  by 
cosmosiin  was  not  due  to  enhanced  transcriptional  activation  of 
the  reporter  (Fig.  5C).  Despite  limitations  of  the  reporter  system, 
the  subtype  selectivity  of  cosmosiin  could  be  characterized  by 
assessing  target  gene  regulation  in  Hs578T-ERa  and  Hs578T-ER(3 
cells.  While  DPN  and  liquiritigenin  maintained  similar  extents  of 


selectivity  as  measured  by  reporter  assays,  cosmosiin  activated 
both  ERa  and  ER[3  as  measured  by  endogenous  gene  regulation 
(Fig.  6).  Cosmosiin  and  liquiritigenin  induced  similar  growth 
inhibitory  effects  as  E2  in  Hs578T-ER(3  cells,  indicating  the 
phytoestrogens  could  elicit  ERp  activation  to  a  similar  extent  as 
E2  (Fig.  7). 

Hs578T-ERaLuc  and  Hs578T-ER[3Luc  cells  have  several  advan¬ 
tages  for  identifying  ERP  selective  agonists  in  comparison  to 
available  mammalian  reporter  cell  lines.  First,  the  Hs578T  reporter 
cell  lines  have  inducible  expression  of  ERa  and  ERp,  allowing 
determination  of  off-target  reporter  activation  by  assessing 
reporter  expression  in  the  absence  of  Dox.  Second,  Hs578T-ERaLuc 
and  Hs578T-ERPLuc  cells  are  highly  sensitive  to  estrogenic  ligands. 
Third,  endogenous  gene  regulation  can  be  used  to  validate  subtype 
selectivity.  Finally,  growth  inhibition  assays  using  Hs578T-ERP 
cells  in  the  presence  and  absence  of  Dox  can  be  used  to  determine 
the  biological  endpoint  of  ERP  activation  and  validate  specificity  of 
ligands  to  ensure  they  do  not  have  off-target  cytotoxic  effects.  High 
throughput  screening  may  be  possible  using  Hs578T-ERaLuc  and 
Hs578T-ERPLuc  cells,  and  luciferase  assay  optimization  using 
Hs578T-ERPLuc  cells  has  shown  a  Z  factor  of  0.5  (data  not  shown), 
an  acceptable  range  for  high  throughput  screening  [41  ].  Therefore, 
Hs578T-ERaLuc  and  Hs578T-ERpLuc  cells  are  useful  for  the 
identification  and  characterization  of  ER  subtype  selective  ligands 
that  may  hold  therapeutic  promise. 
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Abstract 

Estrogens  play  essential  roles  in  the  progression  of  mammary  and  prostatic  diseases.  The  transcriptional  effects  of  estrogens 
are  transduced  by  two  estrogen  receptors,  ERa  and  ERfi,  which  elicit  opposing  roles  in  regulating  proliferation:  ERa  is 
proliferative  while  ERfi  is  anti-proliferative.  Exogenous  expression  of  ERfi  in  ERa-positive  cancer  cell  lines  inhibits  cell 
proliferation  in  response  to  estrogen  and  reduces  xenografted  tumor  growth  in  vivo,  suggesting  that  ER(3  might  oppose 
ERa's  proliferative  effects  via  formation  of  ERa/p  heterodimers.  Despite  biochemical  and  cellular  evidence  of  ERa/p 
heterodimer  formation  in  cells  co-expressing  both  receptors,  the  biological  roles  of  the  ERa/p  heterodimer  remain  to  be 
elucidated.  Here  we  report  the  identification  of  two  phytoestrogens  that  selectively  activate  ERa/p  heterodimers  at  specific 
concentrations  using  a  cell-based,  two-step  high  throughput  small  molecule  screen  for  ER  transcriptional  activity  and  ER 
dimer  selectivity.  Using  ERa/p  heterodimer-selective  ligands  at  defined  concentrations,  we  demonstrate  that  ERa/p 
heterodimers  are  growth  inhibitory  in  breast  and  prostate  cells  which  co-express  the  two  ER  isoforms.  Furthermore,  using 
Automated  Quantitative  Analysis  (AQUA)  to  examine  nuclear  expression  of  ERa  and  ERp  in  human  breast  tissue  microarrays, 
we  demonstrate  that  ERa  and  ERp  are  co-expressed  in  the  same  cells  in  breast  tumors.  The  co-expression  of  ERa  and  ERP  in 
the  same  cells  supports  the  possibility  of  ERa/p  heterodimer  formation  at  physio-  and  pathological  conditions,  further 
suggesting  that  targeting  ERa/p  heterodimers  might  be  a  novel  therapeutic  approach  to  the  treatment  of  cancers  which  co¬ 
express  ERa  and  ERp. 
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Introduction 

Estrogens  exert  their  biological  effects  via  interaction  with  two 
estrogen  receptors  (ERs),  ERa  and  ERfi  [1,2].  ERs  regulate  key 
physiological  functions  in  the  reproductive  tract,  breast,  prostate, 
bone,  brain  and  the  cardiovascular  system  [1,2].  In  some  organs, 
ERa  and  ERfi  are  expressed  at  similar  levels  but  in  different  cell 
types  [3].  For  example,  in  the  prostate,  ERa  is  predominately 
expressed  in  stroma  while  ERfi  is  expressed  in  the  epithelium. 
Both  receptors  are  expressed  in  normal  mammary  epithelial  cells 
[4],  Studies  with  ERa  knockout  mice  (aERKO)  demonstrate 
that  ERa  is  essential  for  ductal  formation  and  mammary  gland 
development  [5].  Although  ERfi  knockout  mice  (fiERKO) 
generate  mild  mammary  phenotypes,  Ki-67  expression  is 
increased  in  luminal  mammary  epithelial  cells  of  fiERKO  mice 
[6],  suggesting  that  ERfi  may  be  important  for  terminal 
differentiation  of  mammary  epithelial  cells.  ERa  and  ERfi  are 
also  involved  in  growth  and  differentiation  of  the  prostate  gland 


and  progression  of  prostate  disease  [7,8] .  A  recent  study  showed 
that  stromal  ERa  promotes  prostatic  carcinogenesis  [9].  More¬ 
over,  hyperplasia  was  observed  in  the  prostates  of  fiERKO  mice 
[10]  and  ERfi  expression  was  silenced  in  a  subset  of  malignant 
human  breast  and  prostate  cancers  [11,12],  suggesting  that  ERfi 
plays  protective  roles  in  diese  diseases. 

The  classic  mechanism  through  which  the  ERs  modulate  gene 
expression  is  a  cascade  of  events:  ligand  binding  to  ERa  or  ERfi 
induces  receptor  dimerization,  either  as  homodimers  (ERa/ERa 
or  ERfS/ERf!)  or  heterodimers  (ERa/ERfi),  translocation  of 
dimers  to  the  nucleus,  and  recognition  of  Estrogen  Response 
Elements  (EREs)  on  DNA.  The  target  genes  activated  by  these 
events,  and  hence  the  physiological  responses,  depend  on  the 
dimer  pair  activated  by  the  ligand.  Indeed,  several  studies  have 
shown  that  ERa  and  ERfi  exhibit  opposing  roles  in  cellular 
proliferation  and  apoptosis,  with  ERa  inducing  the  transcription  of 
pro-proliferative  and  anti-apoptotic  target  genes,  and  ERfi  being 
anti-proliferative  and  pro-apoptotic  [13,14,15],  In  accordance 
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with  this  notion,  target  gene  studies  reveal  that  ERa  and  ER(3  may 
have  distinct  biological  functions;  it  is  believed  that  ERa  promotes 
cell  growth,  while  ERfl  inhibits  it  in  breast  and  prostate  cancer 
cells  [1 1,14,16, 17,18,19].  It  has  thus  been  deduced  that  the  role  of 
the  ERa/ a  homodimer  is  to  accelerate  cellular  proliferation,  thus 
lending  to  carcinogenesis  and  tumor  progression,  while  conversely 
the  transcriptional  activation  from  ERp/p  homodimers  is  diought 
to  be  protective  against  hormone-dependent  diseases  including 
breast  and  prostate  cancers  [13,14,15]. 

ERP  has  well  known  growth  modulatory  activity  in  ERa- 
positive  breast  cancer  cells.  Compared  with  tumors  expressing 
ERa  alone,  the  co-expression  of  ERP  has  been  correlated  with  a 
more  favorable  prognosis  [20]  and  decreased  biological  aggres¬ 
siveness  [11,21,22,23,24],  Moreover,  ERp  has  been  shown  to 
modulate  the  proliferative  actions  of  estrogens  when  co-expressed 
with  ERa  [13,19,25,26]  and  can  be  considered  an  endogenous 
partial  dominant  negative  receptor  [27,28],  ERP  is  thought  to 
counteract  the  stimulatory  effects  of  ERa  through  heterodimer¬ 
ization  of  the  two  receptors  [29,30].  Indeed,  these  heterodimers 
have  been  shown  to  form  and  maintain  function  [31],  and  they 
have  been  suggested  to  be  responsible  for  the  activation  of  target 
genes  which  are  distinct  from  those  induced  by  either  homodimer 
[32,33].  The  co-expression  of  ERP  with  ERa  results  in  reduced 
ERa-mediated  proliferation  and  invasion  of  breast  cancer  cells 
[1 1,16,17,18, 19],  at  least  in  part  due  to  ERP's  inhibition  of  ERa 
selective  target  gene  expression.  Furthermore,  in  the  ERa/ERP- 
positive  mouse  mammary  epithelial  cell  line  HC 1 1 ,  ERa  drives 
cellular  proliferation  whereas  ERP  contributes  to  growth  inhibi¬ 
tion  and  apoptosis  in  response  to  1 7  P-estradiol;  (E2);  the  loss  of 
ERP  in  diis  cell  line  results  in  cellular  transformation  [14].  Thus, 
the  ERa:ERP  ratio  determines  whether  E2  will  induce  cellular 
proliferation.  Despite  the  fact  that  the  ERa/p  heterodimer  has 
been  proposed  to  have  a  biological  role  that  is  unique  from  that  of 
either  homodimer,  the  biological  function  of  these  heterodimers  in 
vivo  has  until  now  remained  elusive,  at  least  in  part  due  to  the 
heterogeneous  population  of  dimers  existent  upon  the  co¬ 
expression  of  ERa  and  ERP  and  the  lack  of  heterodimer-specific 
compounds  to  elucidate  their  functions. 

To  circumvent  this  issue,  the  identification  of  ERa/p 
heterodimer-selective  ligands  that  activate  the  transcriptional 
effects  of  ERa/ P  heterodimers,  but  not  that  of  either  homodimer, 
were  sought  in  order  to  shed  light  upon  the  transcriptional 
outcomes  and  biological  roles  of  these  heterodimers.  To  this  end,  a 
multi-step  high  throughput  small  molecule  screen  for  ER 
transcriptional  activation  and  dimer  selectivity  was  developed 
(Figure  1).  This  screening  resulted  in  the  identification  of  two 
phytoestrogens  that  are  transcriptionally  active  and  ERa/p 
heterodimer-selective  at  specific  concentrations.  These  compounds 
were  rigorously  characterized  for  their  biological  activity  in  cell- 
based  assays  (Figure  1).  The  results  of  these  studies  suggest  that  the 
ERa/p  heterodimer  exerts  growth  inhibitory  effects  in  breast  and 
prostate  epithelial  cells.  These  compounds  may  setve  not  only  as 
tools  for  deciphering  the  biological  functions  of  the  ERa/p 
heterodimer,  but  also  potentially  as  a  means  for  therapeutically 
targeting  ERa/p  heterodimers  in  hormone-dependent  diseases 
including  breast  and  prostate  cancers. 

Results 

Characterization  of  Lead  Compounds  Cosmosiin  and 
Angolensin  Using  Bioluminescence  Resonance  Energy 
Transfer  (BRET)  and  Reporter  Assays 

We  developed  two-step  high  throughput  screening  (HTS)  for 
identification  of  ER  dimer-selective  ligands  (unpublished).  The 


primary  screening  and  counter-screening  in  the  presence  of  the 
antagonist  ICI  182,780  (Fulvestrant)  for  ER-specific  transcription¬ 
al  activity  was  performed  in  T47D-KBLuc  as  described  in  the 
Methods  section.  ER  dimer  selectivity  of  the  primary  hits  was 
assessed  in  secondary  HTS  BRET  assays  as  described  in  the 
Methods  section  and  in  [34],  Several  compounds  with  dimer 
selectivity  were  identified  after  performing  two-step  HTS  on 
>5200  compounds  at  the  UWCCC  Small  Molecule  Screening 
Facility  (unpublished  results).  Two  phytoestrogens,  cosmosiin 
(apigenin-7-glucoside)  and  angolensin  (R)  (Fig.  2),  were  identified 
in  HTS  as  ER  dimer  selective  ligands.  Angolensin  exists  in  two 
enantiomeric  forms;  only  the  R  form  was  identified  and  used  in 
this  study  and  is  thus  abbreviated  as  angolensin  hereafter.  To 
determine  if  they  bind  the  same  ligand  binding  pocket  as 
1 7  P-estradiol  and  to  measure  their  binding  affinity  to  recom¬ 
binant  ERs,  we  employed  in  vitro  Fluorescence  Polarization  (FP) 
competition  binding  assays  [35].  The  IC50  values  for  cosmosiin 
binding  to  ERa  and  ERP  were  15.9  pM  and  3.3  pM,  respectively 
(Fig.  2A).  The  IC50  values  for  angolensin  binding  to  ERa  and  ERP 
were  2.2  pM  and  4.7  pM,  respectively  (Fig.  2B). 

The  ER  dimer  selectivity  was  validated  in  BRET  and  reporter 
assays  in  ER-negative  HEK293  cells  as  described  [35].  While 
cosmosiin  exhibits  preference  for  inducing  both  ERp/p  homodi- 
mers  and  ERa/p  heterodimers  (Fig.  3A),  angolensin  exhibits 
ERa/p  heterodimer  selectivity  (Fig.  3B).  Neither  compound  shows 
preference  for  inducing  ERa/ a  homodimers.  Because  the  lower 
limit  of  detection  for  these  compounds  was  1  pM,  concentrations 
lower  than  1  pM  are  not  shown  in  this  figure,  although  they  were 
tested  in  a  range  from  1  nM  to  10  pM;  below  1  pM,  the  BRET 
ratios  were  the  same  as  vehicle-treated.  Furthermore,  the  ability  of 
these  lead  compounds  to  induce  the  transcriptional  activity  of  ERa 
alone,  ERP  alone,  or  ERa  in  combination  with  ERP  was  tested  at 
a  range  of  concentrations  using  the  HEK293  ERE-luciferase 
reporter  assays  (Fig.  3C  and  3D).  Although  these  reporter  assays 
do  not  directly  examine  ERa/p  heterodimerization,  the  condition 
in  which  ERa  and  ERP  are  cotransfected  can  be  compared  with 
each  receptor  transfected  alone. 

As  shown  in  Figure  3B,  BRET  assays  reveal  that  angolensin  is 
capable  of  efficiently  inducing  the  formation  of  ERa/  P  heterodimers 
at  1  pM  and  10  pM,  while  not  inducing  ERa/a  or  ERp/p 
homodimers.  ERa/p  heterodimerization  appears  to  be  favored  in 
the  presence  of  angolensin  ,  and  in  the  condition  in  which  ERa 
and  ERP  are  co-transfected  for  luciferase  reporter  assays,  the 
highest  fold  induction  of  transcriptional  activity  relative  to  DMSO 
vehicle  is  obseived  (Fig.  3D).  Thus,  angolensin  (R)  appears  to  be 
an  ERa/p  heterodimer-selective  ligand  at  10  pM.  Cosmosiin 
appears  to  be  less  selective  in  terms  of  its  ability  to  induce  ERa/  P 
heterodimers,  as  ERp/p  homodimers  are  also  induced  in  BRET 
assays;  however,  ERa/ a  homodimers  are  not  induced  by 
cosmosiin  (Fig.  3A).  Cosmosiin  at  1  pM  appears  to  transcription¬ 
ally  activate  ERp/p  homodimers  and  ERa/p  heterodimers 
(Fig.  3C).  At  10  pM  cosmosiin,  while  ERa/a  and  ERp/p 
homodimers  were  slightly  activated,  co-transfecting  ERp  with 
ERa  exhibited  much  stronger  transcriptional  activity  (Fig.  3C). 
Thus,  cosmosiin  appears  to  be  ERp/p  homodimer-  and  ERa/p 
heterodimer-selective  at  1  pM. 

The  transcriptional  activity  of  ERa/ a  homodimers  treated  with 
1 0  pM  cosmosiin  is  despite  the  finding  that  the  BRET  assay  does 
not  show  statistically  significant  ERa/a  homodimerization 
(Fig.  3A).  The  most  likely  explanation  for  this  discrepancy  is 
differences  in  sensitivity  between  BRET  and  the  luciferase  reporter 
assays.  These  BRET  assays  and  luciferase  reporter  assays  are 
performed  under  different  conditions  and  measure  different  signal 
outputs:  BRET  captures  a  single  moment  in  time  in  which  ERa 
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Figure  1 .  Flow  scheme  of  high  throughput  screening  and  characterization  of  compounds  with  selectivity  for  ERa/ERp  heterodimers. 

A  library  of  >5200  small  molecules  was  screened  ER  transcriptional  activity  using  T47D-KBLuc  cells.  Molecules  with  transcriptional  activity  were  then 
screened  for  ERot/ot,  ERa/(S,  or  ER  p/p  dimerization  potential  using  BRET  assays.  Two  phytoestrogens,  angolensin  and  cosmosiin,  were  identified  as  ER 
dimer  selective  ligands.  These  molecules  were  validated  using  in  vitro  binding  assays  and  BRET  and  ERE-luciferase  reporter  assays.  Heterodimer 
selective  concentrations  were  identified  as  10  pM  angolensin  and  1  pM  cosmosiin.  The  cellular  effects  of  these  two  heterodimer-selective 
concentrations  were  characterized  using  cell  migration  and  proliferation  assays. 
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and  ERP  may  or  may  not  be  dimerized.  This  moment  in  time  was 
observed  after  1  hour  incubation  with  ligand.  Conversely,  the 
luciferase  reporter  assay  measures  an  accumulation  of  transcrip¬ 
tional  output  signal  (the  transcribed  luciferase  protein)  over  18— 
24  hours.  Consequently,  the  dimerization  ratios  obtained  via  the 
BRET  assay  do  not  always  completely  agree  with  the  transcrip¬ 
tional  profiles  obtained  in  the  luciferase  reporter  assays  for  a  given 
ligand.  Therefore,  it  is  important  to  consider  the  direct  dimerization 
of  ERa  and  ERp  in  conjunction  with  the  transcriptional  output  of 
these  diverse  dimer  pairs. 

Selection  and  generation  of  cell  lines  expressing  different 
amounts  of  ERa  and  ERp 

In  order  to  characterize  the  cellular  effects  of  cosmosiin  and 
angolensin,  we  surveyed  a  variety  of  breast  and  prostate  cell  lines 
for  co-expression  of  ERa  and  ERp.  As  shown  in  Fig.  4,  the  non- 
tumorigenic  mammary  epithelial  cell  HC 1 1  and  prostate  cancer 
cell  line  PC  3  were  found  to  express  both  receptors  (Lanes  1  and  2) 
as  reported  by  others  [14,36];  in  contrast,  DU-145  expresses  only 
ERp  [36]  (lane  6)  and  MDA-MB-231  is  negative  for  both  ERa 
and  ERP  (lane  5).  To  delineate  the  functions  of  ERa/p 
heterodimers,  we  knocked  down  ERa  and  ERP  transcript  levels 


in  PC  3  cells  by  means  of  stable  transfection  with  specific  shRNA 
plasmids  targeting  ERa  and  ERP,  respectively.  Western  blotting 
results  showed  that  ERa  is  selectively  silenced  in  PC3-shERa  cells 
and  ERP  is  selectively  silenced  in  PC3-shERP  cells  (Fig.  4A,  lanes 
3  and  4).  The  silencing  of  one  ER  did  not  influence  the  expression 
of  the  other.  All  of  these  characterized  cell  lines  were  subsequently 
used  for  determination  of  compounds’  cellular  effects. 

Cosmosiin  and  angolensin  inhibit  cell  motility  and 
migration  but  not  apoptosis  in  PC3 

In  order  to  examine  the  influences  of  these  ERa/  P  heterodimer¬ 
activating  compounds  on  cell  migration,  wound  healing  assays 
were  employed  using  migratory  PC3  cells.  This  assay  gives  a 
qualitative  measure  of  a  compound’s  ability  to  inhibit  cell 
migration.  For  these  assays,  1  pM  cosmosiin  and  10  pM 
angolensin  were  utilized  because  these  are  the  concentrations  at 
which  ERa/  P  heterodimers  are  most  highly  selectively  induced  by 
each  respective  compound.  As  shown  in  Figure  5A,  the  vehicle 
DMSO  (0.1%)  was  unable  to  inhibit  the  migration  of  PC3  cells  in 
scratch  wound  healing  assays:  cells  can  be  seen  infiltrating  the 
wound  24  hours  after  scraping,  and  the  wounds  are  completely 
filled  72  hours  after  scraping.  Conversely,  both  10  pM  angolensin 
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Figure  2.  Fluorescence  polarization  competition  binding  assays  for  ERa  and  ER0.  Cosmosiin  (A)  and  angolensin  (B)  bind  to  recombinant 

ERa  and  ER|3  with  jaM  affinities. 
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and  1  pM  cosmosiin  are  able  to  inhibit  the  ability  of  PC  3  cells  to 
infiltrate  the  wounds,  indicating  that  these  compounds  can  hinder 
cell  motility. 

To  quantitatively  measure  the  ability  of  cosmosiin  and  angolensin 
to  inhibit  cell  migration,  transwell  assays  were  employed.  Figure  5C 
shows  that  10  pM  angolensin  can  inhibit  the  ability  of  PC  3  cells  to 
migrate  through  the  pore,  and  this  inhibition  of  migration  is  ablated 
by  the  ER  antagonist  ICI  182,780.  1  pM  angolensin,  a  concentra¬ 
tion  at  which  ERa/ P  heterodimers  are  not  transcriptionally  active 
(Fig.  3D),  has  a  negligible  effect  on  cell  migration.  Both  1  pM  and 
10  pM  cosmosiin  can  inhibit  cell  migration  through  the  pore,  and 
this  inhibition  of  migration  is  ablated  by  ICI  182,780  (Fig.  5B). 
These  results  are  recapitulated  when  the  transwell  is  coated  with 
matrigel  (data  not  shown),  indicating  that  in  addition  to  dampening 
the  ability  of  PC  3  cells  to  migrate,  these  compounds  are  able  to 
dampen  the  ability  of  PC3  cells  to  invade. 

The  abilities  of  these  lead  compounds  to  influence  apoptosis  in 
PC3  cells  were  next  evaluated  using  caspase  3/7  assays.  PC3  cells 
were  incubated  with  the  indicated  concentrations  of  DMSO 
vehicle  (0.1%),  the  indicated  concentrations  of  cosmosiin  or 
angolensin  (Fig.  SI  A  and  SIB),  or  the  positive  control  cisplatin 
(10  pg/mL)  for  24,  48,  and  72  hours.  Cisplatin  did  not  activate  the 
caspases  3/7  pathway  at  24  hours  and  48  hours  (data  not  shown); 
only  at  72  hours  was  a  weak  induction  of  the  caspases  3/7 


observed  (Fig.  SI).  At  no  time  point  did  these  compounds  reveal 
any  activation  of  the  caspase  3/7  pathway.  Thus,  it  appears  that 
cosmosiin  and  angolensin  are  not  strong  inducers  of  apoptosis,  at 
least  through  the  caspase  3/7  pathway. 

Determination  of  the  growth  effects  of  compounds  in 
PC3,  PC3-shERa,  PC3-shERp  cells 

To  determine  if  these  compounds  also  inhibit  cell  proliferation 
in  addition  to  migration,  MTT  assays  were  employed.  This  assay 
measures  mitochondrial  activity  when  yellow  MTT  (3-(4,5- 
Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium  bromide)  is  reduced 
to  its  purple  formazan  metabolic  product  [37],  Thus,  the  ability  of 
a  cell  to  metabolize  MTT  to  formazan  is  correlated  to  its 
metabolic  activity  and  cellular  growth.  To  show  that  PC 3  cells 
express  functional  ERs  and  that  E2’s  cellular  effects  are  ER- 
dependent,  we  compared  E2’s  growth  effects  in  PC3,  PC3-shERa, 
PC3-shER(3  cell  lines.  As  shown  in  other  ERa  and  ERP  co¬ 
expressing  cell  lines  [14],  E2  exhibits  no  effects  in  proliferation  of 
PC3  (Fig.  S2A).  However,  when  ER|3  expression  was  blocked,  E2 
induced  proliferation  (Fig.  S2C)  and  E2’s  proliferative  effects  were 
completely  abrogated  by  the  pure  ER  antagonist  ICI  182,780  and 
the  ERa  selective  antagonist  MPP  dihydrochloride  (Fig.  S2C, 
middle  and  right  panels).  This  result  recapitulates  the  previous 
finding  in  HC 1 1  mammary  epithelial  cells  that  ERa  drives 
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Figure  3.  The  dimer  selectivity  for  cosmosiin  and  angolensin  was  demonstrated  in  dose-response  BRET  assays  (A  and  B)  and 
reporter  assays  (C  and  D)  in  HEK293  cells.  ER  dimer-specific  BRET  assays  were  performed  over  a  range  of  compound  concentrations  of 
cosmosiin  (A)  and  angolensin  (B).  HEK293T  ERE-luciferase  transcriptional  assays  reveal  each  compound's  ability  to  transcriptionally  activate  various 
dimer  pairs  (C  and  D).  ERa  alone,  ER(3  alone,  or  ERa+ER(3  was  transfected  along  with  an  ERE-luciferase  element  in  order  to  test  the  ability  of  cosmosiin 
(C)  and  angolensin  (D)  to  transcriptionally  activate  these  various  ER  dimer  pairs.  RLU,  relative  luciferase  units.  Error  bars  represent  standard  deviations 
from  the  mean  of  triplicate  samples.  In  BRET  (A),  p  values  indicate  all  pairs  with  statistical  significance  by  the  Student's  T-Test. 
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proliferation  in  response  to  E2  [14].  It  appears  that  silencing  ER(3 
in  PC  3  cells  causes  the  cells  to  respond  to  E2  with  increased 
proliferation,  similar  to  breast  cancer  cells  [19,28].  In  contrast  to 
HC 1 1  where  ERfl  is  growth  inhibitory,  knockdown  of  ERa  did 
not  result  in  E2-dependent  growth  inhibition  (Fig.  S2B).  The 
discrepancy  might  be  due  to  cell  line  specific  effects. 

The  cellular  effects  of  cosmosiin  and  angolensin  were  deter¬ 
mined  in  PC3,  PC3-shERa  or  PC3-shER|3  cells  at  concentrations 
that  display  ER  selectivity.  As  shown  in  Figure  6,  both  1  pM  and 
10  pM  cosmosiin  (Fig.  6A)  and  10  pM  angolensin  (Fig.  6B)  were 
able  to  inhibit  the  growth  of  PC  3  cells  compared  to  the  vehicle 
DMSO.  The  inhibition  of  growdi  due  to  10  pM  angolensin  was 
ablated  by  the  antagonists  ICI  182,780,  indicating  that  this 
response  is  ER-specific.  The  inhibition  of  growth  by  1  pM 
cosmosiin  was  also  ER-specific.  However,  the  inhibition  of  growth 
by  10  pM  cosmosiin  was  not  ablated  by  the  antagonist,  indicating 
that  this  response  is  not  ER-specific  in  PC  3  cells  and  is  likely  due  to 
off-target  effects  or  non-genomic  ER  signaling.  Cell  counting  and 
viability  assays  with  Trypan  blue  staining  ruled  out  the  possibility 
of  general  cytotoxicity  due  to  these  compounds  (Fig.  S3A,  S3B). 


In  PC3-shERa  cells,  ERfl  is  the  only  functional  ER  present; 
thus,  ERfl/p  homodimers  are  the  only  ER  dimers  capable  of 
forming  and  activating  transcription.  The  growth  inhibition 
obsetved  by  1  pM  cosmosiin  (Fig.  6C)  and  10  pM  angolensin 
(Fig.  6D)  in  the  parent  PC  3  cells  is  ablated  with  the  loss  of  ERa  in 
PC3-shERa  cells.  The  addition  of  the  ER  antagonist  ICI  182,780 
had  no  effect  on  this  cell  line  in  the  presence  of  1  pM  cosmosiin 
and  10  pM  angolensin  compared  to  these  ligands  alone.  However, 
10  pM  cosmosiin  was  still  able  to  inhibit  the  growth  of  these  PC3- 
shERa  cells  in  both  the  absence  and  presence  of  the  ER  antagonist 
ICI  182,780  (Fig.  6C).  Thus,  10  pM  cosmosiin  is  confirmed  to 
have  off-target,  ER  non-specific  influences  on  growth  regulation. 

In  PC3-shERP  cells,  ERa  is  the  only  functional  ER  present; 
thus,  ERa/ a  homodimers  are  the  only  ER  dimers  capable  of 
forming  and  activating  transcription.  As  shown  in  Figure  6F, 
angolensin  has  a  negligible  effect  in  this  cell  line,  and  treatment 
with  the  ER  antagonist  ICI  182,780  completely  ablates  any 
growth  effects  observed  in  the  presence  of  this  compound.  This 
finding  is  consistent  with  angolensin’s  high  degree  of  ERa/p 
heterodimer  selectivity.  However,  contrary  to  observations  in  PC  3 
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Figure  4.  Determination  of  relative  expression  levels  of  ERa 

and  ERp  in  various  cell  lines.  Western  blotting  analyses  of  ERor  and 
ERp  expression  in  HC11  (lane  1),  PC3  (lane  2),  PC3shERa  (lane  3), 
PC3shERp  (lane  4),  MDA-MB-231  (lane  5),  and  DU145  (lane  6). 
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cells  and  PC3shERa  cells,  cosmosiin  increases  the  growth  of  PC3- 
shERp  cells  at  both  1  |tM  and  10  pM  (Fig.  6E).  The 
transcriptional  activation  of  ERa/ a  homodimers  is  induced  with 
10  pM  cosmosiin  in  HEK293  ERE-luciferase  assays  (Fig.  3C), 
which  is  in  keeping  with  its  ability  to  increase  the  growth  of  PC  3- 
shERp  cells  at  this  concentration.  However,  the  increase  in 
growth  due  to  1  pM  cosmosiin  is  not  predicted  by  the  HEK293 
ERE-luciferase  assay  (Fig.  3C).  These  data  were  confirmed  with 
cell  counting  and  viability  assays  with  Trypan  blue  staining  (data 
not  shown).  These  growth  increases  in  PC3-shERP  cells  due  to 
1  pM  cosmosiin  are  ablated  by  the  antagonist  ICI  182,780 
(Fig.  6E),  suggesting  that  these  growth  increases  are  due  to  ERa/ 
a  homodimers.  Intriguingly,  treatment  with  diese  antagonists  in 
the  presence  of  10  pM  cosmosiin  not  only  ablates  the  growth 
increases  observed  at  this  concentration,  but  actually  results  in 
decreased  growth  (Fig.  6E).  This  inhibited  growth  in  PC3-shERP 
cells  when  ERa  is  antagonized  may  be  explained  by  the  off-target 
effects  mediated  by  this  compound:  when  ERa  is  the  only  ER 
present,  it  is  not  damped  by  heterodimerization  with  ERP  and  is 
instead  able  to  bind  cosmosiin  to  increase  cellular  growth; 
however,  when  ERa  is  antagonized  in  this  cell  line,  cosmosiin  is 
free  to  mediate  its  off-target  growth  inhibitory  effects,  resulting  in 
decreased  growth. 

The  ERa/p  heterodimers  were  found  to  be  growth 
inhibitory  using  PC3  derived  cell  lines  and  ERa/p  heterodi¬ 
mer-selective  ligands  at  concentrations  determined  to  be 
heterodimer-selective.  The  effects  of  ERa/p  heterodimer- 
selective  ligands  in  PC  3  cells  suggest  that  1  pM  cosmosiin 
and  10  pM  angolensin  are  responsible  for  mediating  the 
physiological  responses  of  ERa/p  heterodimers  on  a  cellular 
level  since  loss  of  either  ERa  or  ERp  abrogates  growth 
inhibition  at  these  concentrations  (Fig.  6C— F).  10  pM  cosmo¬ 
siin  mediates  growth  inhibitory  effects  via  ERp/p  homodimer¬ 
ization  and  off-target  effects  when  ERa  is  lost,  and  both 
concentrations  of  cosmosiin  increase  growth  via  ERa  homo¬ 
dimers  when  ERP  is  lost.  Therefore,  the  expression  levels  of 
ERs  appear  to  be  important  to  the  physiological  outcome  of 
these  ligands  at  cellular  levels. 


The  growth  effects  of  cosmosiin  and  angolensin  on 
additional  cell  lines  with  differing  ERocERp  expression 
ratios 

The  differing  cellular  effects  in  PC3,  PC3-shERa,  and  PC3- 
shERP  suggest  that  die  ratio  of  ERa:ERP  may  be  a  determinant  for 
the  ability  of  these  dimer-selective  ligands  to  act  in  a  proliferative  or 
anti-proliferative  manner.  To  address  this,  growth  and  viability 
assays  in  several  cell  lines  with  differing  expression  levels  of  ERa  and 
ERP  were  compared.  HC 1 1  is  a  normal  mouse  mammary  cell  line 
that  expresses  both  ERa  and  ERp  (Fig.  4A  and  [14]).  As  shown  in 
Fig.  S3,  cosmosiin  and  angolensin  are  both  able  to  inhibit  the 
growth  of  this  cell  line.  Specifically,  1  pM  angolensin,  a  concentra¬ 
tion  at  which  ERa/ P  heterodimers  are  not  predicted  to  be  activated 
(Figures  3B,  D)  has  no  effect  on  the  growth  of  this  cell  line,  whereas 
1 0  pM  angolensin  inhibits  HG 1 1  ’s  growth  by  ~  10%  compared  to 
the  vehicle  DMSO  (Fig.  S3D),  and  this  inhibition  is  ablated  by  the 
antagonist  ICI  182,780,  which  suggests  that  this  inhibition  is  ER- 
specific.  Cosmosiin  is  also  able  to  inhibit  the  growth  of  HC  1 1  cells 
at  1  pM  and  10  pM  (Fig.  S3C).  The  ~15%  inhibition  of  growth 
resulting  from  1  pM  cosmosiin  treatment  is  ablated  by  the 
antagonist  ICI  182,780,  indicating  that  this  response  is  ER-specific. 
10  pM  cosmosiin  inhibits  the  growth  of  HC11  cells  by  ~25% 
compared  to  the  vehicle  DMSO,  and  this  response  is  not  completely 
ablated  by  the  antagonist  ICI  1 82,780,  indicating  that  the  inhibition 
of  proliferation  by  1 0  pM  cosmosiin  is  not  ER-specific  in  agreement 
with  earlier  findings  (Fig.  6A  and  6D).  Cell  counting  and  viability 
assays  with  Trypan  blue  staining  confirmed  these  findings  of  growth 
inhibition  and  indicated  that  they  were  not  due  to  general 
cytotoxicity  (Fig.  S3A  and  S3B).  The  growth  inhibitory  effects  of 
1  pM  cosmosiin  and  1 0  pM  angolensin  in  HC  1 1  cells  support  the 
notion  that  ERa/p  heterodimers  are  growth  inhibitory.  Further¬ 
more,  we  examined  the  compounds’  effects  on  ERa— /ERp—  cell 
line  MDA-MB-231  and  ERa— /ERP+  DU- 145.  Neither  compound 
has  any  effect  on  cell  growth  at  all  tested  concentrations  in  MDA- 
MB-231  breast  cancer  cells  (Fig.  S4)  nor  DU- 145  prostate  cancer 
cells  (Fig.  S5).  This  result  suggests  that  the  growth  effects  exerted  by 
compounds  are  ERa  and  ERP-dependent  in  breast  and  prostate 
epithelial  cells.  This  conclusion  is  supported  by  the  findings  that 
growth  effects  elicited  by  1  pM  cosmosiin  and  10  pM  angolensin 
could  be  completely  antagonized  by  ER  antagonist  in  PC3,  PC3- 
shERa,  and  PC3-shERP  cells  (Fig.  6). 

Nuclear  co-localization  of  ERa  and  ERp  in  human  breast 
tumor  specimen 

Our  studies  indicate  that  cosmosiin  and  angolensin  could  be 
therapeutically  useful  for  inhibiting  the  growth  of  breast  cancer 
cells  that  co-express  ERa  and  ERp.  Although  previous  studies 
have  shown  60%  of  ERa-positive  breast  tumors  express  ERP 
[11,21],  in  order  for  ERa/p  heterodimerization  to  occur,  ERa 
and  ERP  must  be  co-expressed  in  the  same  cell.  To  investigate  the 
co-expression  of  ERa  and  ERP  in  breast  tumor  samples,  we 
analyzed  a  breast  cancer  tissue  microarray  (TMA)  using  the 
quantitative  immunofluorescence  AQUA®  technology  (HistoRx) 
that  allows  the  quantitative  measurement  of  proteins  of  interest 
within  subcellular  location  of  tissue  samples  by  calculation  of  an 
AQUA®  score.  Such  precision  is  not  possible  with  conventional 
testing  methods,  such  as  standard  immunohistochemistry  (IHC). 
This  TMA  was  purchased  from  US  Biomax  (BR2082)  and 
contained  32  cases  of  metastatic  carcinoma,  68  cases  of  invasive 
ductal  carcinoma,  22  cases  of  invasive  lobular  carcinoma,  22  cases 
of  intraductal  carcinoma,  4  cases  each  of  squamous  cell  carcinoma 
and  lobular  carcinoma  in  situ,  8  cases  of  fibroadenoma,  16  cases 
each  of  hyperplasia  and  inflammation,  10  cases  of  cancer  adjacent 
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Figure  5.  Cosmosiin  and  angolensin  inhibit  PC3  cell  motility  and  migration.  (A)  Wound  healing  assays  showing  the  effect  of  1  pM 
cosmosiin  and  10  pM  angolensin  on  ERoc,p-positive  PC3  cells.  Vehicle  (DMSO)  treatment  resulted  in  cell  motility  to  fill  the  wound  (top  panels)  that  was 
inhibited  by  1  pM  cosmosiin  ( middle  panels)  and  10  pM  angolensin  ( bottom  panels).  (B)  Transwell  migration  assays  measured  the  ability  of  cosmosiin 
(B)  and  angolensin  (C)  to  inhibit  cellular  migration  of  PC3  cells  toward  a  chemoattractant.  Cosmosiin  (B)  and  angolensin  (C)  decreased  the  ability  of 
PC3  cells  to  migrate  through  the  pore,  and  this  decreased  ability  was  ablated  by  the  antagonist  ICI  182,780  at  100  nM. 
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normal  breast  tissue  (NAT)  and  6  cases  of  normal  tissue.  A  total  of 
208  cores  were  analyzed  for  nuclear  ERa  and  ERP  intensity  with 
DAPI  and  P-actin  staining  as  references.  As  shown  in  Fig.  7A,  the 
ERa/ERP  ratio  increases  throughout  the  stages  of  carcinogenesis 
and  progression.  Pairwise  analysis  with  two-sample  t-tests  of  benign 
tissue  versus  hyperplasia  (p-value  =  0.0039)  and  versus  carcinoma 
(in  situ,  inflammation,  metastatic,  and  malignant  cases  with  p-values 
0.0092,  0.0035,  0.0042,  respectively)  indicate  that  this  ratio  is 
significantly  higher  in  cases  of  hyperplasia  and  carcinoma  compared 
to  benign  tissue.  Figure  7B  shows  that  ERa  and  ERP  colocalize 
within  the  nucleus  of  the  same  cell  in  tissue  samples.  Overlaying  the 
high  resolution  images  (Fig.  7B,  right)  for  ERa  staining  with  those  for 
ERP  staining  shows  that  ERa  and  ERP  co-localize  to  the  same  spots 
within  the  same  nucleus,  rendering  the  possibility  that  ERa/p 
heterodimerization  is  feasible  in  these  tissues. 

.  PLoS  ONE  |  www.plosone.org 


Discussion 

While  the  roles  of  ERa  and  ERP  in  hormone-dependent  diseases 
such  as  breast  and  prostate  cancers  are  becoming  increasingly 
elucidated,  with  ERa  having  a  proliferative  and  ERP  having  an  anti¬ 
proliferative  role,  the  mechanism  by  which  these  two  receptors 
interact  with  each  other  in  both  normal  and  diseased  states  has 
remained  elusive.  Because  the  co-expression  of  ERP  along  with  ERa 
dampens  the  proliferative  action  of  ERa,  direct  interaction  of  ERa 
and  ERP  is  thought  to  convey  growth  inhibitory  effects,  and  the 
ERa/p  heterodimer  has  been  proposed  to  activate  target  genes 
mediating  these  anti-proliferative  effects  [38,39].  However,  the 
heterogeneous  population  of  dimer  pairs  present  when  ERa  and 
ERP  are  co-expressed  and  the  lack  of  full  length  heterodimerized  ER 
structures  prevent  a  clear  understanding  of  their  biological  function. 
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Figure  6.  ER  dimer-selective  compounds  influence  cell  growth  in  an  ER-dependent  and  dose-dependent  manner.  Cosmosiin  (A)  and 
angolensin  (B)  decrease  the  growth  of  PC3  prostate  cancer  cells  in  a  dose-dependent  manner.  These  decreases  are  ER-specific  for  1  pM  cosmosiin  and 
1 0  pM  angolensin,  since  the  growth  decreases  are  ablated  by  the  ER  antagonist  ICI  1 82,780.  These  ER-specific  effects  by  1  pM  cosmosiin  (C)  and  1 0  pM 
angolensin  (D)  are  lost  with  the  silencing  of  ERa  in  PC3-shERoc  cells,  while  ER  non-specific  effects  due  to  10  pM  cosmosiin  are  retained.  Silencing  ER(3  in 
PC3-shER|3  results  in  cosmosiin-dependent  increases  in  cell  growth  (E)  that  are  ablated  in  the  presence  of  the  antagonist,  and  furthermore,  when  ERs  are 
antagonized,  ER  non-specific  growth  inhibition  in  PC3-shER|3  is  retained.  Angolensin  has  no  statistically  significant  effects  in  PC3-shER(3  (F).  Error  bars 
represent  standard  deviations  from  the  mean  of  triplicate  samples.  *  indicates  statistical  significance  by  the  Student's  T-Test. 
doi:1 0.1 371  /journal.pone.0030993.g006 


Thus,  in  order  to  shed  light  upon  the  biological  action  of  these  ERa/ 
P  heterodimers,  we  sought  to  identify  small  molecule  ligands  capable 
of  specifically  inducing  heterodimers  while  not  inducing  ERa/a 
homodimers  or  ERp/p  homodimers  with  the  rationale  that  diese 
ligands  could  be  used  to  decipher  the  biological  action  of  ERa/p 
heterodimers. 

The  BRET  technology  developed  in  our  lab  [34,35]  allowed  the 
examination  of  each  ER  dimer  pair  (ERa/ a  homodimers,  ERp/ P 
homodimers,  and  ERa/p  heterodimers)  in  isolation.  This 
segregation  was  especially  essential  in  the  case  of  the  ERa/p 
heterodimer,  as  the  co-expression  of  ERa  and  ERP  leads  to  the 
formation  of  all  three  dimer  forms  and  prevents  separation  of  the 
action  of  each  individual  dimer  pair  as  they  function  in  concert  in 


vivo.  However,  the  BRET  assay  allows  the  examination  of  ERa/  P 
without  observing  homodimer  formation.  Specifically,  we  have 
previously  shown  that  two  phytoestrogens,  genistein  and  liquir- 
itigenin,  preferentially  induce  different  ER  dimers  [34],  Liquir- 
itigenin  selectively  induced  formation  of  ERp/ P  homodimers  and 
ERa/p  heterodimers  but  not  ERa/a  homodimers  at  1  pM  [34], 
which  provides  proof-of-principle  that  small  molecule  compounds 
which  preferentially  induce  ERa/p  heterodimers  over  ERa/a 
homodimers  do  indeed  exist.  We  had  further  shown  that  BRET 
assays  can  be  optimized  for  HTS  [35] .  The  goal  of  secondary  HTS 
BRET  screening  in  this  study  was  to  find  a  compound  with  similar 
characteristics  to  liquiritigenin  but  with  greater  ERa/p  heterodi- 
mer  selectivity.  If  a  library  compound  was  able  to  induce  ERa/  P 
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Figure  7.  Automated  quantitative  measurement  of  ERa  and  ERp  expression  in  breast  cancer  tissue  microarrays.  (A)  Tissue  microarray 
analysis  with  the  AQUA®  technology  shows  that  the  ERocER(3  ratio  increases  from  benign  samples  throughout  various  stages  of  malignancy.  Pairwise 
two-sample  t-tests  between  the  benign  and  malignant  samples  showed  a  statistically  significant  difference  (p-values<0,01).  (B)  AQUA®  analysis 
indicates  that  ERa  and  ERp  colocalize  to  the  nucleus  within  the  same  cell  in  human  breast  tumors.  Figures  are  shown  at  200  x  magnification,  and 
scale  bars  are  present  in  the  lower  left  corners  of  each  image.  The  blow-up  picture  is  an  amplified  section  of  AQUA  staining  showing  co-localization  of 
ERa  and  ERp  to  the  nucleus.  *  indicates  statistical  significance  by  the  Student's  T-Test.  NAT  =  cancer  adjacent  normal  breast  tissue. 
doi:10.1371/journal.pone.0030993.g007 


heterodimerization  while  not  inducing  ER  homodimerization,  the 
ligand  could  be  used  in  biological  systems  to  determine  the 
function  of  these  heterodimers  with  minimal  interference  from 
either  homodimer. 

Two  lead  compounds  were  successfully  identified  in  BRET 
screening.  The  two  lead  compounds  are  flavonoids,  a  group  of 
potentially  chemoprotective  compounds  widely  distributed  in  fruit, 
vegetables,  and  beverages  of  plant  origin  including  tea  and  wine,  and 
have  similar  structures  that  consist  of  two  phenolic  benzene  rings 
linked  to  a  heterocyclic  pyre  or  pyrone  [40] .  Isoflavones  represent  an 
important  group  of  phytoestrogens  and  are  found  mainly  in  plants 
belonging  to  the  Leguminosae  family.  Angolensin  ( Trifolium  pretense, 
2l,4'-dihydroxy-4"-methoxy-a-methyldeoxybenzoin,  l-(2,4-dihydrox- 
yphenyl)-2-(4-methoxyphenyl)propan-l-one;  CAS  642-39-7),  is  an 
isoflavone  that  was  first  isolated  from  the  wood  of  Pterocarpus  angolensis 
and  later  from  the  wood  and  bark  of  Pterocarpus  indicus.  Angolensin  is  a 
metabolite  of  Biochanin  A  and  formononetin,  which  are  present  in  red 
clover  [4 1 ,42] .  Dietary  supplements  manufactured  from  red  clover  are 
widely  marketed  to  provide  beneficial  health  effects  of  isoflavones 
without  dietary  changes.  Specifically,  red  clover  supplements  are  often 
consumed  for  the  purported  alleviation  of  post-menopausal  symp¬ 
toms.  Cosmosiin  (apigenin  7-O-beta-glucoside;  apigenin-7-D-gluco- 
side;  apigenin-7-O-beta-D-glucopyranoside;  apigenin-7-glucoside; 
cosmetin,  Cosmosiine,  Apigetrin,5-hydroxy-2-(4-hydroxyphenyl)-7- 
[(2S,3R,4S,5S,6R)-3,4,5-trihydroxy-6-(hydroxymethyl)oxan-2-yl]oxy- 
chromen-4-one;  CAS  578-74-5)  is  a  flavonoid  present  in  chamomile 
flowers  which  are  used  pharmaceutically  and  cosmetically  for  their 
anti-spasmodic,  anti-inflammatory  and  antimicrobial  properties  and 


also  as  a  natural  hair  dye  and  fragrance.  Cosmosiin  has  also  been 
isolated  from  Veratrum  grandiflorum  (white  hellebore)  and  Kummerouia 
striata  (Korean  clover).  Cosmosiin  has  been  shown  to  exhibit  anti¬ 
inflammatory  properties  [43]  and  has  been  shown  to  exhibit  HIV  anti¬ 
viral  properties  [44],  although  it  has  not  received  FDA  approval  for 
these  purposes.  The  direct  binding  of  angolensin  and  cosmosiin  to  the 
E2-binding  pocket  of  ERs  are  observed  (Fig.  2).  To  our  knowledge, 
this  is  the  first  demonstration  of  cosmosiin  as  an  estrogenic  compound. 
Furthermore,  we  validated  ERa/p-heterodimer  specificity  using 
BRET  and  reporter  assays  and  showed  that  1  pM  cosmosiin  and 
10  pM  angolensin  are  specific  to  ERa/  P-heterodimers  (Fig.  3). 

Using  ERa/p-heterodimer  selective  compounds  at  specific 
concentrations,  we  are  able  to  show  that  the  ERa/ P-heterodimer 
is  growth  inhibitory.  These  compounds  inhibit  cell  proliferation  in 
HC1 1  and  PC3  cells  which  co-express  ERa  and  ERp.  Inhibition 
of  cell  growth  (Fig.  6)  and  migration  (Fig.  5)  due  to  1  pM 
cosmosiin  and  10  pM  angolensin  is  ablated  with  treatment  of  ICI 
182,780  or  the  silencing  of  either  ERa  or  ERP  in  PC3-shERa  and 
PC3-shERP,  respectively.  These  compounds,  however,  did  not 
have  an  effect  on  ER-negative  MDA-MB-231  and  ERa-negative/ 
ERP-positive  DU- 145  cells,  further  supporting  that  the  growth 
inhibitory  effects  observed  with  these  compounds  were  dependent 
on  expression  of  both  ERa  and  ERp.  While  these  compounds 
appear  to  have  little  or  no  effect  on  ERa/ a  homodimerization  and 
transcriptional  activation  in  HEK293  BRET  and  ERE-luciferase 
assays  employing  exogenous  ERs  (Fig.  3),  treatment  of  breast  and 
prostate  cancer  cells  expressing  ERa  at  a  much  higher  level  than 
ERp  (PC3-shERP)  results  in  ERa-dependent  growth  increases 
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(Fig.  6E).  This  result  is  in  agreement  with  the  common  theme  that 
ERa  is  a  major  growth  driver,  and  it  also  implicates  the 
dependence  of  these  compounds’  growth  effects  on  the  relative 
expression  ratio  of  ERa:ER(3,  as  these  compounds  ablate  growth 
increases  in  PC3  and  HC1 1,  in  which  ERa  and  ERP  expression 
levels  are  relatively  similar  and  heterodimerization  may  be  favored 
[31].  Taken  together,  these  data  suggest  that  the  ratio  ofERa:ERP 
in  the  same  tumor  cells  is  extremely  important  for  physiological 
effects  of  these  compounds.  While  the  data  presented  herein 
provide  initial  evidence  for  a  growth-inhibitory  function  of  the 
ERa/ P  heterodimer,  identification  of  higher  affinity  compounds 
with  greater  ERa/p  heterodimer  selectivity  will  be  needed  to 
validate  our  findings  since  both  compounds  are  weak  agonists,  and 
cosmosiin  at  10  pM  appears  to  have  off-target  effects. 

Compounds  exhibiting  ERa/p  heterodimer-selectivity  may 
have  therapeutic  or  preventive  efficacy  in  hormone-dependent 
diseases.  A  recent  study  shows  that  the  tamoxifen  metabolite 
endoxifen  is  capable  of  degrading  ERa  [45] ,  stabilizing  ERP,  and 
inducing  ERa/  P  heterodimerization  in  a  concentration  dependent 
manner  [46].  Tamoxifen  is  a  widely-utilized  FDA- approved  breast 
cancer  treatment  and  prevention  drug.  This  finding  suggests  that 
tamoxifen’s  cancer  preventive  effects  may  be  mediated  by 
stimulation  of  ERa/p  heterodimer  formation.  The  possibility  is 
supported  by  the  fact  that  both  ERs  are  expressed  in  normal 
mammary  epithelial  cells  [4],  Similarly,  naturally-occurring 
estrogen-like  compounds  such  as  phytoestrogens,  a  group  of 
plant-derived  compounds  with  estrogenic  and/ or  antiestrogenic 
activities  hold  promise  for  action  as  preventive  or  therapeutic  ER- 
regulators  via  their  abilities  to  mediate  estrogenic  responses  tissue- 
specifically.  Indeed,  consumption  of  soy  phytoestrogens  has  been 
correlated  with  decreased  breast  cancer  risk  [47],  although  these 
data  remain  somewhat  controversial  [48].  Furthermore,  consump¬ 
tion  of  genistein  [49],  resveratrol  [50],  and  soy  [51]  has  been 
inversely  correlated  with  prostate  cancer  risk.  Although  these 
compounds  may  stimulate  the  proliferative  action  of  ERa  when 
ERP  is  lost  in  tumors,  they  may  have  preventative  effects  under 
normal  physiological  conditions  when  both  ERs  are  expressed. 

Furthermore,  our  examination  of  nuclear  co-localization  of 
ERa  and  ERP  within  the  same  tumor  cell  using  the  AQUA® 
technology  (Fig.  7)  support  that  ERa/p  heterodimerization  could 
potentially  occur  within  tumor  cells.  Prior  to  these  studies,  the  co¬ 
localization  of  ERa  and  ERP  within  the  same  cell  had  not  been 
examined.  The  punctate  staining  pattern  suggests  that  ERa  and 
ERP  are  co-localized  on  DNA,  and  therefore  may  be  transcrip¬ 
tionally  active  in  these  cells  as  ERa/p  heterodimers.  Furthermore, 
AQUA®  analysis  showed  that  the  ERa:ERp  ratio  is  higher  in 
malignant  states  compared  to  benign  tissue  samples,  in  agreement 
with  the  finding  that  ERP  levels  often  decrease  in  malignant  breast 
cancers  [52].  The  growth  inhibitory  effects  of  ERa/p  heterodi¬ 
mers  might  due  to  their  activation  of  different  target  genes  from 
their  respective  homodimers.  Recently,  global  ChIP-Seq  analyses 
of  ERa  and  ERP  target  genes  show  that  perfectly  or  imperfectly 
palindromic  EREs  are  preferential  binding  sites  for  ERa/p 
heterodimers  as  compared  to  ERa/ a  or  ERp/ P  homodimers 
which  are  more  flexible  in  DNA  recognition  [53].  This  is 
consistent  with  other  reports  that  ERa/p  heterodimers  might 
regulate  distinct  genes  [32,33].  The  ERa/p  heterodimer-selective 
ligands  identified  in  this  study  will  allow  identification  of 
heterodimer  target  genes  in  cells  co-expressing  ERa  and  ERP 
(e.g.  PC3).  While  our  findings  implicate  the  ERa/p  heterodimer  as 
a  putative  preventative  and  therapeutic  target  for  hormone- 
responsive  cancers,  this  example  highlights  the  imminent  need  to 
decipher  the  role  diese  heterodimers  in  breast  and  prostate 
cancers. 


In  conclusion,  these  data  provide  a  proof-of-principle  that  ERa/ 
P  heterodimer-selective  ligands  can  inhibit  cell  growth  and 
migration  in  ERa/ERp-positive  cells  such  as  PC3  and  HG1 1  when 
ERa  and  ERp  are  expressed  at  similar  levels.  We  also  found  that  the 
compounds’  growth  effects  depend  on  the  relative  expression  levels 
of  ERa  and  ERp.  Upon  knockdown  of  ERp  in  PC3  cells,  cosmosiin 
increases  PC  3  cell  growth  in  an  ERa-dependent  manner.  Thus, 
more  heterodimer  selective  ligands  need  to  be  identified  to  clarify 
whether  the  heterodimer-selective  ligands  become  growth  stimula¬ 
tory  when  ERP  expression  is  lost  in  human  tumors.  Although  more 
studies  are  needed  to  demonstrate  the  ERa/p  heterodimer  as  a 
therapeutic  target,  the  concept  of  inducing  ERP  to  pair  with  ERa, 
thus  antagonizing  ERa’s  proliferative  function,  is  distinct  from 
existing  breast  cancer  therapeutic  strategies  of  targeting  ERa  alone. 
We  also  suggest  that  the  relative  ERa  and  ERP  expression  levels  in 
patient  tumors  should  be  carefully  evaluated  to  better  understand 
the  ER-targeted  drugs’  therapeutic  performance,  as  many  of  these 
drugs  have  not  been  evaluated  for  their  dimer  selectivity,  and  ERp 
expression  in  patient  tumors  is  not  routinely  evaluated. 

Materials  and  Methods 

High  Throughput  Screening  Methods 

All  primary  and  secondary  screens  were  performed  at  the 
University  of  Wisconsin  Carbone  Cancer  Center  (UWCCC)  Keck 
Small  Molecule  Screening  Facility  (SMSF).  Ten  thousand  T47D- 
KBLuc  cells  [54]  were  seeded  into  384-well  plates  and  allowed  to 
attach  overnight.  The  next  day,  0.5  pi  of  1  mM  compound  was 
added  to  a  final  concentration  of  10  pM  using  an  automated 
robotic  system  (Beckman  Biomek  FX).  10  nM  E2  and  1%  (0.5  pi) 
DMSO  were  used  as  positive  and  negative  controls,  respectively. 
Cells  were  incubated  with  compound  for  18  hrs  at  37°C  in  5% 
C02  in  a  cell  culture  incubator.  On  day  3,  media  were  removed, 
and  25  pi  lysis  buffer  (Promega,  cat#  E2661)  was  added  to  each 
well  using  the  robot.  Cells  were  allowed  to  lyse  for  10  min  with 
constant  agitation,  and  lysis  was  confirmed  by  microscopically 
viewing  a  clear-bottom  384-well  plate  maintained  in  parallel  under 
identical  conditions.  25  pi  luciferase  substrate  (Promega,  Cat# 
E2620)  was  then  added,  mixed  for  30  seconds,  and  luciferase 
emission  was  immediately  detected  on  a  Tecan  Safire  2  plate 
reader  at  0. 1  seconds  per  well.  Counter-screening  was  performed 
in  a  similar  fashion  in  the  presence  and  absence  of  the  ER 
antagonist  ICI  182,780.  Secondary  Bioluminescence  Resonance 
Energy  Transfer  (BRET)  screening  was  performed  in  transiently 
transfected  HEK293  cells  (ATCC,  CRL-1573).  DNA  encoding 
BRET  fusions  were  transfected  as  described  in  [34],  Following 
24  hours  of  protein  expression  after  transfection,  cells  were 
trypsinized  and  counted  using  a  Nexcelcom  Cellometer,  and  cell 
viability  was  determined  to  be  >95%  in  each  condition.  Cells  were 
seeded  at  11,000  cells  per  well  of  384-well  white-walled  white- 
bottom  plates  in  40  pL  PBS.  0.2  pL  of  1  mM  library  compound 
was  then  added  to  each  well  using  the  Biomek  FX  Robot  such  that 
the  final  concentration  per  well  was  5  pM.  Cell  suspensions  were 
incubated  with  library  compounds  for  1  hour  in  a  dark  cabinet  at 
room  temperature,  at  which  point  1 0  pL  of  the  Renilla  Luciferase 
(RLuc)  substrate  coelenterazine  h  was  added  to  a  final  concentra¬ 
tion  of  5  pM.  Plates  were  then  gently  shaken  on  a  plate  shaker  for 
10  seconds  at  300  rpm,  and  RLuc  emission  was  read  at  460  nm 
followed  immediately  by  YFP  emission  at  535  nm  at  0.1  second 
per  wavelength  read  per  well.  Each  RLuc  and  YFP  emission 
measurement  was  taken  consecutively  per  well  before  moving  to 
the  next  well.  Emission  values  were  used  to  calculate  the  BRET 
ratio  as  described  in  [34],  Additional  details  for  BRET  screening 
were  described  in  [35], 
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In  vivo  BRET  assays  to  monitor  ER  dimer  formation  in 
living  cells 

HEK293  cells  (ATCC,  CRL-1573)  were  either  transfected  with 
a  single  BRET  fusion  plasmid  (pCMX-ERa-RLuc  or  pCMX- 
RLuc-ERP)  or  co-transfected  with  RLuc  and  YEP  BRET  fusions 
(pCMX-ERa-RLuc+pCMX-YFP-ERP  for  ERa/ERp  heterodi¬ 
mers;  pCMX-ERa-RLuc+pCMX-ERa-YFP  for  ERa  homodimers; 
or  pCMX-RLuc-ERP+pCMX-YFP-ERP  for  ERp  homodimers) 
[34],  “Empty”  expression  vector  pCMX-pL2  was  used  to  keep  the 
total  amount  of  transfected  DNA  constant.  24  hr  post-transfection, 
cells  were  trypsinized,  counted,  and  resuspended  in  PBS  in 
quadruplicate  at  ~50,000  cells  per  well  of  a  96-well  white-bottom 
microplate.  Cells  were  incubated  with  ligands  for  1  hour.  Coelenter- 
azine  h  (Promega,  Madison,  WI)  was  added  in  PBS  at  a  final 
concentration  of  5  pM,  and  460  nm  and  530  nm  emission  detection 
measurements  were  immediately  taken  at  0. 1  second  per  wavelength 
read  per  well  on  a  Perkin  Elmer  Victor  3-V  plate  reader. 

Immunofluorescence  Staining 

Deparaffinization  and  heat  induced  epitope  retrieval  were 
performed  simultaneously  using  the  Lab  Vision  PT  module 
(Thermo  Fisher  Scientific,  Fremont,  CA)  with  Lab  Vision  citrate 
buffer  pFI  8.0  at  98°C  for  20  minutes.  All  staining  was  performed 
at  room  temperature  using  the  Lab  Vision  360  automated  staining 
system.  Endogenous  peroxidase  was  blocked  for  5  minutes  with 
Peroxidazed- 1  (Cat.No.  PX968,  Biocare  Medical).  Non-specific 
protein  binding  was  eliminated  via  a  60  minute  block  with  Biocare 
Medical  Sniper,  and  non-specific  avidin  was  blocked  using  Biocare 
Medical  Avidin  Biotin  kit,  incubating  15  minutes.  DaVinci  Green 
Antibody  Diluent  (Cat.No.  PD900L,  Biocare  Medical)  was  used 
for  antibody  dilution.  Breast  TMA  BR2082  containing  208  cores 
was  purchased  from  US  Biomax  Inc.  (http://www.biomax.us/ 
tissue-arrays/Breast/BR2082).  ERa  was  detected  using  ERa 
rabbit  mAb  SP1  (1:50,  1  hr)  (Thermo  Fisher)  and  visualized  with 
goat  anti-rabbit  conjugated  with  Alexa  Fluor  555  (Invitrogen) 
secondary  antibody.  ER|3  was  detected  with  mouse  mAb  14C8 
(Abeam,  1:1600,  1  hr)  and  visualized  with  Alexa  Fluor  647 
conjugated  Tyramide  Signal  Amplification  system  (Invitrogen), 
which  included  biotinylated  goat  anti-mouse  immunoglobulin, 
streptavidin-horseradish  peroxidase  and  Alexa  Fluor  647-Tyra- 
mide.  Breast  epithelial  nuclei  were  masked  using  ProLong  Gold 
Antifade  Reagent  with  DAPI  mounting  medium  (Invitrogen). 

Automated  Image  Acquisition 

Automated  image  capture  was  performed  by  the  HistoRx  PM- 
2000  using  the  AQUAsition  software  package  (New  Haven,  CT). 
High-resolution  (2048_2048  pixel,  7.4  mm),  8-bit  grayscale  digital 
images  are  obtained  for  each  area  of  interest  resulting  in  256 
discrete  intensity  values  per  pixel  of  an  acquired  image  [55].  The 
breast  epithelial  nuclear  compartment  was  defined  with  DAPI 
(blue).  The  target  markers  (ERa  and  ERP)  were  visualized  with 
Alexa  Fluor  555  (green)  and  647  (red),  respectively. 

AQUA  ®  Score  Generation 

Since  the  distributions  of  the  original  AQUA®  scores  exhibited 
deviation  from  the  normal  distribution,  we  took  the  natural  log 
transformation  of  the  original  scores  and  then  performed  two 
sample  t-tests  for  pairwise  comparisons  among  different  samples. 
Results  from  these  tests  were  consistent  with  a  Wilcoxon  rank  sum 
test  on  the  original  scores.  Images  were  evaluated  before  scoring. 
Histospots  showing  <5%  tumor  area,  tissue  folding,  too  much 
debris,  and  those  that  were  out  of  focus  were  disqualified  from 
scoring.  Nuclear  AQUA®  scores  for  ERa  and  ERP  for  each 


histospot  were  generated  based  on  the  unsupervised  pixel-cased 
clustering  algorithm  for  optimal  image  segmentation  for  use  in  the 
pixel-based  locale  assignment  for  compartmentalization  of  expres¬ 
sion  algorithm  as  described  previously  [56],  Pixels  that  could  not 
accurately  be  assigned  to  a  compartment  were  discarded.  The  data 
were  saved  and  subsequently  expressed  as  the  average  signal 
intensity  per  unit  of  compartment  area.  All  the  signals  in  each 
compartment  were  dien  added.  The  AQUA®  score  is  expressed  as 
target  signal  intensity  divided  by  the  compartment  pixel  area  and 
is  expressed  on  a  scale  of  0  to  33333  (AQUA_2.0,  HistoRx).  The 
resultant  AQUA®  score  is  continuous  and  directly  proportional  to 
the  number  of  molecules  per  unit  area. 

Additional  descriptions  of  cell  culture,  TMA  and  experimental 
procedures  can  be  found  in  Methods  SI. 

Supporting  Information 

Figure  SI  Caspase  3/7  apoptosis  assays  showed  that 
cosmosiin  and  angolensin  exhibited  no  apoptotic  effect 
via  caspases  3  or  7  in  PC3  cells  at  96  hours.  Cosmosiin  (A) 
and  angolensin  (B)  modesdy  increase  apoptosis  through  this 
pathway  in  PC  3  cells  to  a  statistically  non-significant  level  compared 
to  the  strong  apoptotic  inducer  cisplatin,  which  served  as  a  positive 
control.  Statistical  analysis  method:  Students  T-Tests.  Error  bars 
represent  standard  deviations  from  the  mean  of  triplicate  samples. 
(TIF) 

Figure  S2  MTT  assays  showing  the  effect  of  10  nM  17P- 
estradiol  in  ER«, [1-positive  PC3  cells  and  variants  of 
these  cells  in  which  ERa  has  been  silenced  (PC3-shERa) 
or  ERp  has  been  silenced  (PC3-shERP).  E2  has  no  effect  on 
the  proliferation  of  PC3  cells  (A)  or  PC3-shERa  cells  (B);  however, 
the  silencing  of  ER[I  in  this  cell  line  allows  E2  to  increase  cellular 
growth  (C,  left  panel)  by  binding  to  ERa,  since  the  presence  of  the 
antagonists  ICI  182,780  (C,  middle  panel)  and  MPP  Dihydro¬ 
chloride  (C,  right  panel)  ablates  this  increase.  Statistical  analysis 
method:  Students  T-Test;  *  indicates  p<0.05.  Error  bars  represent 
standard  deviations  from  the  mean  of  triplicate  samples. 

(TIF) 

Figure  S3  Cosmosiin  at  1  pM  and  angolensin  at  10  pM 
inhibited  the  growth  of  ERa/ERp  positive  HC11  cells  in 
an  ER-dependent  manner.  Cosmosiin  (A)  and  angolensin  (B) 
had  no  cytotoxic  effects  at  all  tested  concentrations.  The  growth 
inhibitory  effects  of  cosmosiin  at  1  pM  (C)  and  angolensin  at  1 0  pM 
(D)  were  ablated  by  pure  ER  antagonist  ICI  182,780,  suggesting  the 
growth  inhibitory  effects  are  ER-dependent.  These  decreases  due  to 
1 0  pM  cosmosiin  are  ER-independent  since  they  are  retained  in  the 
presence  of  the  antagonist  ICI  182,780  (C).  Statistical  analysis 
method:  Students  T-Test;  *  indicates  p<0.05.  Error  bars  represent 
standard  deviations  from  the  mean  of  triplicate  samples. 

(TIF) 

Figure  S4  Neither  cosmosiin  (A)  nor  angolensin  (B) 
influenced  the  growth  of  ER-negative  MDA-MB- 
231breast  cancer  cells.  Error  bars  represent  standard 
deviations  from  the  mean  of  triplicate  samples. 

(TIF) 

Figure  S5  Neither  cosmosiin  (A)  nor  angolensin  (B) 
influenced  the  growth  of  ERa-negative,  ERP-positive  DU- 
145  prostate  cancer  ceUs.  Error  bars  represent  standard 
deviations  from  the  mean  of  triplicate  samples. 

(TIF) 

Methods  SI  Supplemental  Methods  file. 

(DOC) 
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The  histone  coactivator-associated  arginine  methyltransferase  1 
(CAR Ml)  is  a  coactivator  for  a  number  of  transcription  factors, 
including  nuclear  receptors.  Although  CARM1  and  its  asymmetrically 
deposited  dimethylation  at  histone  H3  arginine  17  (H3R17me2a)  are 
associated  with  transcription  activation,  the  mechanism  by  which 
CARM1  activates  transcription  remains  unclear.  Using  an  unbiased 
biochemical  approach,  we  discovered  that  the  transcription  elon¬ 
gation-associated  PAF1  complex  (PAFIc)  directly  interacts  with 
H3R17me2a.  PAFIc  binds  to  histone  FH3  tails  harboring  dimethyla¬ 
tion  at  R17  in  CARM1 -methylated  histone  octamers.  Knockdown  of 
either  PAFIc  subunits  or  CARM1  affected  transcription  of  CARM1- 
regulated,  estrogen-responsive  genes.  Furthermore,  either  CARM1 
knockdown  or  CARM1  enzyme-deficient  mutant  knockin  resulted  in 
decreased  H3R17me2a  accompanied  by  the  reduction  of  PAFIc  oc¬ 
cupancy  at  the  proximal  promoter  estrogen-responsive  elements.  In 
contrast,  PAFIc  knockdown  elicited  no  effects  on  FI3R17me2a  but 
reduced  the  FI3K4me3  level  at  estrogen-responsive  elements.  These 
observations  suggest  that,  apart  from  PAFIc's  established  roles  in 
directing  histone  modifications,  PAFIc  acts  as  an  arginine  methyl 
histone  effectorthat  is  recruited  to  promoters  and  activates  a  subset 
of  genes,  including  targets  of  estrogen  signaling. 

estrogen  receptor-a  |  histone  arginine  methylation  |  PRMT 

Coactivator-associated  arginine  methyltransferase  1  (CARM1), 
also  known  as  PRMT4,  is  the  first  protein  arginine  methyl¬ 
transferase  identified  as  a  coactivator  for  steroid  receptors  (1).  It  was 
later  found  to  be  a  coactivator  for  a  variety  of  transcription  factors 
and  is  a  multifunctional  protein  engaged  in  a  variety  of  cellular 
processes,  including  gene  expression,  coupling  of  transcription  with 
mRNA  splicing,  regulation  of  protein  stability,  and  functioning  in 
embryonic  development  (2).  Loss  of  CARM1  in  the  mouse  embryo 
leads  to  abrogation  of  the  estrogen  response  and  reduced  expression 
of  some  estrogen  receptor-a  (ERa)  target  genes  (3).  Recently, 
CARM1  was  shown  to  be  a  unique  coactivator  of  ERa  that  can  si¬ 
multaneously  block  cell  proliferation  and  induce  differentiation 
through  global  regulation  of  ERa-regulated  genes  in  breast  cancer 
cells  (4).  The  methyltransferase  activity  of  CARM1  is  required  for 
its  ability  to  regulate  transcription  of  pS2  (TFF1),  a  standard  ERa 
target  gene  (5).  Moreover,  the  enzyme-dead  CARM1  knockin  mice 
have  defects  similar  to  those  seen  in  their  knockout  counterparts  (6). 
These  observations  suggest  that  the  enzymatic  activity  of  CARM1  is 
indispensable  for  the  majority  of  CARMl’s  in  vivo  functions.  Sev¬ 
eral  putative  mechanisms  have  been  proposed  to  explain  CARM1 
activator  function,  including  histone  H3  methylation  at  arginine  17 
(R17),  methylation  of  other  key  transcription  coactivators,  and  re¬ 
cruitment  of  chromatin  remodeling  factors.  However,  it  is  unclear 
which  mechanisms  drive  CARM1  transcription  activation  function. 

CARM1  was  originally  characterized  as  a  histone  arginine 
methyltransferase  because  of  its  ability  to  dimethylate  arginines 
on  histone  H3.  Subsequent  biochemical  analyses  showed  that 
histone  H3  can  be  modified  at  R2,  R17,  and  R26  (7).  A  recent 
study  showed  that  asymmetrically  dimethylated  histone  H3  at  R2 


(H3R2me2a)  was  mainly  deposited  by  PRMT6  (8).  ChIP  analysis 
identified  elevated  levels  of  H3R17me2a  at  the  estrogen-responsive 
pS2  promoter  (5).  Furthermore,  kinetic  ChIP  analysis  revealed  that 
CARM1  was  recruited  in  a  cyclic  manner,  occurring  at  40-min 
intervals  upon  treatment  with  17p-estradiol  (E2),  and  recruitment 
of  CARM1  correlated  with  an  increase  in  the  H3R17me2a  mark 
(9).  These  studies  suggest  that  H3R17me2a  likely  accounts  for 
CARMl’s  coactivator  function  in  transcriptional  regulation.  Be¬ 
cause  many  protein  domains  have  been  characterized  to  specifically 
bind  to  modified  histone  marks,  effector  molecules  may  exist  that 
read  H3R17me2a  marks  to  mediate  transcriptional  activation.  In 
support  of  this  hypothesis,  a  recent  study  identified  one  effector 
protein,  TDRD3,  as  a  “reader”  for  H3R17me2a  and  H4R3me2a 
using  a  protein  domain  microarray  approach  (10).  TDRD3  con¬ 
tains  a  tudor  domain  that  is  typical  for  mediating  methyl-specific 
binding.  Furthermore,  TDRD3  functions  as  a  coactivator  in  es¬ 
trogen-responsive  elements  (ERE)-luciferase  reporter  assays  and 
endogenous  TDRD3  was  detectable  at  the  pS2  promoter  in  ChIP 
assays  (10).  Despite  these  observations,  it  is  unclear  how  TDRD3 
promotes  transcriptional  activation. 

To  further  elucidate  the  mechanism  by  which  CARM1  enzymatic 
activity  and  H3R17me2a  mediate  transcription  activation,  we  took 
an  unbiased  biochemical  approach  to  identify  effector  proteins  that 
directly  interact  with  the  H3R17me2a  mark  using  peptide  pull¬ 
down  assays  from  HeLa  nuclear  extract.  We  unexpectedly  identified 
the  transcription  elongation-associated  PAF1  complex  (PAFIc)  in 
direct  association  with  H3R17me2a  and  found  this  association  is 
required  for  complete  activation  of  CARM1 -regulated  ERa  target 
genes.  Furthermore,  PAFIc  was  detected  at  EREs,  and  its  occu¬ 
pancy  was  abrogated  by  either  knockdown  of  CARM1  or  in¬ 
troduction  of  a  CARM1  enzyme-deficient  mutant,  which  was 
accompanied  by  the  attenuation  of  H3R17me2a.  This  study  pres¬ 
ents  a  unique  mechanism  by  which  H3R17me2a  functions  to  recruit 
PAFIc  to  facilitate  transcription  activation  and  encrypt  a  means  for 
PAFIc’s  recruitment  to  the  promoters  of  target  genes. 

Results 

Identification  of  Human  PAFIc  as  a  Binding  Partner  for  the  Histone  H3 
Tail  Containing  the  Asymmetrically  Dimethylated  Arginine  17.  To 

identify  proteins  that  recognize  asymmetrically  dimethylated  ar¬ 
ginine  17  at  histone  H3  (H3R17me2a),  we  took  an  unbiased  bio¬ 
chemical  approach  using  biotinylated  histone  tail  pull-down  from 
HeLa  nuclear  extract,  according  to  previously  published  methods 
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(11-13).  Two  histone  H3  N-terminal  tail  peptides  encompassing 
amino  acids  1-20,  with  or  without  asymmetrically  dimethylated  R17 
(H3R17me2a),  were  used  for  pull-down  assays.  Equal  amounts  of 
biotinylated  peptides  preimmobilized  on  streptavidin  magnetic 
beads  were  incubated  with  HeLa  nuclear  extract  followed  by  acid 
elution  (12).  Eluted  proteins  were  immediately  neutralized,  loaded 
on  SDS/PAGE,  and  silver-stained.  Fig.  L4  showed  that  three  bands 
were  specifically  associated  with  H3R17me2a  tail  compared  with 
the  nonmodified  H3  tail.  This  experiment  was  repeated  on  a  large 
scale,  and  the  bands  were  Coomassie-stained,  excised,  and  ana¬ 
lyzed  by  MALDI-TOF  MS.  Three  bands  were  identified  as  hCtr9, 
hCdc73,  and  hSki8,  all  of  which  belong  to  the  same  multisubunit 
PAFlc.  The  matched  peptides  sequences  are  shown  in  Fig.  SI.  The 
PAFlc  was  originally  identified  in  yeast  as  a  RNA  polymerase  II 
(RNAPII)-associated  complex  mediating  transcription  elongation 
(14).  The  human  PAFlc  (hPAFlc)  has  been  shown  to  be  involved 
in  multiple  steps  in  transcription,  including  control  of  F12B  ubiq- 
uitination  and  H3K4me3,  transcription  initiation  and  elongation, 
and  pre-mRNA  processing  (15).  Because  hPAFlc  is  composed  of 
five  subunits  (hCtr9,  hLeol,  hPafl,  hCdc73,  and  hSki8),  we  spec¬ 
ulated  that  the  whole  complex  (including  two  missing  subunits, 
hLeol  and  hPafl)  bound  to  the  H3R17me2a  tail.  Their  absence  in 
silver  staining  (Fig.  L4)  could  be  because  of  inefficient  elution  or 
poor  staining  by  silver.  The  antibodies  for  individual  subunits  were 
used  to  detect  their  presence  in  eluted  fractions  from  3  biotinylated 
peptide  pull-downs.  The  results  showed  that  although  all  PAFlc 
subunits  are  strongly  detected  in  H3R17me2a  tail  pull-down,  much 
weaker  PAFlc  binding  was  detected  with  the  nonmodified  F13  tail, 
and  no  binding  was  detected  with  negative  control  peptide  pull- 
downs  (Fig.  IB). 
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Fig.  1.  The  RNAPII-associated  PAFlc  displays  binding  specificity  to  asym¬ 
metrically  dimethylated  histone  H3R17me2a.  (A)  Immobilized  histone  H3  tail 
(amino  acids  1-20,  Left)  and  the  identical  peptide  carrying  asymmetric 
dimethylated  arginine  17  (Right)  were  used  to  identify  specific  effectors  for 
the  H3R17me2a  mark  in  HeLa  cell  nuclear  extract.  Three  proteins  were 
specifically  pulled  down  with  H3R17me2a  and  later  identified  to  be  hCtr9, 
hCdc73,  and  hSki8  by  in-gel  digestion  and  mass  spectrometry  analyses.  (6) 
Western  blots  showed  that  biotinylated  H3R17me2a  peptide  selectively  pull 
down  PAFlc  subunits,  including  Ctr9,  Leol,  Pafl,  Cdc73,  and  Ski8.  (C)  The 
PAFlc  was  affinity-purified  from  HEK293-Flag-Cdc73  cells  and  silver-stained. 
(D)  Coomassie  staining  of  the  reconstituted  histone  octamers.  (£)  PAFlc 
preferentially  interacts  with  CARM1  methylated  histone  octamers,  which 
were  reconstituted  from  bacterially  expressed  histone  H3,  H2A,  H2B,  and 
H4.  Reconstituted  histone  octamers  were  in  vitro  methylated  by  CARM1  or 
PRMT1  in  the  presence  of  3H-AdoMet  and  incubated  with  PAFlc  immobi¬ 
lized  on  anti-Flag  resin.  Equal  volume  of  input,  bound,  and  flow-through  of 
PRMT1 -methylated  3H-histone  H4R3me2  (lanes  1,  3,  and  5)  and  CARM1- 
methylated  3H-histone  H3R17me2  (lanes  2,  4,  and  6)  were  loaded.  Autora¬ 
diography  shows  that  more  CARM1 -methylated  3H-histone  H3R17me2  was 
found  bound  to  PAFlc  on  the  beads  compared  with  PRMT1 -methylated 
3H-histone  H4R3me2. 


To  identify  the  PAFlc  subunit  responsible  for  direct  in¬ 
teraction  with  the  histone  FI3R17me2a  tail,  individual  PAFlc 
subunits  were  recombinantly  expressed  as  Flag-tagged  proteins 
in  Sf9  insect  cells  and  purified  (Fig.  S2/4).  Purified  subunits  were 
individually  used  for  pull-down  assays.  Western  blots  showed 
that  Pafl  subunit  elicited  the  strongest  binding,  whereas  all  other 
subunits  displayed  either  weak  (Cdc73)  or  negligible  binding  to 
histone  H3  tails  (Fig.  S2B).  This  result  suggested  that  the  Pafl 
subunit  is  the  main  component  in  PAFlc  that  interacts  with  the 
histone  H3  tail.  The  observation  that  Pafl  unexpectedly  binds  to 
both  nonmodified  and  H3R17me2a  tails  indicates  that  other 
subunits  in  the  complex  may  contribute  to  the  specificity  of  Pafl. 

The  recently  identified  H3R17me2a  effector  molecule  TDRD3 
was  also  found  to  interact  with  F14R3me2a,  a  site  modified  by 
PRMT1  on  histone  F14  (10).  To  investigate  if  PAFlc  interacts 
with  the  PRMTl-modified  histone  mark,  we  reconstituted  his¬ 
tone  octamers  from  bacterially  expressed  core  histones  (Fig.  ID) 
and  in  vitro  methylated  histone  octamers  by  PRMT1  in  the 
presence  of  3H-AdoMet  for  binding  assays.  Similarly,  histone 
octamers  were  methylated  by  CARM1  as  a  positive  control.  The 
in  vitro  methylation  of  reconstituted  histone  octamers  by  PRMTs 
generated  3H-labeled  octamers  uniquely  carrying  either  asym¬ 
metric  dimethylation  of  R17  on  histone  FI3  (H3R17me2a)  me¬ 
diated  by  CARM1  or  asymmetric  dimethylation  of  R3  on  histone 
F14  (H4R3me2a)  mediated  by  PRMT1.  PAFlc  was  purified  from 
Flag-Cdc73  stably-expressing  F1EK293  cells,  as  shown  by  silver 
staining  (Fig.  1C),  and  preimmobilized  on  anti-Flag  M2  resin. 
Equal  amounts  of  preimmobilized  PAFlc  were  incubated  with 
JH-H3R17me2a  or  ’H-H4R3me2a  histone  octamers  in  parallel. 
The  input,  flow-through,  and  bead-bound  samples  were  resolved 
on  SDS-PAGE,  and  the  dried  gel  was  exposed  to  an  X-ray  film  to 
detect  3H-labeled  methylated  octamers  in  each  fraction.  The 
autoradiography  showed  that  the  efficiency  of  octamer  methyl¬ 
ation  by  CARM1  and  PRMT1  were  similar,  as  indicated  by  the 
input  samples  containing  equivalent  radioactivity  (Fig.  IE, 
compare  lane  1  with  lane  2).  Flowever,  PAFlc  strongly  inter¬ 
acted  with  the  3H-H3R17me2a  histone  octamers  methylated  by 
CARM1  (Fig.  1 E,  compare  lane  4  with  lane  3).  Consistently, 
more  PRMT1  -methylated  3H-F14R3me2a  histone  octamers  were 
detected  in  the  flow-through  (Fig.  IE,  compare  lane  5  with  lane 
6).  Our  results  demonstrated  that,  unlike  TDRD3,  PAFlc  pref¬ 
erentially  binds  to  CARMl-methylated  histone  octamers  over 
those  methylated  by  PRMT1.  Moreover,  PAFlc  not  only  inter¬ 
acts  with  histone  H3  tail  peptides  but  also  histone  octamers 
carrying  asymmetric  dimethylation  at  R17  of  histone  H3. 

Estrogen-Enhanced  Cooccupancy  of  PAFlc  and  H3R17me2a  at  EREs  of 
ERa  Target  Genes.  Given  that  H3R17me2a  peptide  selectively 
binds  to  PAFlc  in  vitro,  we  next  examined  if  PAF I  c  is  recruited 
to  ERa  target  genes  that  harbor  the  FI3R17me2a  mark  in  vivo. 
pS2,  PTGES,  and  IGFBP4  genes  were  selected  for  ChIP  analyses 
because  they  were  previously  shown  as  CARMl-regulated  ER 
target  genes  (1,  16)  and  contain  defined  EREs.  In  agreement 
with  past  studies,  H3R17me2a  was  enriched  at  EREs  of  selected 
genes  upon  17p-estradiol  (E2)  treatment  (Fig.  2  A,  E,  and  I). 
Notably,  Pafl  (Fig.  2  B,  F,  and  J)  and  Cdc73  (Fig.  2  C,  G,  and  K) 
occupancies  at  EREs  were  also  enhanced  by  E2  treatment.  To 
validate  the  specificity  of  Cdc73  antibody  used  in  ChIP  experi¬ 
ments,  a  MCF7-tet-on-shCdc73  cell  line  was  constructed  and 
used  for  ChIP  assays  using  anti-Cdc73  antibody  (Fig.  2 M).  E2- 
induced  recruitment  of  Cdc73  to  EREs  was  found  decreased 
after  knockdown  of  Cdc73  in  MCF7  cells  (Fig.  2N).  Taken  to¬ 
gether,  these  results  suggest  that  intact  PAFlc  is  recruited  to 
EREs  in  an  E2-enhanced  manner,  and  the  recruitment  of  PAFlc 
is  accompanied  by  the  increase  of  H3R17me2a  mark. 

Recruitment  of  hPAFlc  to  EREs  Is  CARM1 -Dependent.  To  determine 
if  PAFlc  is  recruited  to  EREs  via  recognition  of  H3R17me2a 
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Fig.  2.  E2  induces  cooccupancy  of  PAFIc  and  H3R17me2a  at  the  EREs  of 
several  ERa-target  genes:  pS2,  PTGES,  and  IGFBP4.  Cells  were  treated  with  or 
without  E2  for  45  min  followed  by  ChIP  and  quantitative  PCR  assays  using 
the  indicated  antibodies.  PAFIc  occupancy  at  the  ERE  regions  (B,  F,  and  J  for 
Pafl;  C,  G,  and  K  for  Cdc73)  of  ERa  target  genes  coincides  with  the  presence 
of  the  H3R17me2a  mark  (A,  E,  and  /).  ( D ,  H,  and  L)  IgG  control  for  pS2,  PTGES, 
and  IGFBP4  ERE  regions.  ( M )  MCF7-tet-on-shCdc73  cell  line  was  generated 
and  knockdown  efficiency  of  Cdc73  was  shown  in  Western  blots.  (A/)  Cdc73 
occupancy  at  pS2  ERE  region  was  decreased  after  knockdown  of  Cdc73  in 
MCF7-tet-on-shCdc73  cells.  Error  bars  represent  SD  (n  =  3). 
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Fig.  3.  The  recruitment  of  PAFIc  to  proximal  promoter  EREs  of  ERa  target 
genes  is  CARM1-  and  H3R17me2a-dependent.  MCF7-tet-on-shCARM1  cells 
were  pretreated  with  or  without  0.5  (ig/mL  Dox  for  4  d  before  E2  or  vehicle 
treatment  for  45  min.  The  occupancy  of  PAFIc  subunit  Cdc73  at  EREs  of  pS2 
(A)  and  PTGES  (B)  was  analyzed  using  ChIP  assays.  When  CARM1  is  knocked 
down  with  Dox  treatment,  E2  stimulated  FI3R17me2a  and  PAFIc  occupancy 
at  EREs  of  pS2  and  PTGES  were  greatly  reduced  (compare  lane  4  with  lane  2). 
(C)  Western  blots  showed  the  expression  levels  of  CARM1WT  and  CARM1VLD 
were  comparable  to  the  endogenous  CARMI  in  MCF7  cells  when  endoge¬ 
nous  CARM1  was  efficiently  silenced.  (D  and  E)  ChIP  analyses  showed  that 
expression  of  the  enzyme-defective  CARM1VLD  mutant  was  insufficient  to 
restore  FI3R17me2a  level  (D,  *P  <  0.01)  and  recruit  PAFIc  to  pS2  ERE  (£,  *P  < 
0.01).  In  contrast,  CARM1WT  can  fully  rescue  the  function  of  endogenous 
CARMI  (histograms  3  and  4).  Error  bars  represent  SD  (n  =  3). 


mark,  we  used  a  CARMI  inducible  knockdown  cell  line  (4)  for 
ChIP  assays  because  CARMI  is  the  only  reported  enzyme  to 
create  H3R17me2a  mark.  MCF7-tet-on-shCARMl  cells  were 
treated  with  DMSO,  E2,  doxycycline  (Dox),  or  Dox+E2  and  sub¬ 
jected  to  ChIP  assays.  E2  treatment  induced  H3R17me2a  accu¬ 
mulation  at  the  ERE  of  pS2,  whereas  E2-induced  H3R17me2a  level 
was  greatly  reduced  when  CARMI  was  knocked  down  by  Dox 
treatment  (Fig.  3 A,  compare  lane  2  with  lane  4).  A  similar  obser¬ 
vation  was  made  with  the  PTGES  gene  under  identical  treatment 
conditions  (Fig.  3 B).  These  results  confirm  that  CARMI  is  re¬ 
sponsible  for  writing  the  H3R17me2a  mark  in  response  to  E2.  Using 
this  cell-line  model,  we  found  that  Cdc73  E2-induced  occupancy  at 
the  EREs  of  pS2  or  PTGES  was  greatly  abrogated  after  C  ARM  I 
knockdown  (Fig.  3  A  and  B,  compare  lane  2  with  lane  4).  Because 
(  ARM  I  has  been  implicated  in  transcription  elongation  and 
splicing  processes,  complete  loss  of  CARMI  may  indirectly  affect 
transcription  regardless  of  histone  H3R17  dimethylation.  To  ad¬ 
dress  this  possibility,  we  replaced  wild-type  CARMI  with  an  en¬ 
zyme-deficient  CARMI  mutant  by  coinfecting  MCF7  cells  with 
lentiviruses  encoding  CARMlshRNA  and  FlagCARMI VLD,  a  well- 
characterized  CARMI  mutant  defective  for  methyltransferase  ac¬ 
tivity  (1).  As  shown  in  Fig.  3C,  endogenous  CARMI  was  successfully 
knocked  down  and  the  CARMULD  mutant  was  expressed.  As 
a  control,  we  also  restored  wild-type  CARMI  under  the  endogenous 
CARMI  knockdown  background.  ChIP  assays  were  performed  5 
d  after  viral  infection  and  E2  or  vehicle  treatment  of  cells  for  45  min. 
The  results  showed  that  E2-enhanced  H3R17me2a  level  (Fig.  3D) 
and  Cdc73  recruitment  at  ERE  of  pS2  (Fig.  3 E)  were  significantly 
impaired  in  CARMVLD  compared  with  CARMl'^-expressing  cells. 


These  data  suggest  that  PAFIc  occupancy  at  EREs  is  (  ARM  I -de¬ 
pendent,  most  likely  mediated  by  the  H3R17me2a  mark. 

PAFIc  and  CARMI  Coordinately  Regulate  a  Common  Set  of  ERa 
Target  Genes.  We  had  previously  identified  CARMl-regulated 
ERa  target  genes  using  microarrays  (4).  To  establish  the  func¬ 
tional  role  of  PAFIc  in  CARMl-regulated  transcription,  we 
attempted  to  knockdown  PAFIc  components  and  measure  the 
effects  on  mRNA  levels  of  CARMl-regulated  ER  targets. 
Knocking  down  Ctr9  was  previously  reported  to  affect  expression 
of  other  PAFIc  subunits  in  yeast  and  human  cells  (17,  18),  and 
Ctr9  and  Pafl  were  shown  as  scaffold  proteins  for  the  formation 
of  PAFIc  (19).  Therefore,  we  generated  a  Dox-inducible  Ctr9 
knockdown  cell  line,  MCF7-tet-on-shCtr9,  which  allows  transient 
knockdown  of  Ctr9  to  minimize  cellular  toxicity  associated  with 
the  loss  of  functional  RAF!  c.  To  examine  if  Ctr9  knockdown 
interferes  with  PAFIc  integrity,  MCF7-tet-on-shCtr9  cells  were 
infected  with  lentiviral  Flag-tagged  Ski8  to  immunoprecipitate 
PAFIc  with  anti-Flag  affinity  resin.  This  cell  line  was  treated  with 
Dox  or  vehicle  for  5  d  to  knock  down  Ctr9,  and  then  both  whole¬ 
cell  lysates  and  anti-Flag  immunoprecipitated  PAFIc  were  sub¬ 
jected  to  Western  blots.  In  accompaniment  with  Ctr9  knock¬ 
down,  the  endogenous  Pafl  level  was  significantly  reduced  in 
whole  lysate  (Fig.  S3,  lanes  1  and  2)  as  well  as  in  Flag-Ski8 
immunoprecipitated  PAFIc  (Fig.  S3,  lanes  3  and  4).  Inter¬ 
estingly,  despite  no  detectable  loss  of  endogenous  Leol  expres¬ 
sion,  the  presence  of  Leol  in  PAFIc  was  diminished,  whereas 
Cdc73  level  was  moderately  affected  (Fig.  S3).  Our  findings  are 
in  agreement  with  previous  reports  that  Ctr9  plays  a  pivotal  role 
in  PAFIc  assembly  (19),  and  further  indicate  that  Leol  might 
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directly  interact  with  Pafl  in  PAFlc.  Thus,  the  MCF7-tet-on- 
shCtr9  cell  line  is  validated  as  a  good  model  for  elucidating  the 
function  of  PAFlc  in  E2-dependent  gene  expression. 

Because  the  MCF7-tet-on-shCARMl  cell  line  (4)  was  estab¬ 
lished  from  the  same  parental  MCF7-tet-on  clone  as  MCF7-tet- 
on-shCtr9,  we  compared  expression  of  ERa  target  genes  in  both 
cell  lines  simultaneously.  The  Ctr9  and  CARM1  knockdown  ef¬ 
ficiencies  were  confirmed  by  Western  blots  showing  that  5  d  of  Dox 
treatment  was  sufficient  to  robustly  decrease  both  protein  levels  by 
90%  (Fig.  4  A  and  B ).  To  demonstrate  that  loss  of  CARM1  does 
not  affect  steady-state  levels  of  any  of  the  PAFlc  subunits,  we 
performed  Western  blots  of  PAFlc  subunits  in  MCF7-tet-on- 
shCARMl  cell  line.  The  results  showed  that  PAFlc  subunit  levels 
were  not  altered  by  the  loss  of  CARM1  (Fig.  S4).  The  Ctr9  and 
CARM1  inducible  knockdown  cells  were  treated  under  four 
conditions  and  harvested  for  mRNA,  followed  by  quantitative  RT- 
PCR  (qRT-PCR)  analysis.  pS2,  PTGES,  IGFBP4,  EGR3,  and 
MY C  genes  were  selected  for  analyses  because  they  were  found  to 
be  CARMl-dependent  E2  target  genes  (4).  Fig.  4  showed  that  E2 
induced  expression  of  all  these  genes  by  more  than  twofold,  and 
induction  was  inhibited  over  50%  by  either  Ctr9  or  CARM1 
knockdown  (Fig.  4  C-G).  Moreover,  the  trends  for  these  ERa 
target  genes  in  MCF7-tet-on-shCARMl  and  MCF7-tet-on-shCtr9 
cell  lines  under  the  four  treatment  conditions  were  similar  (Fig.  4 
C-G).  These  findings  were  further  confirmed  in  MCF7  cells 
transiently  transfected  with  siRNAs  targeting  Cdc73  or  Pafl  sub¬ 
units  (Fig.  S5).  We  also  examined  PAFlc  dependency  of  three  E2- 
induced,  CARM1  nonregulated  genes  (4)  by  qRT-PCR  in  two 
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Fig.  4.  Expression  of  several  selected  ERa  target  genes  were  similarly  af¬ 
fected  by  knockdown  of  CARM1  or  Ctr9.  MCF-tet-on-shCtr9  and  MCF7-tet- 
on-shCARMI  cell  lines  were  treated  with  vehicle,  E2  (5  h),  Dox  (5  d),  or  Dox 
(5  d)  +  E2  (5  h)  in  parallel.  Endogenous  Ctr9  (A)  and  CARM1  ( B )  were  knocked 
down  by  90%  after  5  d  of  Dox  treatment.  The  mRNA  levels  of  pS2  (C),  PTGES 
(D),  IGFBP4  (£),  MYC  (F),  and  EGR3  (G)  genes  were  analyzed  by  qRT-PCR  in 
both  cell  lines  after  the  indicated  treatment.  The  attenuation  of  estrogen 
response  was  observed  with  all  examined  ER  target  genes  in  both  Ctr9  and 
CARM1  knockdown  cells  after  Dox  induced  silencing  of  Ctr9  or  CARM1,  re¬ 
spectively.  Error  bars  represent  SD  (n  =  3). 


inducible  cell  lines.  Slc7a5  and  Cyclin  D1  were  found  not  to  be 
regulated  by  either  CARM1  or  PAFlc  (Fig.  S6  A  and  B),  whereas 
TGF-a  was  found  to  be  regulated  by  PAFlc  but  not  by  CARM1 
(Fig.  S6C).  These  data  imply  that  PAElc  shares  both  common  and 
distinct  sets  of  genes  from  those  of  CARM1.  Collectively,  our  data 
strongly  support  that  PAFlc  and  CARM1  coregulate  at  least 
a  subset  of  genes  in  the  estrogen  pathway. 

PAFlc  Is  Recruited  to  EREs  of  ERa  Target  Genes  Downstream  of  the 
H3R17me2a  Mark.  We  have  shown  that  CARM1  knockdown 
affects  H3R17me2a  and  occupancy  of  PAFlc  at  the  EREs  of 
ERa  target  genes,  indicating  that  PAFlc  binding  is  a  down¬ 
stream  event  of  CARM1  and  H3R17me2a.  The  MCF7-tet-on- 
shCtr9  cell  line  allows  us  to  further  dissect  the  consequence  of 
PAFlc’s  association  with  H3R17me2a.  We  first  comprehensively 
examined  the  association  of  PAFlc  along  the  pS2  gene  during 
E2-stimulated  transcriptional  activation.  The  distribution  pat¬ 
terns  of  Cdc73,  Pafl,  RNAPII,  C-terminal  domain  S2-phos- 
phorylation  (CTD  S2-P),  and  H3R17me2a  on  the  pS2  gene  were 
determined  by  ChIP  assays  using  primers  amplifying  different 
regions  of  the  pS2  gene:  -2,844  (upstream  of  transcription  start 
site),  -192  (promoter  ERE),  +879  (upstream  coding  region), 
+  1,846  (downstream  coding  region).  CTD  S2-P  was  primarily 
found  at  the  elongation  region.  In  contrast,  the  H3R17me2a 
mark  was  only  detected  at  the  ERE  but  not  the  elongation  re¬ 
gion.  Consistent  with  the  established  role  of  PAFlc  as  a  RNAPII- 
associated  elongation  factor,  Cdc73  and  Pafl  were  detected  at  both 
the  promoter  and  coding  regions  coinciding  with  total  RNAPII 
binding  patterns  (Fig.  5/4).  These  data  suggest  that  PAFlc  is  in¬ 
volved  in  transcriptional  elongation  of  ERa  target  genes. 

To  delineate  if  H3R17me2a  could  serve  as  an  upstream  sig¬ 
nal  for  recruitment  of  PAFlc,  we  measured  H3R17me2a  and 
H3K4me3  levels  at  the  ERE  of  pS2  upon  Ctr9  knockdown  (Fig.  5 
B  and  C).  PAFlc  was  previously  shown  to  recruit  Setl,  a  histone 
lysine  methyltransferase  responsible  for  H3K4me3,  and  thus 
H3K4me3  is  downstream  of  PAFlc  and  serves  as  a  positive  con¬ 
trol  to  show  the  effects  of  Ctr9  knockdown  on  transcription. 
MCF7-tet-on-shCtr9  cells  were  pretreated  with  or  without  Dox  for 
5  d,  then  treated  with  20  nM  E2  or  vehicle  for  45  min  before  ChIP 
analyses  followed  by  qRT-PCR.  Fig.  5 B  shows  that  Ctr9  knock¬ 
down  indeed  reduced  H3K4me3  at  the  proximal  ERE  of  pS2, 
consistent  with  a  previous  report  (20).  In  contrast,  H3R17me2a 
level  was  not  affected  (Fig.  5C).  Although  H3R17me2a  level  at 
the  ERE  of  pS2  was  unchanged,  pS2  transcription  was  inhibited 
by  Ctr9  knockdown  (Fig.  4C).  This  result  suggests  that  although 
H3R17me2a  is  typically  defined  as  a  mark  of  transcriptional 
activation,  this  mark  alone  is  insufficient  to  activate  transcription 
without  PAFlc.  PAFlc  was  previously  reported  to  mediate 
H3K4  trimethylation  by  recruitment  of  the  mixed  lineage  leu¬ 
kemia  (MLL)  complex  via  direct  interaction  with  Ash2L,  a 
common  subunit  in  MLL-1  and  -2  complexes  (21,  22).  We  thus 
speculate  that  H3R17me2a  recruits  PAFlc,  which  mediates 
downstream  H3K4  trimethylation  by  interaction  with  the  MLL 
complex  in  ER-mediated  gene  activation.  Inducible  knockdown 
of  Ash2L  cells  were  constructed  by  infecting  MCF7  cells  with 
pTRIPZ-Ash2LshRNA  (Fig.  5 D).  Western  blots  results  showed 
efficient  knockdown  of  Ash2L  by  Dox  treatment  of  MCF7-tet- 
on-shAsh2L  cells  for  6  d  (Fig.  5/)).  As  expected,  ChIP  assays 
showed  that  knockdown  of  Ash2L  significantly  decreased  the 
level  of  H3K4me3  at  the  ERE  of  pS2  (Fig.  5 E),  whereas  the  level 
of  H3R17me2a  mark  and  occupancy  of  Cdc73  at  ERE  of  pS2 
were  not  significantly  affected  (Fig.  5  F  and  G ).  These  results 
suggest  that  H3K4me3  is  downstream  event  of  H3R17me2a  and 
PAFlc  recruitment  in  ER-mediated  gene  activation. 

Given  that  PAFlc  and  H3R17me2a  occupancies  at  EREs 
were  inhibited  by  CARM1  knockdown  (Fig.  3  D  and  E),  whereas 
Ctr9  knockdown  did  not  affect  H3R17me2a  (Fig.  5C),  our  data 
support  the  sequential  transcriptional  activation  model  (Fig.  S7) 
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Fig.  5.  Knockdown  of  Ctr9  did  not  alter  the  H3R17me2a  level  at  the  ERE  of 
pS2;  however,  H3K4me3  level  was  reduced.  ( A )  ChIP  assays  showed  the  oc¬ 
cupancies  of  Cdc73,  Pafl,  total  RNAPII,  and  CTD  S2-phosphorylation  at  distal 
and  proximal  promoter  and  encoding  regions  of  the  pS2  gene.  ( B  and  C) 
Quantitative  PCR  of  ChIP  products  showed  that  knockdown  Ctr9  did  not 
affect  H3R17me2a  level  at  the  ERE  region  of  the  pS2  gene  (C);  however, 
H3K4me3  level  at  the  ERE  was  significantly  decreased  ( B ,  *P  <  0.02).  (D) 
Ash2L  is  effectively  knocked  down  after  6-d  Dox  treatment  in  MCF7-tet-on- 
shAsh2L  cells.  ( E-G )  ChIP  quantitative  PCR  results  showed  knockdown  of 
Ash2L  decreased  H3K4me3  level  at  the  ERE  of  pS2  (£,  *P  <  0.002),  but  had 
negligible  effect  on  H3R17me2a  level  ( F )  and  recruitment  of  Cdc73  (G).  Error 
bars  represent  SD  (n  =  3). 


in  which  CARM1  deposits  H3R17me2a  at  the  ERE,  leading  to 
the  recruitment  of  PAFlc,  which  results  in  H3K4me3  deposition 
and  transcriptional  elongation  via  RNAPII. 

Discussion 

CARM1  and  H3R17me2a  are  linked  to  transcription  activation, 
yet  little  is  known  about  how  H3R17me2a  signals  transcription 
activation.  Using  an  unbiased  biochemical  approach,  here  we 
report  the  identification  of  PAFlc  as  a  unique  transcription 
mediator  of  this  active  histone  mark.  We  showed  that  PAFlc  is 
an  effector  complex  that  binds  to  H3R17me2a  peptide  and  his¬ 
tone  octamers  methylated  by  CARM1  (Fig.  1).  The  specificity  of 
PAFlc  was  demonstrated  by  its  inability  to  bind  to  H4R3me2a 


containing  histone  octamers  modified  by  PRMT1  (Fig.  IE)  (23). 
Nonetheless,  our  in  vivo  ChIP  assays  demonstrated  that  PAFlc 
was  recruited  to  EREs  of  CARMl-regulated  ERa  target  genes 
(Fig.  2),  and  the  occupancy  of  PAFlc  on  the  chromatin  was  de¬ 
pendent  on  CARM1  and  the  H3R17me2a  mark  (Fig.  3).  In  con¬ 
trast,  knockdown  of  PAFlc  subunits  did  not  alter  the  occupancy  of 
CARM1  and  H3R17me2a  at  EREs,  whereas  the  known  down¬ 
stream  H3K4me3  mark  was  attenuated  (Fig.  5 B).  These  results 
strongly  suggest  that  PAFlc  occupancy  at  EREs  is  mediated  by 
H3R17me2a,  which  functions  as  an  activation  mark  because  of  its 
ability  to  stimulate  H3K4me3  through  PAFlc  recruitment. 

PAFlc  is  known  to  participate  in  multiple  steps  of  transcrip¬ 
tion,  particularly  in  controlling  histone  methylation  in  later 
steps  of  transcriptional  regulation.  PAFlc  is  required  for  both 
H3K4me3  and  H3K36me3  during  transcription  elongation  via 
interaction  with  the  histone  H2B  ubiquitination  complex,  the 
Setl  methyl  transferase-containing  COMPASS  complex,  and  the 
Set  2  methylase  (15,  20).  Despite  the  established  roles  of  PAFlc 
in  directing  histone  modifications,  it  remains  unclear  how  PAFlc 
is  recruited  to  its  target  genes.  In  yeast,  genes  encoding  Pafl  and 
Cdc73  are  not  essential,  and  their  loss  results  in  change  of  only 
a  small  subset  of  transcripts  (14,  24).  Similarly,  only  several 
hundred  genes  are  affected,  either  positively  or  negatively,  by 
the  loss  of  PAFlc  in  Cdc73-depleted  HeLa  cells  (25).  Thus,  the 
target-gene  spectrum  of  metazoan  PAFlc  remains  to  be  deter¬ 
mined,  and  the  roles  of  PAFlc  in  gene  expression  are  speculated, 
at  least  in  part,  to  be  determined  by  the  mechanism  of  PAFlc 
recruitment  to  specific  promoters.  Studies  suggest  that  some 
transcription  factors  directly  recruit  PAFlc  to  their  target  genes. 
For  example,  Cdc73  directly  binds  p-catenin  in  both  Drosophila 
and  mammalian  cells  to  activate  transcription  of  genes  in  the 
Wnt  pathway  (26).  Recently,  the  prototypical  transcriptional 
activator  GAL4-VP16  was  shown  to  directly  bind  Pafl  and  re¬ 
cruit  PAFlc  to  the  DNA  template  (27).  Our  finding  that  PAFlc 
regulates  ERa  target  genes  expands  the  role  of  PAFlc  in  regu¬ 
lating  gene  expression.  Furthermore,  we  show  that  the  recruitment 
of  PAFlc  to  EREs  depends  on  CARM1  and  H3R17me2a,  con¬ 
stituting  a  unique  mechanism  for  PAFlc’s  recruitment  to  target 
genes.  Moreover,  the  newly  defined  role  of  PAFlc  as  an  arginine 
methyl-histone  effector  distinguishes  PAFlc  from  its  established 
role  in  controlling  downstream  histone  ubiquitination  and  lysine 
methylation. 

Although  we  demonstrated  the  functional  connection  between 
PAFlc  and  CARMl-mediated  H3R17me2a  in  estrogen-stimu¬ 
lated  transcription,  PAFlc  may  activate  other  transcription  fac¬ 
tors  regulated  by  CARM1  and  functionally  overlap  with  CARM1 
in  other  biological  processes.  In  support  of  this  idea,  PAFlc  and 
CARM1  are  known  to  regulate  common  target  genes.  One  such 
gene,  p-catenin,  is  a  gene  coregulated  by  PAFlc  (26)  and  CARM1 
(28).  Additionally,  both  CARM1  and  the  PAFlc  subunits  are 
essential  for  mouse  embryonic  development  and  maintenance  of 
ES  cell  pluripotency.  CARMl-null  mice  die  at  birth  and  are 
smaller  than  their  wild-type  littermates  (3).  CARM1  and  H3R17/ 
R26me2  are  associated  with  Oct4  and  Sox2  promoters  in  ES  cells; 
moreover,  CARM1  is  required  for  maintaining  pluripotency  of 
ES  cells  (29).  Coincidently,  Cdc73  null  mice  are  embryonic-lethal 
(30).  Microarray  analysis  of  Cdc73  knockout  mouse  embryonic 
fibroblast  cells  demonstrated  that  PAFlc  directly  regulates  many 
essential  genes  involved  in  cell  growth  and  survival  (30).  More¬ 
over,  a  recent  genome-scale  RNAi  screen  identified  PAFlc  as 
a  regulator  of  Oct4  expression,  demonstrating  an  important  role 
for  PAFlc  in  maintaining  ES  cell  identity  (20).  Finally,  both 
PAFlc  and  CARM1  are  involved  in  transcriptional  regulation 
at  multiple  levels.  PAFlc  functions  in  transcription  initiation, 
elongation  (15),  mRNA  3'  end  cleavage  (25),  polyadenylation  of 
mRNA  precursors  (27),  and  small  nucleolar  RNA  (snoRNA)  3' 
end  formation  (31).  In  addition  to  transcriptional  activation, 
CARM1  methylates  splicing  factors  (32),  modulates  recognition 
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of  5'  splicing  sites  (33),  and  regulates  the  coupling  of  transcription 
and  splicing  (34).  Furthermore,  a  recent  report  showed  that 
RNAPII  carboxyl-terminal  domain  is  methylated  by  CARM1 
(35).  Mutating  the  CARM1  methylation  site  in  RNAPII  carboxyl- 
terminal  domain  affects  snoRNA  expression  similar  to  the  effect 
observed  in  CARM1-'-  mouse  embryonic  fibroblast  cells  (35). 
Coincidently,  3'  end  formation  of  mRNAs  and  snoRNAs  are 
dependent  on  PAFlc  (31).  Collectively,  PAFlc  and  CARM1 
share  overlapping  genes  and  cellular  processes,  although  the 
processes  that  require  CARM1  methyltransferase  activity  and 
H3R17me2a  for  PAFlc  function  await  elucidation. 

Arginine  methylation  of  histone  tails  correlates  with  either 
transcriptional  activation  or  repression.  The  specific  marks,  as 
part  of  histone  codes,  are  thought  to  promote  or  prevent  the 
docking  of  key  transcriptional  effector  molecules  (36,  37).  PAFlc 
subunits  are  nuclear  proteins  directly  associated  with  RNAPII 
and  can  regulate  multiple  steps  of  transcription.  Recently, 
TDRD3  was  found  to  bind  both  the  H3R17me2a  and  H4R3me2a 
marks  (10).  It  is  plausible  that  multiple  effector  proteins  exist  for 
the  H3R17me2a  mark  that  direct  this  mark  to  divergent  bi¬ 
ological  processes.  Future  work  should  determine  if  these  two 
H3R17me2a  effectors  directly  interact  and  whether  they  regulate 
different  CARMl-coupled  biological  processes  in  cell-type  and 
tissue-specific  manners.  In  conclusion,  we  identified  PAFlc  as  an 
effector  complex  for  H3R17me2a  and  linked  its  function  to  the 
regulation  of  estrogen-responsive  genes.  This  study  presents  a 
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unique  mechanism  by  which  the  CARMl-catalyzed  H3R17me2a 
mark  is  linked  to  transcriptional  activation  via  recruitment  of 
PAFlc  (Fig.  S7). 

Materials  and  Methods 

Cell  Culture  and  Estrogen  Treatment.  MCF7  and  HEK293  cells  were  maintained 
in  DMEM  supplemented  with  10%  (vol/vol)  FBS.  Three  days  before  E2 
treatment,  the  cells  were  transferred  to  phenol  red-free  DMEM  containing 
6xSFBS  (six-time  stripped  FBS).  The  working  concentration  for  E2  was  20  nM. 

Peptide  Pull-Down  Assays.  Peptide  pull-down  was  performed  according  to 
previously  published  methods  (12).  For  each  assay,  biotinylated  peptide  (1 
pig)  were  immobilized  on  200  [ig  streptavidin  coated  DynaBeads  (Invitrogen). 
HeLa  nuclear  extract  or  purified  recombinant  proteins  were  incubated  with 
the  magnetic  beads  immobilized  with  indicated  peptide  for  2  h  at  4  °C.  The 
peptide  interaction  protein  was  detected  by  Western  blot  or  silver  staining. 

Description  on  cell  lines,  antibodies,  peptide  pull-down,  qRT-PCR,  ChIP 
assay,  transfection  and  in  vitro  methylation  of  histone  octamer  are  available 
in  SI  Materials  and  Methods.  All  of  the  primers  and  oligo  DNA  are  listed  in 
Table  SI. 
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